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FOREWORD 


The balloon has a long history of use as a platform for upper atmospheric 
research. In fact for many years it was the only vehicle capable of obtaining 
in situ data on the lower stratosphere. The development of the various poly- 
mers which could be extruded into uniform thin films dramatically increased the 
capability of the balloon as a lower stratospheric research vehicle. The 
materials not only made it possible to fabricate much larger balloons at low 
cost but also to construct the units so that rather than bursting they remained 
at float altitude for extended periods of time. Initially the units were 
designed with modest capabilities (100 kgs to 30 km). The existence of this 
capability increased the number of investigators and the range of research being 
performed with them. This, in turn, led to interest in increasing the altitude 
and load capability of the vehicle. The size of the balloons available 
increased and thickness of the material required to fabricate the balloon 
decreased rather rapidly so it became possible to get light (<1000 kgs) payloads 
to 35 km. Heavy payloads (>1000 kgs) were a problem and most heavy payloads 
were flown on balloons constructed of a mylar scrim. Since these were very 
expensive their use was limited. Another problem occurred as the larger units 
were used at low latitudes and this was the failure of the units at the tropo- 
pause . It was recognized that the failure was due to the material becoming 
brittle at the cold temperatures associated with the low latitude tropopause. 
This led to the development of a material designed with better cold brittleness 
characteristics for fabricating balloons. With the additional development of 
load tapes and caps to handle heavy loads, both problems were overcome, and the 
balloon has become a reliable heavy load platform. It is interesting to note 
the recent records achieved with large balloons. These include carrying a 4536 
kg payload to 34 km and a 232 kg payload to 51.8 km. This capability has been 
used by many investigators to obtain data on the middle atmosphere. 

In this volume we discuss some of the techniques employed by investigators 
using balloons to obtain data on the properties of the middle atmosphere. The 
balloon provides an excellent platform for obtaining such data and the number of 
investigators using it is large. As a result it is difficult for a volume such 
as this to be complete and several techniques are not included. This lack of 
completeness is also partially because manuscripts from several authors were not 
received in time for inclusion in this volume. The editor apologizes for the 
omissions and emphasizes that their omission does not imply that the techniques 
are unimportant. 

The first stratospheric parameter to be measured with balloon-borne instru- 
ments was temperature. Temperature continues to be an important parameter and 
most of the data concerning the temperature and dynamics of the lower strato- 
sphere continue to be obtained by means of balloon-borne sensors. This is due 
to the development of the small radiosondes which can be flown to 30 km on small 
inexpensive balloons. The importance of measurements of the variable minor 
species (0. and HjO) has long been recognized. A great deal of effort has 
gone into developing instruments which could be used on small balloons to 
measure these constituents. These are discussed in detail in the first chapter. 

The question of whether man' s activities could affect the composition of 
the lower stratosphere (particularly the ozone layer), which arose during the 
controversy over the building of the Supersonic Transport aircraft, increased 
the interest in the measurement of stratospheric composition. This interest 
carried over into the study of other possible sources of pollution including the 
possible effects of chlorof luormethanes on the ozone layer. This in turn added 
to the list of compounds for which measurements of the stratospheric concentra- 
tion were badly needed. Many of these compounds were predicted to be present in 
the ppt to ppb range and instruments were required that could measure these 
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compounds from a balloon platform. Several different techniques have been used 
to obtain such information and rather than arranging the chapters according to 
compounds, they have been arranged according to the measuring technique. In 
chapters 2 through 4 various remote sensing techniques used to obtain data on 
atmospheric composition are described. These are arranged according to the 
wavelength regions used to make the measurement. 

Several in situ techniques have been developed for use on balloon plat- 
forms. In addition to the measurements of neutral constituents, a number of 
groups have developed instrumentation required for measurement of the ion 
composition from balloon platforms. In many respects, the balloon measurement 
has proven to be much more difficult than measurements at the higher altitudes. 
These are discussed in chapter 5. 

Interest in the stratospheric aerosol seems to wax and wane with volcanic 
activity and the last few years has seen a tremendous increase in this interest. 
The balloon has proven to be an excellent vehicle for studying both the physical 
and chemical properties of the aerosol. The final chapter in this volume 
reviews the balloon measurements in this area. 


D. G. Murcray 
Editor 
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1. THE MEASUREMENT OF TEMPERATURE, OZONE AND WATER VAPOUR 


T. McElroy 

Atmospheric Environment Services 
Downsview, Ontario, Canada M3H 574 


INTRODUCTION 

In 1904 Teisserenc de Bort published the results of 581 unmanned balloon 
ascents during which recording instruments measured temperature and pressure. 

As a number of these flights reached 14 kilometers or more, de Bort was able to 
observe the change from the normal tropospheric adiabatic lapse rate to the near 
isothermal conditions of the stratosphere. In fact, de Bort named the strato- 
sphere and troposphere to indicate the difference in static stability which is 
typical of these tvo regions of the atmosphere. The determination of the 
vertical, temperature and pressure profiles of the atmosphere is an important 
task of every stratospheric research flight. 

Ozone was clearly identified in the atmosphere when FABRY and BUISSON 
(1921) found characteristic absorption features in spectrographic measurements 
of the sun. DOBSON (1930) reported the development of a spectrophotometer 
designed specifically for the ground-based measurement of ozone. The instrument 
measures the differential absorption of light, by o zone, at two wavelengths in 
the near ultraviolet. The observed absorption is converted to an amount of 
ozone based on laboratory measurements of the absorption coefficients of ozone. 
The Dobson spectrophotometer still serves as the primary standard of measurement 
for the amount of ozone in the atmosphere. 

Because of the importance of knowing the temperature and humidity structure 
in the troposphere for the purpose of forecasting the weather, temperature and 
humidity measurements as a function of pressure have been made operationally for 
many decades. It is only recently, however, that studies of the structure of 
the upper atmosphere has become of general interest. 

During the second world war the British Ministry of Defense began to 
measure the humidity structure at high altitude (20,000 - 40,000 feet). DOBSON 
(1946) and BREWER (1949) undertook the measurement of the very low values of 
stratospheric humidity using an airborne, manually operated frost point hygro- 
meter. The results of a recent program of water vapour measurements which grew 
out of the wartime effort are reported by CLULEY (1978). 

Since these early observations, many methods have been employed to make 
measurements of these three important atmospheric parameters. Some of these 
were developed, or conducted through the use of the stratospheric research 
balloon, or high altitude sounding balloon. A sampling of these measurement 
techniques will be presented in more detail in the pages to follow. 

PRESSURE, TEMPERATURE AND HUMIDITY 

In order for an operational weather forecast to be prepared, it is 
necessary to have a continually updated knowledge of the curent state of the 
atmosphere. This current "picture" or analysis of the atmosphere includes as 
one of its fundamental elements a measure, over the surface of the earth, of the 
height of standard pressure levels. The determination of these constant 
pressure surfaces is accomplished by the release of radiosonde sounding 
instruments to determine the pressure, temperature and humidity of the 
atmosphere. The interrelation of these three atmospheric variables with height 
and air density is governed by two important physical relations. One of these 
is the perfect gas law, which may be expressed as: 




p = pKT (1) 

3 

where p is pressure (in Pa), p density (molecules/m ), T temperature (°K), 
and K is Boltzmann's constant (1.3806 E -23 J/ °K/molecule). 

This relation will apply instantaneously at each point in the atmosphere 
over the range of pressures and temperatures encountered there. The second 
relation is the equation of hydrostatic balance for a parcel of air in the 
atmosphere, which expresses the fact that the pressure difference across the 
parcel must support its weight due to the local gravitational attraction. The 
relation may be written as: 3 

-J 

dp = -gpMdz (2) 

where dp is pressure difference across parcel, dz the thickness of the parcel, 
g the acceleration due to gravity (9.8 m/s ), and H the mass of one air mole- 
cule (4.81 E -26 kg). 

Combining equation 1 and 2 to eliminate density yields the usual form of 
the equation of hydrostatic balance: 


dp gM 

P KT 


dz 


(3) 


The logarithmic change of atmospheric pressure with height is immediately 
obvious with equation (3) integrating to: 


In 


PL 

PU 



(4) 


where subscripts L,U represent the lower and upper limits in altitude of the 
range of integration. 


The variables g,T are left inside the integral since they may in general 
both be functions of height, that is g(z), T(z). One useful simplification is 
to take g as a standard, surface value g , and remove it from inside the 
integral. If this simplification is made, the resulting heights, which will not 
quite agree with the actual physical height of a given pressure surface, will be 
expressed in geopotential meters. Since g(z) falls as the inverse square of the 
distance from the center of the earth, and heights within the stratosphere are 
small compared with the earth's radius, the dif ference ' is not great. It may be 
important, however, if measurements from radar heights are to be compared to 
measurements which have been plotted as a function of geopotential height. The 
relation between geopotential and geometric height is tabulated in standard 
atmospheres for particular atmospheric temperature profiles, (e.g.: US Standard 
Atmosphere; ICAO Standard Atmosphere; COSPAR International Reference 
Atmosphere. ) 

If the temperature is also taken to be constant, 'which can be true for 
reasonable distances in the stratosphere, equation 4 becomes particularly 
simple: 




( 5 ) 


or 
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KT 

where H = — is the scale height of the atmosphere (8.3 km for 300°K). 


The scale height of the atmosphere is the pressure e - folding distance 
in the vertical. It is an approximate measure of the expected vertical 
resolution of several remote sensing techniques. Equation 5 also has a simple 
solution for a linear rate of change of temperature with height: 
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where b is (gM)/(Ka). a is the temperature lapse rate with respect to height. 

When an atmospheric sounding is made, the variables actually measured are 
pressure, temperature, and humidity. In order to determine the height at which 
a pressure and temperature measurement was made, the whole sounding from the 
surface pressure up to the pressure of the observation must be numerically inte- 
grated according to equation 3 to yield the height corresponding to each 
pressure. The perfect gas law, equation 1, may be used to calculate the local 
air density at any point from the corresponding pressure and temperature. 

A useful re-expression of the6e equations may be developed by noting that 
the pressure difference across a macroscopic layer: 

*P " P L " Pu (7) 

must be sufficient to support the weight of molecules between those pressure 
levels. This relation yields an accurate mean number density for a particular 
layer since: 

Ap = gp MAz (8) 

where Az^Z^ - is the layer thickness and p is the mean number density 
in the layer. 


The use of this relation prevents the "loss" of molecules which will result 
from the imprecision of numerically integrating the density function through the 
layer. 

Since temperature errors of a few degrees K can result in a height error of 
more than 1000 meters at a float height of 30 km, it is important that accurate 
pressure and temperature data be available for a particular balloon flight 
analysis. It may not be acceptable, in relation to the accuracy of other 
measurements made, to rely on a standard atmosphere in the reduction of data. 
This is particularly so in the case of solar occultation data, and long path 
length data at near horizontal angles of view. For example, LOUISNARD (1983) 
calculates that the air mass traversed by the incident solar beam may be in 
error by as much as 5% from an error in calculated balloon height of 0.5 km when 
a layer at 20 km is observed from a balloon height of 40 km. 


The measurement of temperature then has two important contributions to the 
analysis of a research balloon flight. One of these is the calculation of the 
air density as a function of height for the ambient atmosphere, and the other is 
the provision of local temperature measurements for those instruments, 
principally emission radiometers, which require it for the evaluation of their 
data. 
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THE MEASUREMENT OF TEMPERATURE 

A considerable number of atmospheric sounding devices have been built over 
the years. They all have incorporated some relatively simple device for con- 
verting the local air temperature into an electrical signal, which is then 
transmitted via a radio link to the ground. While some unusual devices such as 
magnetic flux modulating thermometers have been used, (British Kew radiosonde) 
the majority of radiosondes have employed a simple thermistor, which is a 
temperature variable resistance. 

Since the VIZ radiosonde has been reported as being used in 20 countries 
around the world (RICHNER, 1981), and it is also the most." likely choice to be 
used in support of a large balloon payload, it is the only device which will be 
described here. Similar comments to those made here would be, in a large part, 
appropriate in the case of most radiosondes currently in use around the world. 
The global temperature data obtained by radiosondes forms the standard of 
evaluation of performance for other techniques, such as temperatures deduced 
from satellite atmospheric radiance measurements. 

The standard VIZ radiosonde, of United States manufacture, is a very simple 
instrument. It has an aneroid driven pressure contact switch which selects 
between the output of the humidity sensor, the thermistor output, and internal 
reference signals. Each time a pressure contact changes from one signal source 
to another it indicates that the instrument has passed through the corresponding 
pressure level. Each baroswitch is calibrated after construction and the 
calibration data are supplied along with the instrument. 

Many details must be properly looked after if the instrument is to provide 
optimum performance. To list just a few of these: (1) -The temperature and 
humidity sensors must be properly checked on the ground with a baseline testing 
unit; (2) The pressure aneroid must be adjusted to reflect the surface pressure 
at launch; (3) The humidity sensor must be handled carefully to avoid contamina- 
tion; (4) The thermistor must be clean and properly rigged to avoid a large 
radiation correction. In short the radiosonde must be carefully used by know- 
ledgeable operators, 

A number of studies of the performance of the radiosonde system have been 
conducted. The most useful of these are the two-sonde balloon tests, since 
these demonstrate random errors in the sondes of the absence of a highly 
accurate knowledge of the atmospheric state when the tests are done. Notable 
examples of this type of inter comparison can be found in the work of RICHNER 
(1981), LENHARD (1970, 1973), PHILLIPS (1983), SCHMIDLIN (1982) and ANTIKAINEN 
(1983), Generally the results suggest a random error of approximately 1 to 2 mb 
in pressure and 0.4 to 0.5 °K in temperature. 

The most recent work, that of Phillips, suggests an expected error of 
approximately 1 mb in pressure. With the exception of Schmidlin, no attempt is 
made to discuss absolute accuracy in the determination ."of atmospheric 
temperature and pressure. Rather the random differences between particular 
instruments is measured, and the instrumental bias is taken to be zero. 

An important point is made by Richner, who suggests a radiational heating 
correction of 1 - 2°K should be applied to the temperatures measuring during the 
day by the VIZ radiosonde. It should be noted that the correction is pressure 
dependent, and is quite important above an altitude of approximately 20 km. 

This effect is also noted by PHILLIPS (1983) who has used the correction cited 
by Richner in compiling his tables of relative performance among a group of 
radiosondes of diverse manufacture. 
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Schmidlin* s paper asserts that radar is the only solution to the height 
accuracy problem in the use of a radiosonde. Be suggests that the error of 
±2 mb quoted from the results of the study will cause an uncertainty of close 
to 1 km in geopotential height calculated using radiosonde ascent data. As 
with the other authors mentioned the study reported showed rather small random 
temperature effects compared to the influence of pressure errors. It is not 
clear from the paper, however, whether a radiation correction was applied to 
the measured temperatures, and whether that would have influenced the final 
conclusion. 

In addition, Schmidlin' s paper does not include the analysis of another 
option open to the stratospheric investigator. When using a radiosonde to the 
maximum limit involving altitudes near 40 km, the inclusion of a hypsometer is a 
sensible alternative to the acquisition of a radar set. Performance data for 
the hypsometer sondes reported by LENHARD (1970, 1973) clearly indicates the 
potential for the radiosonde to provide useful pressure data, and geopotential 
heights to more than 35 km. 

It would appear that the radiosonde instrument, properly handled, with 
either a hypsometer or radar height recording, will provide adequate vertical 
temperature, pressure, and density information for the support of research 
balloons. A factor which influences the accuracy of a constituent profile is 
the quality of the balloon position data. If the radiosonde is supplying direct 
support to an in situ sensor, then there is no dislocation between the position 
of the experiment and that of the radiosonde. If, however, the sonde is 
supporting a payload with solar absorption instrumentation, then the actual 
position of the payload in latitude and longitude is required. In addition, if 
the height of the payload is needed and the radar is not available, small errors 

may develop in the height determined by the pressure data, since the payload may 

not be contained in the same region of the atmosphere through which the radio- 
sonde climbed to determine the temperature-pressure profile. Errors of this type 
are harder to quantify, although some discussion of closely spaced balloon 
launches is presented by LENHARD (1970, 1973). 

A useful guide to the general reliability of the atmospheric profile 
estimate for a given balloon flight may be made by the release of a series of 
radiosondes. If several ascents can be made within ±6 hours of the period of 
interest, the data set may be considered as a whole, and a better assessment of 

the state of the atmosphere at a particular time and place may be made. In 

addition the magnitude of the random component in the error in determining the 
atmospheric variables may be estimated by the judicious comparison of profile 
differences observed in the time series. 

For example, a series of atmospheric soundings taken with VIZ radiosondes 
in support of the Canadian Stratoprobe program were published by BAIN (1976). 

Radiosondes are not the only method by which the atmospheric temperature 
profile may be measured using a balloon. The emission of radiation by the 
molecules which make up the atmosphere may be detected by a spectro-radiometer 
and the temperature profile derived therefrom. SHAW (1970) and McCLATCHEY 
(1969) report the successful determination of a vertical temperature profile in 
this way. 

WATER VAPOUR 

ELLSAESSER (1983), in a recent review article, has observed that in the 20 
years or so which have elapsed since GUTNICK's (1961) review of water vapour 
measurenents, the variability of the measurements has not significantly 
decreased. Both the variability among the various measurement techniques and 
the variability from occasion to occasion for the use of one technique, have 
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provided such a range of values as to leave some fundamental questions 
unanswered. 

NEWELL (1981) also has made a similar point, although has gone on to assume 
that the largest bodies of data available, the flights of MASTENBROOK (1974a) 
and the Meteorological Research Flights of the British Meteorological Office are 
substantially correct. This is to assume that the mixing ratio values at the 
extreme low end of the scale are closest to the correct values. The extreme 
dryness which this implies for the stratosphere is then extremely difficult to 
explain from a theoretical standpoint (NEWELL, 1981; ELLSAESSER, 1983; 

DANIELSON, 1982). 

Another unsolved problem in the stratosphere is the reconciliation of the 
mixing ratio behaviour with height, with the expected HjO source from methane 
oxidation in the stratosphere. ELLSAESSER (1983) suggests that methane 
oxidation should provide a source of water vapour of up to 2 ppmm at heights 
of 30-40 km. The Mastenbrook results are substantially constant with height 
in the stratosphere, but some Murcray results show an increase with height which 
is greater than that expected if it were due entirely to^the oxidation of 
methane. It is instructive to consider the methods and instruments with which 
these conflicting data have been gathered over the years. 

In the following short descriptions of instruments and techniques, no 
particular attempt has been made to create a chronological review. Rather the 
techniques selected for elaboration have been chosen according to the following 
rather loose criteria. The selected techniques are usually the best known and 
most recently used of their type. Furthermore any researcher whose results are, 
in the light of current knowledge, rather extreme when considered as a part of 
the extant data set, are not used as examples of measurement technique, even 
though their method has been successfully employed since. 

The current concept of the stratosphere as a dry region of the atmosphere 
with mixing ratios on the order of 1 to 2 ppm by mass, owes its existence to the 
work of Dobson and Brewer. Many measurement flights on aircraft using manually 
operated frost point hygrometers were conducted both during and after the second 
world war by the British Meteorological Office. Data from these flights are 
reported by DOBSON (1946) and by BREWER (1949), and for the later Meteorological 
Research Flights by CLULEY (1978). While these flights produced the first large 
data set which accurately represented the humidity conditions in the strato- 
sphere, they were limited, even with the use of higher flying aircraft in later 
years, to the lower stratosphere. 

An obvious solution to the altitude limitation was the use of stratospheric 
balloons for sampling the water vapour content of the stratosphere. The many 
and varied techniques for the measurement of water vapour form the basis of the 
following section of this chapter. 

THE FROST POINT HYGROMETER 

The principle of the frost point hygrometer has been employed by a number 
of investigators over the years to make measurements in the stratosphere. 

Indeed, as mentioned in the introduction to this section, the first measurements 
indicating the existence of a "dry" stratosphere were made on board aircraft by 
DOBSON (1946), using a manually operated frost point hygrometer. The most 
extensive set of balloon-based observations of water vapour using frost point 
hygrometers has been contributed to the literature by MASTENBROOK (1961, 1965, 
1968, 1971, 1974a, 1974b, 1979). 

The Mastenbrook frost point hygrometer was described in a review delivered 
by MASTQIBROOK (1974b) at the Australian conference on the structure and 



composition of the stratosphere, and in more detail in MASTQIBROOK and DINGER 
(1961). 


7 


The instrument consists of a condensing mirror surface on which frost is 
allowed to form. Air is drawn into the measurement cavity of the instrument by 
a blower and passed over the mirror surface. The mirror is connected to a heat 
sink which consists of freely boiling Freon 13, and induction heated, so that 
its surface may be maintained at any tsnperature above the boiling point of the 
Freon 13 by adjusting the heating rate. 

In order to detect the formation of frost, an optical sensing element views 
light scattered from frost which has condensed on the cold mirror. The optical 
signal is detected and used to control a servo loop through the induction 
heater. The system may be balanced with a small, equilibrium amount of frost on 
the reflecting surface. The temperature of the surface is detected with a 
thermistor, and the measured temperature is taken to be the frost point 
temperature of the air sampled. 

In order to avoid contamination problems from the presence of entrained 
water in the vicinity of the balloon and payload, measurements of stratospheric 
water vapour amounts are made only in the descent phase of the balloon flights. 
In addition the mirror is temporarily heated at the top of the flight profile to 
ensure good sensitivity to the small vapour pressures observed in the strato- 
sphere. 

MASTENBROOK (1974a) has conducted a large number of balloon ascents since 
the beginning of his measurement program at the National Research Laboratory in 
the early sixties and has discussed them, and their associated error analysis in 
considerable detail. It would seen appropriate, however, to mention briefly 
some of the systematic error sources which he has identified in the course of 
his research. 

The problem of contamination, which has been discussed by other authors 
as well, (ZANDER, 1966) is probably the single biggest difficulty encountered in 
the in situ measurement of the low water vapour mixing ratios observed in the 
stratosphere. Mastenbrook's early work identifies two related, but distinct 
methods by which some of the early ascents were contaminated. The first of 
these is the balloon/payload water vapour cloud, which is carried aloft into the 
stratosphere in the form of water vapour absorbed onto the balloon and payload 
surfaces during the ascent through the wet troposphere. Water may also be 
evolved from the interior of the materials used in the payload construction at 
the lower pressures encountered in the stratosphere. This effect was eliminated 
by the use of a descent measurement so that the vertical velocity of the 
instrument package presented uncontaminated air to the instrument, without 
allowing water vapour from the balloon or instrument package to mix into the air 
sample, and contaminate the measurement. 

The second problan encountered was internal contamination in the hygrometer 
sampling system itself, which was caused by an inadequate airflow through the 
instrument. The inadequate airflow resulted from a malfunctioning sampling 
blower; a problem which was not encountered during the later flights on which 
Mastenbrook placed his greatest emphasis. 

The final conclusions drawn from the series of measurements reported by 
Mastenbrook are: (1) The stratosphere is dry with a mixing ratio of 2.3 ppmm 

in the low stratosphere; and (2) That there is no substantial gradient in water 
vapour mixing ratio between this level and 25-30 km altitude. 

It is tempting to accept the study of Mastenbrook as the final assessment 
of water vapour in the stratosphere for two reasons. First, the low altitude 
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measurements are in reasonable agreement with extmsive sets of data collected 
by the Meteorological Research Flight of the Meteorological Office in Bracknell, 
England. (When a correction for frost point temperature measurement is made 
(OLIVER, 1978).) The aircraft measurements are unlikely to be contaminated with 
water vapour since the air collection is done at very high ventilation rates, 
and the aircraft continuously moves at high speed into new stratospheric air. 
Second, Mastenbrook has achieved these results after grappling extensively with 
the problan of balloon contamination, and has assembled a large enough data set 
that temporal variations are not likely to bias his mean results to a great 
extent. It must be remembered, however, that at least one flight represents a 
reliable measurement of water vapour in the stratosphere substantially higher 
than the average value which Mastenbrook reports. 

It is probably safe to conclude as Mastenbrook has done, that other 
reported frost point hygrometer flights with mixing ratios an order of magnitude 
or more larger (BARCLAY, 1960; BARRETT, 1950) were in fact the result of the 
local contamination. Even more recent flights (BROWN, 1964) have also been the 
victim of contamination. As has been expressed elsewhere (ELLSAESSER, 1983; 
GUTNICK, 1961) it is extremely difficult to sort the "good" results from the 
"bad" with high reliability. 

AN ALUMINUM OXIDE HUMIDITY SENSOR 

The useful range of conventional radiosonde humidity sensors is such that 
they are not of any use in measuring the water vapour content of the strato- 
sphere. For example, the Canadian radiosonde network does not regard results 
obtained in ambient temperatures below -40 °C as useful. In order to provide an 
inexpensive and easily employed detector for stratospheric humidity sensing, it 
would be desirable to have a simple device like that of the standard radiosonde 
hygristor, available for use on small and simple radio sonde- type instrument 
packages. Farametrics Corporation has attempted to develop such a sensor, and 
some results from its use are available in the literature (HOFMANN, 1972; 

CHLECK, 1979). 

CHLECK (1979, 1966) describes the development and use of an aluminum oxide 
water vapour sensor (Parametrics Aquamax II) capable of the measurement of very 
low levels of humidity in the atmosphere. The sensor is a thin film fabrication 
incorporating an integral heater, temperature monitor and humidity sensor. The 
sensor element is a porous dielectric layer formed on a layer of aluminum, which 
is deposited on a silicon substrate. The dielectric layer is overcoated with a 
layer of porous gold which forms one electrode of the sensor. The aluminum 
beneath the dielectric forms the other. The water vapour pressure in the 
vicinity of the sensor alters the impedance measured between the electrodes as 
water is taken up or released from the dielectric layer. The complete sensor 
occupies a 5 mm square area on the thin film substrate, with a 1 mm total 
thickness. While measuring, the temperature of the sensor is stabilized by 
ohmic heating with the integral heater, to which current is supplied by an 
external control circuit. The control circuit monitors the temperature of the 
substrate with the integral diode temperature monitor. 

The temperature of the device is stabilized in order to eliminate the 
temperature corrections which were required in earlier forms of the device 
(CHLECK, 1966). During radiosonde flights the temperature of the sensor is 
telemetered to the ground to verify the correct operation of the heater, and to 
ensure that the temperature remains within acceptable limits. 

The sensor was tested in a series of flights reported by CHLECK (1979) and 
HOFMANN (1972) in which they describe both the calibration procedures used with 
the sensor, and the results obtained in atmospheric tests. The sensors were 
individually calibrated before flying, since the sensor does not have an 
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absolute physical response to water vapour pressure. 

Two calibration procedures were employed. One used the humidification of a 
variable stream of dry, cold nitrogen by a fixed rate of evolution of water 
vapour from a temperature stabilized permeation tube. The total mass of water 
evolved from the tube was determined by the weight loss experienced throughout 
the calibration procedure. Various flow rates of nitrogen were humidified by 
the fixed water vapour evolution rate to yield the varying humidity levels to 
which the sensor was exposed. 

In addition to the permeation tube calibration a more conventional humidity 
source, a bubbler, was used to produce air at various humidity levels. The gas 
thus produced as passed over the sensor, and measured by a frost point hygro- 
meter which was used to measure the humidity down to a frost point temperature 
of -89 °C. 

The observations made with these methods were intercompared to yield an 
empirical relation for the sensor response. This relation may be used to extra- 
polate the sensor response down to a frost point of -100°C. The final accuracy 
of the calibrated sensor is estimated to be x3°C at a frost point of -65°C. 

In addition to the humidity calibration, the response of the sensor to step 
changes in the humidity was determined. As an example the sensor responded to a 
frost point change from -90°C to -75°C in 70 seconds, and to the opposite change 
in 30 seconds. 

The sensor was tested on a series of balloon flights while connected to a 
standard radiosonde by a 7.6 m cable hanging below the radiosonde. Measurements 
were made both on the ascent and descent legs of the flight. A clear difference 
in the profiles, presumably due to water vapour outgassed by the balloon, is 
seen in the profiles. The true response of the detector can be seen in the low 
water vapour amounts recorded in the descent phase of the flight. The measure- 
ments made included some as low as 1 ppmm, and are in good agreement with those 
of other investigators as stated in the paper (CHLECK, 1979). 

CRYOGENIC SAMPLING 

Several investigators have employed cryogenic sampling instruments to 
collect samples of stratospheric air for subsequent analysis in the laboratory. 
An early attempt was made by BARCLAY (1960) but the large amounts of water 
vapour measured (on the order of 37 ppmm) suggests that the measurement was 
contaminated by water carried into the stratosphere by the sampling system. 

More recent papers by POLLOCK (1980) and LUEB (1975) describe a cryogenic 
sampling system for the collection of stratospheric air samples, which will be 
described in some detail here. 

The collection system consists of a set of stainless steel cylinders 70 cm 
long by 3.8 cm diameter designed to collect cryogenical ly cooled samples of 10 
litres STP volume at all flight levels in the stratosphere. The collected 
gasses, which include all atmospheric constituents except molecular hydrogen, 
neon and helium are collected in the cylinders, which are partially immersed in 
liquid neon which maintains the collecting surface of the cylinder at approxi- 
mately 30°K. 

Up to sixteen cylinders are connected together to form the sampling system 
which is illustrated in Figure 1, taken from LUEB (1975). The individual 
cylinders each have motor operated valves which can be opened or closed in 
flight by command. These valves lead out to a common manifold, which in turn 
connects to the sampling inlet tube. The gas samples are taken at a position 
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Figure 1. Schematic diagram of the balloon- 
borne low temperature sampler. Sample 
cylinders are arranged in concentric circles 
of eight cylinders with the purge pump in 
the centre. The sample cylinders are sub- 
merged in liquid neon. Reprinted from LUEB 
(1975). 


6 m below the gondola which is slung 75 to 150 m below the balloon. These 
precautions are intended to reduce the risk of contaminating the samples 
collected. 

The whole system is baked out at 350°C to AQO^C at high vacuum to remove 
water which might be absorbed onto any of the surfaces with which the sample 
gases might come in contact. The stainless steel is exposed at moderate 
pressure to a source of O2 to passivate it and the system re-evacuated. The 
inlet tube is sealed with a rupture diaphragm, which is blown out in flight by a 
high pressure supply of gaseous nitrogen. The rupture of the disc allows an 
extension, which had been sealed inside the collection tube, to add 0.7 meters 
to the length of the collection tube. In addition to extending the distance 
which the collection inlet is below the gondola, it also ensures that the sample 
taken is always at least 0.7 m away from any unprepared surface. 

Before the start of each sample collection the inlet tube and manifold is 
purged by collecting gas in an extra sample cylinder which is provided for that 
purpose . 

The flight profile is arranged so that the sampling system is held at float 
altitude for approximately two hours before sampling is started. After the 
purge has been done, a sample taken at float altitude, to be analyzed for the 
amount of contamination which collection in the balloon environment has 
produced. The payload is then allowed to descend while the main collection 
procedure is carried out. With a slow descent in the early phases of the 
process and a more rapid one later on, all samples may be easily collected 
within a vertical range of 0.5 km each. This allows good height resolution in 





11 


the final distributions of gases which are analyzed. The collection procedure 
is carried out from a float altitude of 35 km down to about 12 km, where the 
balloon is normally cut down. 

The water vapour contamination problem is illustrated by the comparison of 
the mass mixing ratio observed in the float altitude sample and the first 
descent sample. On the flight reported by LUEB (1975) the float sample 
exhibited a mixing ratio of 21 ppmm, while on the first descent sample a mixing 
ratio of only 5.7 ppmm was found. 

When the samples are returned to the laboratory they are in the form of 
high pressure gas samples. The water vapour is separated in a cryogenic trap, 
converted to hydrogen gas, and the volume of the hydrogen is measured. 

It is immediately obvious that the operational procedures and design 
elements described by Pollock and Lueb are thoughtfully aimed at controlling 
those uncertainties which have caused errors in measurements of trace gases by 
other methods. The procedures for handling the equipment in the laboratory 
control the internal contamination problem, while the balloon contamination 
problem is addressed by the use of an appropriate flight profile. Also, the 
large volume of the sample reduces the contamination problem. 

The authors point out a number of significant, advantageous features of 
their collection technique. These include the fact that the method allows the 
simultaneous sampling of a large number of interesting stratospheric trace gases 
(CH^, H 2 , CO, R.O, CFCl^ , CF 2 C1 2 ) in addition to water vapour. 

Because of the large volume of gas collected, the measurement of HDO and HTO 
i6 also possible. 

The technique is expensive, the equipment is complicated and cumbersome, 
and the samples require laboratory analysis after the flight; however, for the 
simultaneous collection of many stratospheric species it is a very effective 
technique . 


ULTRA VIOLET FLUORESCENCE MEASUREMENTS 


BERTAUX (1978) and KLEY (1978, 1979) describe two independently developed 
instruments for the in situ measurement of water vapour in the atmosphere. 
These instruments are similar in operation; and that of Kley will be presented 
here in detail. The method is based on the detection of fluorescent radiation 
from the decay of OH(A l ) to OH(X i + ). The excited form of OH is produced 
in the photodissociation reaction: 


H 2 0 + hv f OH(A 2 E + ) + H 

where the radiation involved is in the vacuum ultraviolet. 


(9) 


In the instrument, shown schematically in Figure 2, taken from KLEY (1978), 
light from lamp L which emits Lyman alpha radiation at wavelength 121.6 nm, is 
incident on an air sample. The fluorescence of OH as a result of the afore- 
mentioned reaction is observed by a photomultiplier detector oriented at right 
angles to the path of the dissociating radiation at a distance L g from the 
lamp. The intensity of the beam at a distance L^ is monitored by an NO cell 
(NOC^) and the intensity of the beam at the point where the fluorescence is 
detected, is monitored by an NO cell (NOC g ) at distance B’ + B" (=L g ) from 
the lamp. 


The presence of two detectors, NOC^ and NOC g , allows the development of 
relations for the water vapour mixing ratio which are independent of the 
absolute light intensity supplied by the source lamp. 
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NOC b 


Figure 2. Schematic of the OH fluorescence 
flight instrument. L - Lyman-alpha lamp; 
NOCa, NOCp - NO ionization cell 
detectors; PM - photomultiplier; M - mirror. 
1*A “ 18.5 cm; Lg = 9.4 cm. The distance 
B ' + B" is equal to Lg. After KLEY (1979). 


In the troposphere where water vapour concentrations are large, the 
instrument provides a measure of the amount of water in an absorption mode. 

The depletion of the beam along a path equivalent to La - L fi is used to 
calculate the number of water vapour molecules in the path of the beam. A 
correction to the observed absorption is made for the influence of oxygen mole- 
cules, the density of which is determined from the measured ambient air 
pressure. At an intermediate height in the atmosphere (8-14 km) where the water 
vapour mixing ratio is beginning to fall toward the low values observed in the 
stratosphere, the fluorescence detector output is compared to the absorption 
measurement, to provide an in-flight calibration for the detector sensitivity 
and geometry, for the fluorescence measurement. The fluorescence technique is 
used to provide the measurement of water vapour in the stratosphere since the 
remaining absorption in the light path of the instrument at those altitudes is 
not large enough to detect. 

A microprocessor in the instrument provides control and data collection 
functions. The system can make a useful measurement in 1 to 5 seconds, which 
provides good height resolution for both ascents, at approximately 5 m/s, and 
descents, at speeds up to 50 m/s. 

A number of possible sources of instrumental error have been addressed by 

the authors. Because the instrument is mounted with a minimum of obstructions 

between the outside air and the fluorescence area, there is the possibility of 
starlight entering the detector. In addition there will be a small Rayleigh 
scattering signal from the light source. These two effects are measured in 
flight and corrected. By observing the fluorescence region periodically with 
the light source turned off, the stray light from outside may be measured by the 
detector. For another interval a quartz shutter is moved in front of the light 
source to ensure that light from longer wavelengths is not making a significant 
contribution to the detected light level by scattering from the air in the 
fluorescence region. The quartz is transparent in the near ultra violet region 

but does not transmit the dissociating radiation which causes the hydroxyl 

fluorescence signal. 
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Observations have been made on both the ascent and descent phases of 
balloon flights. Analysis of the results have lead the authors to an under- 
standing of both the sample contamination problan, and the problem of outgassing 
from the instrument itself. Suitable strategies have been devised to overcome 
these problems, and the later results provide useful measurenents of the strato- 
spheric water vapour content. The flight the authors refer to as MW-5 (KLEY, 
1979) made from Laramie, Wyoming, shows a water vapour maximum of 2.4 ppmm at 
18 km, which is approximately 2 km above the tropopause, rising to 4.8 ppmm at 
30 km. 


In summarizing the features of thi6 instrument, it may be stated that it 
has a major advantage compared to those instruments which have nonlinear 
properties at low vapour pressures of water. The absorption cell mode of 
operation will clearly operate well for large water vapour mixing ratios, and 
the transition region to low amounts of water provides a calibration for the 
fluorescence cell mode. This in-flight calibration removes the uncertainties 
that can develop in the use of instruments which must be pre-flight calibrated 
in the laboratory. 

While the researchers have made a thorough attempt to reduce or eliminate 
the effects of contamination, the possibility still ranains that the measure- 
ments at higher levels in the stratosphere may still contain the effects of 
contamination, as indicated in their paper. 

The good vertical profile resolution of which the instrument is capable, is 
one of its most useful characteristics. Obviously the fundamental accuracy of 
the results is tied to the data used to calculate the absorption by water vapour 
and ^olecular oxygen, and the radiative lifetime and quenching rate for the 
OH (A I ) produced in the photodissociation. The authors estimate the overall 
accuracy in their final results to be about 20%, and the precision to be 10%. 

LONG PATH SOLAR INFRARED ABSORPTION MEASUREMENTS BY SPECTROMETER 

Many investigators have chosen to determine the humidity of the strato- 
sphere through the application of the solar occultation technique in the 
infrared region of the spectrum. For example, MURCRAY (1966, 1969a) has 
measured water vapour in solar absorption in the 2.6, 6.3 pm region of the 
spectrum. ZANDER (1973) also measured in the 2.6 pm region, though at higher 
resolution. GIRARD (1963) and LOUISNARD (1980, 1983) have made measurements of 
water vapour using their Grille spectrometer. All three instruments measure the 
presence of water vapour through its absorption of solar radiation. Profile 
information is devised by consideration of the solar occultation geometry, and 
the growth of the observed absorption with solar zenith angle. The column 
amount of water vapour above the balloon may also be determined in regions of 
the spectrum where the water vapour features are not confused with solar lines, 
and the exoatmospheric spectrum may be inferred from the observations. 

While MURCRAY (1962) has reported results of solar observations in the 
1-1^ pm region from balloon platforms using a relatively low resolution (5-10 
cm ) prism spectrometer, the instrument described here is a later development 
of the Denver group (1963-1964) as described by MURCRAY (1967) in Applied 
Optics. 

_Jhe design target of the spectrometer was a resolution on the order of 
1 cm in the 2.7 pm region of the spectrum. The instrument is a double pass 
0.5 m Czemy-Turner grating spectrometer with a prism predisperser to provide 
order selection. The dispersion element is a 10.2 cm square Bausch and Lomb 
replica grating with 75 lines/nm, blazed at 12 pm. The instrument has short, 
straight slits. A signal chopper is mounted to chop the light path between the 
two passes through the instrument, as shown in the diagran in Figure 3. The 
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Figure 3. Optical diagram of the Czerny- Turner grating 
spectrometer and prism predisperser. The mirrors (M) 
are numbered sequentially. The incident radiation 
enters the spectrometer at slit , is predispersed 
with a Littrow prism (P) with an aluminized rear sur- 
face and passes into the grating spectrometer at S. . 
The grating is double passed using a corner reflector 
(My, Mg) at the focal point of the first pass. 

The radiation is mechanically chopped with a tuning 
fork (C) between the two mirrors of the corner re- 
flector. The radiation is then focused onto the 
detector (D) behind the exit slit (S-). From 
MUR CRAY (1967). 


signal resulting from the detection of radiation at the output of the mono- 
chromator is phase detected to separate the thermal background from incident 
radiation. 

With a choice of predisperser prisms and orders of the grating in the main 
spectrometer, a variety of spectral regions may be studied. The instrument's 
optical properties were explored in the laboratory, and then the spectrometer 
was equipped with a sun tracking system and flown on a balloon (MURCRAY, 1966, 
1969a, b). The instrument has flown successfully many times and has been used to 
determine water vapour profiles using lines in the 2.7 pm region and the 6.3 pm 
region of the spectrum. 

The analysis procedures used to convert the spectra gathered, to a vertical 
profile of water vapour, are outlined by MDRCRAY (1973, 1974). The method used 
involves fitting a model atmosphere to the observed data using a line by line 
calculation, using published line positions, line strengths and line width data. 
The profile is adjusted to make the model spectra fit the observed spectra, and 
then the model profile is assumed to be representative of the atmospheric 
profile at flight time. 
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Since the solar occultation process weights very heavily to the altitude 
where the bean is instantaneously tangent during sunset, a considerable amount 
of information may be obtained about the vertical distribution of the absorber 
in the atmosphere. SNIDER (1975a, b) reports the atmospheric model used for 
these calculations. The path lengths are calculated for the U.S. standard 
atmosphere in the program, and the mixing ratio of water vapour found from the 
fitting process, line parameters are taken from the AFCRL (1973) summary as 
reported by MDRCRAY (1974). 

ZANDER (1973) reported measurements of water vapour obtained by long path 
length absorption, using a large Ebert-Fastie spectrometer illustrated in Figure 
4. The instrument is a 2.5 m focal length monochromator operated in the double 
pass mode, with a narro_%^ intermediate slit. The resolution of the spectrometer 
is reported at 0.04 cm 1 , and the data gathered for the above paper were taken 
in the 1.81 to 1.89 pm and 2.46 to 2.83 pm regions of the spectrum. A lead 
sulphide detector was employed for the measurenent set reported in that paper 
(ZANDER, 1973). 

The instrument is pointed at the sun to an accuracy of ±2 min of arc with 
a 40 cm Ritchey-Chretian telescope with a servo controlled mirror. The whole 
payload is roughly oriented in azimuth. Because of the high resolution 
capability of this instrument, and the high pointing accuracy of the sun- 
tracker, the shapes of the lines measured may be determined. This allows some 
attempt at profiling the amount of an absorber above the balloon to be made 
through the difference in collision broadened line widths of the measured lines 
as a function of altitude. Terrestrial lines may be easily separated from solar 
features because of the considerable difference in their widths. 

% 

The results obtained by ZANDER (1973) indicate a mean mixing ratio of 
(3.5 1 1.5) ppmm above the balloon at 25 km. In order to eliminate the effects 
of water vapour inside the spectrometer itself, a tungsten lamp was flown which 



Figure 4. Optical layout of the Ebert spectrometer. 
The mirrors M. through Mt form a servo-controlled 
telescope used to image the sun on the spectrometer 
entrance slit at S. Mirrors to and the 
grating G form an Ebert design spectormeter . The 
incident radiation is passed through a chopper (CB) 
and is detected by the detector (D). After ZANDER 
(1966). 


70° ELEVATION (me.) 


0° ELEVATION (win) 
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could be turned on in flight and scanned to allow the amount of water vapour in 
the instrument to be measured and subtracted from the column amount calculated 
from the solar observations. 

The water vapour measurements were based on the intensities of 80 isolated 
water vapour lines. The author makes the point that the results are totally 
dependent on a good knowledge of spectral data for water vapour. Since the 
column amount determined here is for the path from the instrument upward to the 
sun, there will be a path contribution due to water vapour emitted by the 
balloon and payload. The effects of this contamination problem are discussed by 
ZANDER (1966) in an earlier paper. 

INFRARED EMISSION RADIOMETRY 


The downward flux of radiation in the stratosphere due to the emission of 
radiation from vibrationally and rotationally excited water vapour has been used 
by a number of investigators to measure the stratospheric water vapour mixing 
ratio. The instruments used for this measurement technique range from cooled 
spectro-radiometers and filter radiometers, to submillimeter interferometers. 

BROOKS (1972) describes a small grating spectrometer designed to be cooled 
to liquid nitrogen temperatures and to observe thj spectral region 16-30 pm. 

The instrument has a spectral resolution of 2 cm near 25 ym and utilizes 
a liquid helium-cooled copper-germanium detector. The observed region of the 
spectrum exhibits features due to CO., HNO^ and water vapour. In the range 
24-30 van only water has a significant emission. 


Observations are made during a balloon descent to approximately 30 km. The 
resulting radiance versus height data are fitted in a model atmosphere using 
data from the AFCRL atmospheric transmittance program, in a line-by-line cal- 
culation. The amount of water vapour observed during this particular flight was 
about 0.6 ppmm around 17 km, and 4 ppmm at 25 km. 


WILLIAMSON (1965) described the use of a liquid nitrogen cooled wideband 
radiometer which sensed thermal emission in the spectral range 5. 5-8.0 pm 
(1820-1250 cm A ). A band calculation following the method of GOODY (1964), 
was used to estimate the amount of water vapour in the lower stratosphere at 
3 ppmm. A large residual signal at the highest flight altitude of 25 km was 
inconsistent (by 3 orders of magnitude) with estimates of the overburden of 
water vapour at that height. 


Because of the difficulty in understanding the WILLIAMSON (1965) data, PICK 
(1969) developed a new radiometer. Figure 5 illustrates the principhl elements. 
There are several notable features in the radiometer design. All the optical 
elements in the field of view of the detector, and the detector itself, are 
cooled to nearly the tenperature of liquid air by radiation to the walls, which 
form the liquid air reservoir. In addition the gas evolved as the radiometer 
ascends or is warmed by heat entering the reservoir, is exhausted through the 
opening through which radiation is allowed to enter. This eliminates the need 
for a cooled window to prevent the accumulation of ice inside the instrument. 


The detector surface is imaged by reflecting optics to accept radiation 
entering the instrument through a field of view of 15° half angle. A motor 
driven filter wheel with six interference filters, and a brass blank is situated 
in front of the gold doped germanium detector. The radiation impinging on the 
detector is chopped at a point between the detector and the filter wheel. The 
filters have a spectral resolution of approximately 0.3 pm. The instrument is 
flown so the center line of its field of view is mounted at 45° to the zenith. 


Inflight radiometric calibration is accomplished by periodically moving a 
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Figure 5. The radiometer assembly shown is a cross 
section along a diameter of the cylindrical copper 
can which holds the cryogenic coolant. From PICK 
(1969). 


blackbody at ambient temperature (measured) into the instrument's field of view. 
Because the radiation passes through the interference filters at a range of 
angles up to 15°, the wavelength passband is broadened toward longer wavelengths 
than the nominal one for each filter. This necessitated that the instrument be 
calibrated against a Golay cell detector in monochromatic radiation on the 
ground, to determine its true spectral response. The main systematic error is 
reported as being due to the error in the determination of the inflight black- 
body temperature. A 1% uncertainty in blackbody temperature translates to a 4% 
uncertainty in radiance measured. 

The measurements confirmed Williamson's measuranent of anomalous fluxes in 
two spectral intervals measured. The amount of water vapour was calculated by a 
procedure identical to that of the previous paper, and was reported as 3 ppmm, 
constant up to 25 km in the stratosphere. The anomalous radiation has since 
been identified as due to nitric acid. 

LONG PATH SOLAR ABSORPTION MEASUREMENTS WITH INFRARED INTERFEROMETERS 

The measuranent procedures of FARMER (1980) and NIPLE (1980) provide good 
examples of the use of interferometry in the infrared region of the spectrum. 

The instruments are both Michelson Interferometers, but are rather different in 
design philosophy and execution. Both Murcray and Farmer are currently very 
active in the measuranent of stratospheric species and they are therefore 
particularly good examples to cite as users of the interferometric technique. 

The use of a high resolution Michelson at relatively short wavelengths in 
a balloon-based situation is a challenging task. Each group has solved the 
difficulties of alignment, data collection and analysis in its own way. 

The basic interferometer used by FARMER (1980) is described by SCHINDLER 
(1970). The instrument i6 a modified Michelson interferometer in which the 
usual plane mirrors have been replaced by a pair of cat's-eye retroref lector s. 
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The retroref lector s are used in a double pass mode where a beam incident on one 
side of the cat' s-eye exits on the other, reflects from a stationary plane 
mirror and retraces its path back through the cat' 8-eye. The arrangement is 
illustrated in Figure 6 which was taken from Schindler's paper. 

The use of the cat' s-eye reflector in this mode has some advantages. 

Because of the double pass which occurs in the path to the moving cat' s-eye, the 
path difference will be double that of a conventional Michelson arrangement. 
Because the cat' s-eye is a retroref lector , the arrangement is relatively 




Figure 6. The optical layout of the interferometer (top) and 
the foreoptics arrangements. Taken from SCHINDLER (1970). 



19 


tolerant of small angular misalignments. In addition, since the beam from the 
beam splitter and the beam from the stationary plane mirror are parallel, the 
system is also tolerant of lateral misalignments of the cat' s-eye mirrors as 
well. These features tend to make the complete instrument generally resistant 
to small misalignments. 

In order to accurately measure the path difference through which the moving 
mirror has travelled, a monochromatic light source, a helium-neon laser, is 
imaged through the interferometer onto a photodetector. This signal is 
modulated by the motion of the moving mirror. 

In fact the control system is much more complicated than this. A moving 
coil actuator is used to make the large physical movement needed for the motion 
of the moving cat' s-eye. Fine control on the order of a wavelength is provided 
by a ptezo-electric stack which makes small adjustments to the secondary mirror 
of the cat' s-eye. An elaborate digital control system steps the mirror assembly 
from one laser wavelength position to another in either direction. The 
necessary sensing signals for this process are derived from the phase sensitive 
rectification of a small amplitude modulation of the laser interference signal 
by a high frequency oscillation of the secondary mirror of the stationary cat's- 
eye. This process generates a signal which is 180° out of phase with the 
detected total intensity of the laser signal, since the modulation depends on 
the differential of the laser amplitude, which is a sinusoidal function of the 
moving mirror position. The use of the total intensity and demodulated 
differential signal together allows the electronics to determine the position 
and direction of motion of the moving mirror at any time. It is, therefore, 
possible to command the control system to move to a particular fringe position. 

Finally, to allow the point-by-point normalization of the interf erogram of 
an external source to its integrated intensity, the interferometer is rapidly 
(1000 Hz) chopped between the observed source and a local blackbody reference 
source. While this chopping occurs, complementary chopping is done between the 
same two sources for a total power detector included to provide the normaliza- 
tion signal. The chopper is a tuning fork type and because of its small size, 
it is situated at the focus of collimating optics so the beam that needs to be 
interrupted is also very small. 

The quoted performance fjr this instrument (FARMER, 1980J is as follows: 
spectral range, 1800-7900 cm ; unapodized resolution 0.3 cm ; time for 
one scan, 200 seconds; detector, InS^; liquid nitrogen cooled. 

The water vapour profiles are recovered from measurements of lines in the 
2\>2 band at 3152 cm 

The amount of water vapour in the path to the instrument was determined by 
a line-by-line calculation of synthetic spectra, which were matched to the 
observed data. Twelve lines in the water vapour spectrum were used to do the 
analysis. The spectral line parameters were taken from McCLATCHEY et al. 

(Air Force Cambridge Research Laboratory, 1973). 

The results of aircraft measuranent flights were analyzed and published, 
with figures available for the amount of water vapour in both the Northern and 
Southern Henispheres. The Southern Henisphere mixing ratios was determined as 
1.2 ppmm, at 22 km, rising to 3.7 ppmm at 36 km. In the Northern Hemisphere the 
22 km mixing ratio was 3.1 ppmm, rising to 3.7 ppmm at 36 km. 

The Murcray group (NIPLE; 1980) has reported the acquisition of high 
resolution spectral data obtained by the use of a Michelson interferometer as 
well. They have begun the systematic application of a novel design interfero- 
meter to the measuranent of atmospheric constituents m the solar occultation 
mode. 
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The instrument was designed and constructed by Bomem Company (BUIJS, 1980) 
and is a conventional plane mirror Michelson interferometer in its optical 
design* In order to construct an instrument with dependable alignment 
characteristics, which is essential for remote operation, the interferometer is 
equipped with a laser servo control system which maintains the alignment of the 
interferometer. The layout of the instrument is shown in Figure 7* A small 
portion of the aperture of the interferometer is used for the detection beam. 
This beam is derived from the helium-neon laser which is fitted with a beam 
expander* A lens is placed in the output beam of the interferometer which 
images straight line fringes localized at the stationary mirror position (BORN, 
1975). These fringes have their maximum density along the line of maximum 
relative tilt between the two Michelson mirrors* Since one mirror is moving, 
the fringes travel across the field in the direction of maximum density* 

Three detectors are placed in the field of view and the output of two are 
phase detected relative to the third. A signal which is proportional to the 
amount and direction of the relative mirror tilt is developed* This signal is 
then used to drive a two axi6 servo system, aligned with the detector axes, to 
maintain the ^stationary** mirror accurately in plane alignment with the image of 
the moving mirror* 


TWO AXIS 

ADJUSTING MOTORS 



COLLIMATEO 
He -Ne 

LASER SOURCE 



Figure 7. The Bomem interferometer* The tilt of the "stationary** 
mirror is adjusted using information from the laser fringe 
detection system so that the image of its surface in the beam 
splitter r mains dynamically parallel to the surface of the 
moving mirror. 
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Clearly, if the phase difference (tilt) between the two detectors i6 a 
multiple of one full wavelength of the laser light, the servo-system will be 
nulled. In order to ensure that this occurs for the plane parallel condition, 
and not some multiple of the one wavelength tilt, the amplitude of the laser 
signal is sensed. The system searches in a preprogrammed manner, to make the 
amplitude of the modulation from the interference large enough to satisfy a 
preset criterion for correct alignment. 

Since the wavelength of the laser light is much shorter (0.6 pm compared 
to 6 pm) than the wavelengths measured, this alignment procedure will ensure 
correct alignment of the mirrors over their full field for the infrared, even 
though the laser detectors take up only a fraction of the full instrumental beam 
width. In aligning the instrument, procedures are used to align the infrared 
detection optical axis with that of the laser alignment subsystem. This ensures 
that when fringe detection is used to "measure" the position of the moving 
mirror relative to zero path difference, the wavelength scale of the resultant 
interf erogram transform will be directly in accurate wavelength units. These 
units will be correct, independent of the refractive index of the environment of 
the instrument. 

The use of an automated alignment system results in an instrument of 
unusually small size and light weight (about 50 kg). It is relatively immune to 
isolated shocks since it realigns itself after they have occurred. It is not 
particularly insensitive to vibration while operating, but since balloon pay- 
loads are relatively free of vibration, the instrument provides good performance 
in the balloon environment. The automated alignment system greatly reduces the 
thermal stabilization time, and also provides the important information as to 
whether the instrument is properly aligned at any given time, so that data 
analysis can be properly done. 

The instrument described by NIPLE (1980) is capable of simultaneously 
measuring njtrogen dioxide and water vapour. It has an unapodized resolution 
of 0.02 cm , and can produce a complete interf erogram in approximately 50 
seconds. 

The 1603.0-1607.4 cm * region of the spectrum was observed for the 
detection of both water vapour and NOj. This region includes 232 N0 2 lines 
and 19 HjO lines. 

Because of the optical characteristics of the instrument (channel spectra 
and instrument apodization) it was decided to improve the analysis procedure by 
including the known instrumental properties explicitly in the analysis 
procedure, rather than attempting to remove or explain them in their resultant 
effects on the transformed spectra. 

This is accomplished by determining the infinite resolution spectrum 
appropriate for some first guess model atmosphere, using published spectral 
data for the constituents considered. The high resolution spectrum is Fourier 
transformed, and the transform is modified according to the known properties of 
the observing instrument. Both the observed interf erogram and the simulated one 
are then transformed to spectral space and compared. An interaction scheme is 
operated to correct for differences between the two, using a least squares 
metric to quantify the fit. The model atmosphere producing the spectrum with 
minimal error is taken to describe the real atmosphere at the time of the 
observations. 

The flight described by NIPLE (1980) determined the stratospheric water 
vapour mixing ratio to be 2.4 ppmm at 17 km, rising to 3.9 ppmm at about 30 km. 
This paper references various others which describe in detail the algorithms 
used for the analysis of the interf erogram s obtained. 
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FILTER RADIOMETERS FOR ROTATION BAND MEASUREMENTS 

BREWER (1972) reported the development of a long wavelength filter radio- 
meter designed to observe the emission from water vapour in the 20-60 pm region 
of the spectrum. The instrument design is illustrated in Figure 8. It consists 
of a Barnes thermistor bolometer on which radiation is imaged after passing 
through a waveband selecting filter. The filter was compounded from beryllium 
oxide and carbon black, which were embedded in polyethylene. The filter is 
transparent to wavelengths longer than 20 pm. 


MOTOR 



Figure 8. The Brewer-Thomson radiometer. The detector 
is successively exposed to radiation from the sky, an 
ambient temperature black body, and the liquid nitrogen 
reservoir. The detector is chopped by a salt crystal 
chopper driven by a stepping motor. From BREWER (1972). 

The signal detected by the bolometer is chopped by a 32 Hz stepping motor 
chopper. The chopper blade is made of sodium chloride, since the polyethylene 
filter has a small transmission around 5 pm. Sodium chloride transmits in that 
region of the spectrum, but cuts off rapidly to wavelengths longer than 20 pm. 
Therefore the resultant output of the phase detector is due to radiation of 
wavelengths longer than 20 pm. The signal is then phase detected and tele- 
metered to the ground. The chopper temperature in flight is not known, so the 

collimating mirror which images the sky on the detector is turned periodically 
to sample either a black body at ambient temperature, or a flask of liquid 
nitrogen which is carried along with the radiometer. 

A standard radiosonde is flown along with the radiometer to measure 

temperature and pressure. The instrument is light weight and can be flown on a 

2.5 kg balloon or two radiosonde balloons. 

Three ascents were made with the radiometer and yielded a lower strato- 
sphere water vapour mixing ratio of 2 ppmm. 

HYSON (1974a, b; 1978) describes a small radiometer for the detection of 
water vapour in the stratosphere. The instrument is similar in principle to 
that of the Brewer, but has a few interesting differences. The instrument is 
illustrated in Figure 9 taken from HYSON (1976). 
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ZENITH 



Figure 9. The Golay cell receives radiation 
from one of two sources; the intensity of 
radiation from the sky above the radiometer 
is compared to the (assumed) saturated 
signal from emission by water vapour at 
ambient temperature below the horizontal 
direction. The sources are compared through 
the use of a reflective choppier. After 
HYSON (1974b). 


The instrument uses a Golay cell detector, which is alternately illuminated 
by radiation from the sky above the radiometer, or the atmosphere below at a 
depression angle of 16° relative to the local horizontal. The radiation from 
the sky may be examined at zenith angle of 45° or 70°. 

The instrument uses a polyethylene filter to pass radiation in the 40-200 
ym region of the spiectrum to the detector. The detector beam is choppied with a 
reflecting chopper between the upward looking and downward looking mirrors. The 
chopping rate i6 14 Hz and the Golay cell signal is phase detected. 

Since the amount of water vapour below the instrument is very large, the 
signal arriving from that direction is taken to be representative of a blackbody 
at ambient temperature. 

Water vapour amounts in the Australian stratosphere have been reported by 
Hyson. It was found that the water vapour mixing ratio is between 2.3 and 3.1 
ppimm in the altitude range 15 to 21 km and approximately 5 ppmm at 25 km. 

DIFFERENTIAL ABSORPTION RADIOMETER 

B0RKERT (1974) describes a filter radiometer for the measurement of water 
vapiour in long path solar absorption. The instrument is shown in Figure 10. 

When a gas is measured in absorption, the measuranent depends on the detection 
of the amount of radiation removed from the beam on passing through the 
absorbing gas. Since gases, in general, do their absorbing in very narrow 
lines, the detection of the net depletion of a bean by a broadband sensor will 
show only small changes for substantial changes in the amount of absorber in the 
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Figure 10. Signals are alternately detected by the detector D, 
which originate from two telescopes formed by mirrors 
and M_. The telescopes are fitted with two narrow-band 
interference filters; one chosen to lie on an appreciable 
water vapour absorption feature, and one nearby where the 
absorption from water vapour is negligible. After BURKERT 
(1974). 


beam. This signal may be lost in noise sources such as fluctuations in the 
source intensity, or changes in the detector sensitivity. The signal to noise 
ratio may be improved by using spectrometers or interferometers of higher 
resolution, but with the attendant costs of increased size, complexity and the 
loss of speed. 

A filter radiometer is a very small and simple instrument, but suffers from 
the signal to noise problem discussed above. It is possible to gain some of the 
signal to noise advantage of a high resolution device, while maintaining the 
simplicity of the filter radiometer. With reference to Figure 10 taken from 
BURKERT (1974), it is possible to select narrow band interference filters (50 
cni) chosen to lie on, and just beside a water vapour absorption feature. If 
incident radiation is allowed to reach a detector after being absorbed by these 
two different filters, the difference in the amount of radiation taking these 
two paths will be more sensitive to the absorber than either of the two beams. 

This technique will automatically reduce the dependence of the final 
instrument signal on the source intensity or the detector sensitivity (BURKERT, 
1974). In order to calculate the amount of absorber on the solar bean arriving 
at the instrument, it is necessary to know what the instrument output would be 
if there were no absorbing gas in the atmosphere. By allowing the instrument to 
observe the sun at high sun angles, at float altitude for some time, a Langley 
plot may be prepared (a plot of instrument output versus atmospheric path 
length) and the intercept of this plot will be the value the instrument would 
read with zero absorber amount in the atmosphere. 

Thi6 extraterrestrial value, together with the known instrumental response 
to changes in absorber amount, will allow the determination of the total amount 
of absorber in the path to the sun. 

On the balloon flight reported in BURKERT (1974) two filters wereysed to 
select narrow bands at 2.606 pm and 2.280 pm with a bandwidth of 50 cm • 

The extraterrestrial constant for the instrument was determined during a float 
interval at 29.5 km. This, together with data taken during balloon descent, 
allowed the calculation of the total amount of water vapour in the incident path 
as a function of altitude. The conversion from instrument reading to water 
vapour amount was done by using the laboratory measured spectral response of the 
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instrument in a line by line calculation, assuming a Voigt line profile for the 
absorption by water vapour. By calculating from float height downward a 
vertical profile can be determined. 

The published profile exhibits a minimum mixing ratio of (0.8 ± 0.2) ppmm 
at 20 km, which is about 2 km above the tropopause for that flight. The mixing 
ratio rises to about 8 ppmm by 30 km but the uncertainty due to possible 
instrument out gas sing is approximately ±3 ppmm. 

RESULTS OF WATER VAPOUR INVESTIGATIONS 

It is appropriate to introduce a discussion of the results of water vapour 
investigations by following ELLSAESSER (1983) in quoting GUTNICK (1961): "Lack 

of a tenable explanation of the various reports on stratospheric humidity casts 
doubt on the accuracy of part or all of the data". Figure 11 is replotted from 
ELLSAESSER (1983) and is representative of recent data. The large variability 
of the measurements is difficult to reconcile with the Brewer hypothesis for the 
circulation of air through the stratosphere. Not the least of the difficulties 
is the variation in the vertical gradient of water vapour amount, which ranges 
from near zero to a change of 5 ppmm between the tropopause level and 30 to 40 
km altitude. 

In 1961 Gutnick called upon the scientific community to carry out a co- 
operative research program to assess the size and sources of the differences 
between measurements made using different techniques and different instruments. 
Just twenty years after this call, the Federal Aviation Administration (FAA) of 
the United States and the National Aeronautics and Space Administration (NASA) 



Figure 11. A comparison of some water vapour measurements 
made from stratospheric research balloons. After 
EL1SAESSER (1983). 




26 


agreed to jointly administer a water vapour inter comparison program. The goal 
and details of the program were reported in the UPPER ATMOSPHERIC PROGRAMS 
BULLETIN (1981). Table 1 lists the cooperating investigators. The set of 
instruments chosen to participate in the program provides a cross section of the 
techniques used in the past for water vapour measurement. 


When the preliminary results are examined (UPPER ATMOSPHERIC PROGRAMS 
BULLETIN, 1981, 1982) the agreement between the different techniques is seen to 
be reasonably good (Figure 12). This result helps to clarify some of the 
uncertainties raised by ELLSAESSER (1983) in his review. 



H^O mixing ratio (ppmv) 

Figure 12. Preliminary results from the 1981 
comparison of stratospheric water vapour 
instrumentation. Published in the UPPER 
ATMOSPHERIC PROGRAMS BULLETIN (1981, 1982). 

It will be noted that in situ and remote sensing instruments agree on a 
difference in the vertical of 1.2 to 1.5 ppmm between the minimum mixing ratio 
level near 17 km and an altitude of 30 km. This is in contrast to the long- 
standing difference between past MASTENBR00K (1974b) results taken with the 
frost point hygrometer and those observations made by remote sensing measure- 
ments made by the solar occultation technique (MURCRAY, 1966; FARMER, 1980). 

This would suggest that there is a much larger variability in the water vapour 
mixing ratio profiles with respect to time and place than had previously been 
expected. 

The preliminary results (1.7 + 0.4 ppmm) are compatable with the previously 
accepted levels of water vapour measured by Mastenbrook, which were in the 
region of 2.0 ppmm in the lower stratosphere. This level of humidity is 
comparable to that which NEWELL's (1981) stratospheric fountain would introduce 
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Table 1. A partial list of investigators of the water vapour mixing 
ratio in the stratosphere. 



MEASUREMENT 


ALTITUDE 

MIXING RATIO 

INVESTIGATOR 

TECHNIQUE 

LOCATION 

(km) 

(ppmm) 

BARRETT 

(1972) 

Emission 18-23 pm 
Downward Emission 


low strat. 

2. 4-2 .6 

BREWER 

(1972) 

Emission 20-60 pm 
Downward Emission 

Ontario 

Canada 

(1969-70) 

low strat. 

1. 4-3.0 

Bonne tti 

Emission 


18 

2.5 

(BANGHAM; 1980) 

5-80 cm-1 
Long Path 


36 

5 

GOLDMAN 

Emission 25 ym 

New Mexico 

16 

<1.0 

(1972) 

Ascent 

USA (1971) 

25 

3-4 

HARRIES 

(1973) 

Submillimeter 
Emission 
Long Path 

England 

(1971) 

mid strat 

1-3 

HYSON 

(1974) 

Emission 30 ym 


15-18 

4 

KLEY 

Lyman Alpha 


17 

2.4 

(1979) 

Photodissociation 


30 

4.8 

LOUISNARD 

Absorption 5-6 pm 

CNES 

20 

1.8 

(1983) 

Solar Occultation 

France 

40 

5.5 

MAS TEN BROOK 

Frost Point 

Washington 

18 

2 

(1974) 

Hygrometer 

D.C., USA 
(1960 - ) 

22 

3 

MUR CRAY 

Absorption 2.6 ym 

New Mexico 

16 

<2.0 

(1966) 

Ascent 

USA (1964) 

30 

>10.0 

MURCRAY 

Absorption 2.6 ym 

New Mexico 

25 

2.5 

(1969) 

Solar Occultation 

USA (1964) 

30 

3.0 

MURCRAY 

Emission 2.5 ym 

Hyderabad 

16 

4-5 

(1973) 

Ascent 

India (1972) 

25 

<1.0 

PICK 

Emission 

England 

low strat 

3.0 

(1969) 

5.8-6 .8 ym 
Ascent 

(1967) 

25 

3.0 

WILLIAMSON 

Emission 

England 

low strat 

3.0 

(1965) 

5.5-8 .0 um 
Ascent 

(1963) 

25 

3.0 

ZANDER 

(1973) 

Absorption 
1.8, 2.5 ym 
Solar Occultation 


15-25 

3.5 
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into the stratosphere. The observed amount of water vapour leaves little margin 
for any other sources of water, such as tropical hot tower convection 
(ELLSAESSER, 1983) which might inject additional water into the stratosphere. 

The results of the inter comparison also seem to be consistent with the observa- 
tion by KLEY (1979) of a hygropause, which is a mixing ratio minimum in the 
water vapour profile located a couple of kilometers above the tropopause height. 

There remains the possibility that some other mechanism for drying air, 
both in the troposphere and in the stratosphere, remains to be elucidated 
(ELLSAESSER, 1983). 

OZONE 

The limited extent to which vertical profile information may be obtained 
from ground-based measuranents of ozone (GOTZ, 1934) lead investigators to 
devise methods for the direct sampling of the ozone distribution. REGMER 
(1938) flew an ultra violet spectrograph on board a balloon payload to make a 
determination of the ozone vertical profile. Similar measuranents were reported 
by PAETZOLD (1955). .Since these techniques produce only moderate vertical 
resolution of the ozone profile and are little used, they will not be discussed 
here. 

REGENER OZONE SONDE 

REGENER (1960, 1964) developed a chemiluminescent sonde to make a direct 
determination of ozone concentration through a chemical detection process. This 
sonde depends on the reaction of ozone in the atmospheric sample with 
Rhodamine-B. This substance reacts with ozone with a characteristic 
fluoresence, the intensity of which is related to the ozone concentration. 

The fluorescence signal is detected with a photomultiplier to produce an output 
which is used as an indicator of the ozone concentration in the measurenent 
chamber. 

In use as a stratospheric sounder the measurement device is aspirated with 
a diaphragm pump which draws outside air into the reaction vessel. The 
resulting signal is telemetered to the ground via a standard radiosonde, which 
is flown along with the instrument to measure the background atmosphere. While 
the instrument wa6 used for a number of years by the United States Weather 
Service, it was eventually phased out. While the dry chemistry involved in the 
sonde was an advantageous simplification, the lack of an absolute calibration, 
difficult preparation and poor instrument stability were among its problems. 

The instrument was also difficult to use operationally since a preflight 
calibration with an accurately determined ozone source was required. 

BREWER OZONE SONDE 

BREWER (1960) reported the development and use of a wet chemical method for 
the detection of ozone. The method depends on the ability of ozone to release 
iodine in a solution of potassium iodide according to the reaction: 

2KI + 0 3 + H 2 0 -*• 2K0H + I 2 + °2 < 10 ) 

Two electrodes are placed in the electrolyte solution. If a current is 
allowed to flow, the following reactions will take place at the anode and 
cathode respectively: 

I 2 + 2e -»■ 21" (11) 


2I - -*- I 2 + 2e 


( 12 ) 
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Actually the I exists in the form 1^ , but the net result of the 
chemistry is as indicated here. If the I. produced is free to recirculate, a 
current may be maintained vithout the continual addition of iodine. In the 
Brewer system this re-circulation is prevented by using one electrode of 
platinum and one of silver. The silver electrode is the anode, and the iodine 
formed at the anode by reaction R3 is immediately converted to Agl which is 
insoluble in the electrolyte solution. The net result is the formation of an 
electrochemical systaa which allows the circulation of a current of two 
electrons in the external circuit, for each molecule of ozone dissolved and 
reacted in the electrolyte according to reaction Rl. 

In the light of recent knowledge of the behaviour of ozonesondes in the 
presence of sulphur dioxide and other atmospheric pollutants, it is interesting 
to note that Brewer observed the presence of an absorber which caused an effect 
in the sense of "negative" ozone anounts in some of his early work, and he 
concluded that the observed effect was "clearly due to negative shifts of the 
telemetry system". In addition he noted that the response time of the sonde 
seemed to be a function of altitude, an effect which was studied in more detail 
by DE MUER (1980). This determination of the sonde response time is not 
compatible with the response times determined in the laboratory. 

The advantage of the Brewer sonde is that it needs no absolute calibration. 
Operationally, however, the flights are usually normalized to a nearby total 
ozone reading made by a Dobson or Brewer spectrophotometer. The magnitude of 
the correction coefficient is taken as indicative of the absolute accuracy of 
the ozonesonde in performing a measurement of ozone. 

THE ELECTROCHEMICAL CONCENTRATION CELL 

KOMHYR (1967, 1969) reported the development of an electrochemical concen- 
tration cell for the measurement of ozone. The cell system is made up of 
several components and is depicted in the sketch in Figure 13. The cell 
consists of anode and cathode containers and a larger reservoir for the cathode 
electro lyte. 


DETAIL OF CONSTRUCTION OF 
ELECTROCHEMICAL CONCENTRATION CELL 
FOR GAS ANALYSIS 



Figure 13. Details of the construction of the electro- 
chemical concentration cell for gas analysis. From 
KOMHYR (1969). 
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The cathode cell consists of a platinum electrode immersed in a dilute, 
buffered solution of potassium iodide. The anode container is connected to the 
cathode by an electrolyte bridge composed of tightly packed asbestos fibers. 

The bridge prevents bulk flow, but allows the exchange of ions between the 
cathode and anode regions. The anode electrolyte is a concentrated potassium 
iodide solution, and the anode is platinum. 

In the cathode region iodine is reduced according to the following 
reaction: 


I 2 + 2e + 21" (13) 

In the anode region oxidation takes place according to the reaction: 

3 1” -*■ I" + 2e (14) 

The electrons produced or consumed are exchanged through an external 
circuit and constitute the measurement of the iodine consumed by the reduction 
reaction in the cathode area. The measurement of ozone is accomplished through 
the production of iodine in the cathode compartment by the reaction of ozone 
with potassium iodide. This reaction proceeds according to the following 
reaction: 


2KI + 0 3 + H 2 0 2K0H + I 2 + 0 2 (15) 

The ozone i6 introduced into the cathode solution by a small teflon pump 
which drives the air sample into the cathode region. The gas is allowed to 
bubble up through the electrolyte, where the ozone goes into solution and reacts 
with the potassium iodide. 

The chemistry of the electrochemical cell is quite complex (see THORNTON, 
1982; TORRES, 1978) and the cell concentrations were adjusted empirically to 
achieve a calibrated response to ozone. In use it is found that the cell has a 
"zero-signal" with air bubbled through it in the absence of ozone. Thornton has 
identified this as due to the renoval of 1^ which is normally present in 
the cathode chamber immediately after the cell is filled. This constant back- 
ground signal is subtracted from the ozone induced response of the cell to yield 
the measured signal. 

In operational use the ECC sonde i6 connected to a standard radiosonde and 
uses the radiosonde's transmitter to telemeter the ozone profile information to 
the receiving station. In practice it is found that it is useful to normalize 
the ozone profile obtained to a nearby Dobson or Brewer spectrophotometer total 
ozone value. 

In the context of research ballooning the measursnent of the ambient ozone 
concentration is usually done in conjunction with a measursnent flight intended 
to measure other stratospheric constituents. The usual end goal of a measure- 
ment flight is that of providing suitable profile information for the develop- 
ment and verification of photochemical models designed to explain the observed 
properties of the stratosphere. In order to accomplish this goal, it is 
essential to measure the state of the stratosphere at the time that any 
particular trace constituent measurement is made. One reason for this is the 
well known large tanporal and spatial variability of the ozone profile. In 
addition some measurements such as long path infrared absorption measurements 
require a detailed knowledge of the concentration profile of background species 
to allow the calculation of their effects on the observed spectra. 

Generally, the performance characteristics of the Brewer and Komhyr ozone- 
sonde systems are quite similar. It is found operationally that the correction 
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factors needed to bring ground-based total ozone measurements and ozonesonde 
measurements into agreement are slightly different for the two instruments. 

Table 2, taken from MATEER (1977), illustrates the typical mean response of 
Brewer-Mast sondes in relation to Dobson total ozone amount. It is found that 
the ozonesonde amounts must be scaled by approximately 1.2 to 1.3 to yield a 
value for the total ozone which is comparable to that measured by a nearby 
spectrophotometer. A similar table for the ECC ozonesonde is presented as Table 
3 which is taken from KIRCHOFF (1983). It can be seen that the correction 
factor for the ECC sonde is closer to unity for the largest set of observations 
reported than that of the Mast sonde, and the range of the correction factors 
observed is smaller. The ECC sonde is a somewhat better choice for the 
investigator who is supporting a balloon flight without the support of a ground- 
based measurement of ozone. 


Table 2. Typical normalization factors for Brewer-Mast ozonesondes 

compared to Dobson total ozone reading. After MATEER (1977). 

TOTAL OZONE 


LOCATION 

SAMPLE SIZE 

CORRECTION FACTOR 

(DOBSON UNITS) 

Churchil 1 

92 

1.309 ± 0.187 

377 ± 58 

Edmonton 

161 

1.236 ± 0.165 

351 ± 50 

Goose 

244 

1.301 i 0.199 

377 ± 54 

Resolute 

159 

1.186 ± 0.113 

412 ± 86 


Table 3. Statistical parameters of ozone total column amount measured 
by ECC and Dobson instruments. From KIRCHHOFF (1983). 




ECC 

DOBSON 



DIFFER® CE 

PERCENT 
DIFFER® CE 



Geraci and Luers 

(1970-1975 

data, 

95 days) 


Mean 

Standard 

Deviation 

0.319 

0.050 

0.320 

0.034 



-0.001 

0.038 

-0.30 

12.19 



Geraci and Luers (1973 data, 

16 

days) 


Mean 

Standard 

Deviation 

0.338 

0.043 

0.320 

0.024 



0.017 

0.026 

5.25 

7.74 



Natal Data 

(1978-1981, 

43 

day 8) 


Mean 

Standard 

Deviation 

0.286 

0.025 

0.276 

0.013 



0.010 

0.027 

3.84 

9.65 
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DASIBI 

The Dasibi ozone measuring instrument is a widely used device for the 
ground-based monitoring of atmospheric ozone (Dasibi operating and instruction 
manual, Model 1003 AH). The Dasibi instrument measures ozone by measuring the 
extinction of 253.7 nm radiation from an internal mercury lamp. The light is 
made to transverse a known length of absorption tube and the measured light 
intensity is indirectly compared to the intensity the same lamp-detector 
contribution would exhibit if there were no ozone in the absorption tube. This 
normalized measurement is accomplished by comparing the reading the detector 
system recorded when the air sample was passed through an ozone destroying 
filter of manganese dioxide before being admitted to the absorption cell, with 
that obtained when the air sample containing ozone is measured. 

In addition to this comparison procedure which eliminates the effect of 
absorption cell changes, the instantaneous fluctuations of the light source are 
also eliminated by the use of a second detector which views the light source 
without an intervening absorption cell. This detector controls the sensing 
electronics in such a way as to remove the dependence of the final measurement 
on the absolute intensity of the mercury line source. 

The principle of the Dasibi laboratory instrument has been adapted to a 
successful instrument for the measurement of stratospheric ozone. Robbins 
(Johnson Space Centre; private communication) reports the details of the balloon 
version of the instrument. 
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ABSTRACT 

Remote sensing from balloon platforms by visible and infrared spectrometry 
has been extensively used for the past 10 years to measure the concentration 
profiles of stratospheric species. A brief review of the instrumentation and 
method of analysis is made. Recent results concerning HC1, CIO, HF, NOx, 

HNOj, CO, CH^, N 2 O, H 2 O and 0^ are presented. 

INTRODUCTION 

Infrared spectroscopy and radiometry are very powerful techniques for 
remote sensing of trace concentrations of atmospheric gases. The unique spectra 
associated with molecular transitions enables gases to be identified and the 
data obtained from the spectra can be used to determine quantitative informa- 
tion. Two spectroscopic techniques are currently available for remote 
atmospheric monitoring from balloons; these are: 

(1) Absorption spectroscopy which measures the transmission of radiation 

through the atmosphere from a distant source, in most cases the sun; 

(2) Emission spectroscopy which measures the atmospheric thermal emission. 

These two techniques are now well established in the field of atmospheric 
monitoring and it is these which we shall review in this chapter. 

(a) The Mid-Infrared Spectrum of the Atmosphere 

Before reviewing the measurement techniques and instrumentation, let us 
first briefly examine the mid-infrared atmospheric spectrum. Langley in 1884 
was the first to record the transmission of sunlight through the atmosphere in 
the infrared. There have been many atmospheric infrared atlases produced during 
the last century; the reader is referred to the recent atlases produced by the 
Denver group (GOLDMAN et al., 1982) which are compiled from high resolution 
solar absorption spectra measured from the ground and from balloons. 

The mid-infrared region is bounded by intense absorption features due to 
carbon dioxide at 2.7, 4.3 and 15 microns. Water vapour has strong rotation- 
vibration bands at 1-9, 2.7 and 6.3 microns. Other bands due to ozone, methane 
and nitrous oxide occur at 9.6, 7.6 and 7.8 microns, respectively. Between 
these absorption bands, there are regions known as "atmospheric windows" where 
the atmosphere is less opaque. Some absorption in these window regions is still 
present in the lower atmosphere due to the wings of lines in the water vapour 
bands and, at sea level, large portions of the infrared spectrum are very 
heavily, if not totally absorbing. However, detailed spectroscopic investiga- 
tions in the window regions from mountain observations, high flying aircraft and 
balloons have led to the identification of a large number of minor stratospheric 
constituents. The position of the major vibration-rotation bands which are 
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known to exist in the mid-infrared are shown in Figure 1 which also shows the 
spectral emittance of a blackbody at typical stratospheric temperatures. 



MAIN ABSORPTION BANDS 
OF ATMOSPHERIC GASES 


Figure 1. Main absorption bands of atmospheric 
gases in the infrared and spectral radiant 
emittance (W/m 1 cm - ' 1 ) for the sun and two 
typical stratospheric temperatures. 


(b) Measuranents from Balloons 

Balloons made from modern plastics are able to carry payloads of several 
tons to altitudes of 35-40 km. From these altitudes balloon-borne spectrometers 
can make observations over a large range of zenith angles from more than 40° to 
96°, limited on one side by the balloon entering the field of view and on the 
other by the earth. The mass path along the line of sight increases with zenith 
angle; for zenith angles greater than 90° the mass path can be orders of 
magnitude greater than that for a vertical column and so it is ideal for the 
detection of molecules with weak emission or absorption features; it may also be 
shown that most of the absorption or anission is located in a rather narrow 
layer, roughly one scale height thick, above the point along the line of sight 
where the latter is tangential to the earth's surface. The vertical concen- 
tration profile of trace species below the spectrometer can be determined from 
spectral measuranents made at zenith angles between 90° and 96°. This 
technique, known as limb scanning, has been discussed in detail by Gille and 
House (GILLE and HOUSE, ’1971) and is the most important technique for the remote 
monitoring of atmospheric species from balloons and satellites. Sane 
information can also be deduced from spectra measured at zenith angles less than 
90°, for example the mean mixing ratio and the total amount of the species above 
the altitude of the balloon. The main disadvantage of the limb-scanning 
technique is that, because the line of sight is near horizontal, a large 
horizontal scale applies to the measurement. 
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THE EQUATION OF RADIATIVE TRANSFER 

In order to interpret atmospheric emission or absorption spectra, it is 
necessary to understand how an electromagnetive wave interacts with the 
atmosphere. The radiative transfer process through the real atmosphere is 
complex; two assumptions are made to simplify the calculation; (1) scattering is 
not appreciable, so that only processes of absorption and emission are 
important; (2) the atmosphere is in local thermodynamic equilibrium. 

For a plane parallel non-scattering atmosphere, the transfer equation is 
(GODY, 1964): 

dl = -k c(I - B ) ds (1) 

V V V V 

where 1 is the intensity of the incident radiation per unit frequency interval 
per uni^ solid angle, B v is the locally defined Planck function, ky is the 
spectral absorption coefficient, c is the density and 8 is the linear dimension 
along the ray path. 

The solution to equation (1) can be written as: 

1^(8) ■=> I v (o) t^(s,o) + f B v (s') t (s,s') k (s') c(s') ds' (2) 

where t (s.s'), the fractional transmittance from s' to s, is given by: 

f S 

t v (s,s') = exp- k (s) c(s) ds (3) 

' s ' V 

The first term on the right hand side of equation (2) represents the 
radiation from an external atmospheric source I(o) modified by its transmission 
through the atmosphere by the fractional transmittance t(s,o). The second term 
represents the emission, dl, from elanents ds 1 : 


dl = B(s') k v (s') c(s') ds' 


but again modified by the fractional transmittance from s' to s. 
Taking the derivative of t along the line of sight, we find: 
dx (s.s 1 ) = -k (s') c(s') t (s, s ' ) ds 

V v V 

hence substituting (5) in equation (2): 

'° dT v (s,s’) 


M s) “ r v (o) T v ( s .°) + [ V 8 ') 

* « 


ds 


ds' 


(4) 


(5) 


( 6 ) 


This is a convenient form of the equation for numerical computation of radiative 
transfer; the radiation flux at s is simply expressed in terms of t, the 
fractional transmittance, and B^ the Planck function. 

The spectral absorption coefficient, k v , is the key parameter for calcula- 
tion of the transmittance function x(s). For a single line the absorption 
coefficient is: 

k = Sf(v) 
v 

where s, the line strength, is equal to the integrated absorption coefficient, 

S ** / k^dv, and f(v) is the line shape function. 
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The pressure broadened line shape, due to Lorentz, is given by: 

f(v) = — [(v - v ) 2 + y. 2 ] 1 
n o 'ti 

where the halfwidth Y^ is related to the value Y Q at a standard pressure 
I* 0 and temperature T 0 by: 

Y» / Y - (P/Pj * (T/T ) n 

la O O O 

where n, usually assumed to be -0.5, is not well known; depending on the mole- 
cule it could vary between -0.5 and -1.0. 

At low pressures (high altitudes) Doppler broadening becomes dominant and 
the Doppler line shape is given by: 

f(v) “ w7 exp(_(v " V i ) 2 /Y D 2) 
t D 

where y d , the Doppler halfwidth, is: 

Y D = 3.58 x 10 -7 v^T/M) 172 

V. is the central frequency of the line and M is the molecular weight of the 
gas at temperature T. 

In the middle stratosphere the line shape can be represented by a 
convolution of Lorentz and Doppler line shapes known as the Voigt line shape. 

In practice the absorption coefficient at a given frequency is made up from the 
contributions of the individual lines close by. 

We can now go into the details of the two main spectroscopic techniques. 

(a) Absorption Spectroscopy 

A distant source like the sun or a star has an intensity I(o) much higher 
than the emission of the atmosphere. Referring to equation (6), the second term 
on the right hand side can be neglected, hence 

I <8 ) “ I(o) t(s,o) . (7) 

Therefore the signal- to-noise ratio which can be achieved using absorption 
spectroscopy is much greater than that which can be achieved using emission 
spectroscopy. This signal- to-noise advantage can, in classical instruments, be 
traded for increased spectral resolution in order to isolate absorption features 
with little or no contamination from other species. 

Since this technique requires a hot source at the limb the measurements 
are restricted to sunrise or sunset when the sun is used as a source, periods 
which are particularly undesirable for the observation of short-lived photo- 
chemical ly produced species. Measurements using stars as a source are more 
flexible, permitting observations to be made over an extended period of time. 
However the intensity of these sources is considerably weaker than the sun, so 
measurements have been made only in the visible and UV spectral regions where 
their emission is greatest. These measurements are only made at night to avoid 
spurious signals from scattered sunlight. 
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(b) Emission Spectroscopy 

Measurements are made against the cold background of space. The first term 
on the right hand side of equation (6) can now be neglected: 

ts dx 

I (s) = B (s,T) “T~ (8,8* ) d8' (8) 

v I , v ds 

1 s 

The radiative flux, determined by the Planck function, is reduced by dt/ds, 
which is referred to as the weighting function since it "weights" the Planck 
function along the path of the ray and indicates from where the emitted 
radiation originates. Figure 2 indicates the weighting functions versus 
altitude for COj. 



Weighting function (km -1 ) 

Figure 2. Weighting function of limb sounding; 
a single spectral C0_ band has been used 
(585-705 cm _i ') and the different curves 
correspond to different observational zenith 
angles (from G1LLE and HOUSE, 1971). 


This function depends on the concentration and molecular transition 
strengths of the gas. For the more abundant stratospheric gases, such as ozone 
and carbon dioxide, this term can approach unity and the observed thermal 
emission is close to the maximum determined by the Planck function. However, 
for less abundant species this term is often much less than unity. 

Because the radiative flux is determined by the Planck function reduced by 
a weighting function, the following constraints are imposed on the emission 
technique : 

(1) Limitation of the spectral range to the 10-2000 cm ^ where the 
atmospheric thermal emission is greatest; 

(2) The need to know accurately the Planck function and hence atmospheric_j 
temperature. This can be a particularly acute problem beyond 1500 cm 
where the Planck function varies rapidly with temperature; 
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(3) Extremely sensitive instrumentation to detect the low levels of radiation. 
The instruments are often cooled to keep other sources of radiation to a 
minimum. Also elaborate baffle systems are required to prevent radiation 
from the earth entering the instrument when operating in the limb 
scanning mode. 

All these constraints are balanced by the fundamental advantage of this 
technique not requiring a hot source such as the sun. Therefore, unlike 
absorption spectroscopy, measursnents made using the limb-scanning technique are 
not restricted to particular times of the day but can be made at any time of day 
and night and along any compass bearing. 

Hence the emission technique is invaluable for making studies of the 
diurnal variation of stratospheric species concentrations and for studying short 
lived photo-chemical ly produced species. 

INSTRUMENTATION 

(a) General Considerations 

We will characterise the performance of instruments designed to measure 
profiles of minor stratospheric constitutents by two basic criteria which 
broadly determine the accuracy of the measurement: (1) signal-to-noise ratio, 
and (2) spectral resolution. These criteria also determine more specific con- 
straints such as instrument weight and ability to work on stratospheric balloon 
platforms. 

Signal-to-noiae ratio . The incident flux on the detector, S, is always 
a linear function of the following parameters: 

- L, luminance of the source, 

- G, instrumental throughput, equal to the product of the entrance aperture and 
the field of view, 

- T , instrumental transmission, 

- At, integration time for one sample, 

- Aa, spectral interval. 

/ 

S = LGT r At Aa 

The dominant noise can be detector noise (independent of the photon flux) 
or photon noise. 

Until recently the performance of infrared instruments has been limited by 
detector noise while that of UV and visible instruments is limited by photon 
noi se . 

There are two fluxes of radiation incident on the detector: 

- Fg is the background flux, 

- F„ is the source flux. According to the type of instrument as will be 
shown in a later section, the source flux can be composed of two components, 

- Fg u is the signal to be analysed, 

- F_ can be considered as parasitic because it will not be discriminated by 
tnf electronics but contributes to the noise (interferometer, grille spectro- 
meter, ...). 

In the detector noise limited case: 

S/ N - (Fg + F B )/I n 
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while in the photon noise limited case: 

S/N = (F Su + F B )/(F Su + F Sp + F B )1/2 

Spectral resolution . The natural width of lines in the stratosphere 
varies with pressure. In the lower stratosphere where the Lorentz profile is 
predominant, a typical halfwidth lies between 0.005 and 0.01 curl. In the 
upper stratosphere the Doppler profile is predominant. Since the mqlecylar 
weight of most stratospheric species falls between 20 and 80 then the Doppler 
halfwidths typically lie between 10“3 and 6 x 10“^ cm - * at 1000 cm~l . 

The visible and infrared spectral regions are rich in spectral lines. 

These lines can only be observed with minimum contamination from adjacent lines 
using an instrument with resolution less than the natural halfwidth of the lines 
(i.e., jClO in the upper stratosphere). The only measurements with such a 
resolution are reported by Menzies (MENZIES, 1979; MENZIES et al., 1981) using a 
laser heterodyne spectrometer. The resolution of other instruments is generally 
limited by the technical problems associated with operating them on balloon 
borne platforms. 

Siyce the spectral resolution can vary from 0.02 cm * to more than 
0.1 cm in instruments using absorption techniques and up to several wave- 
numbers in instruments using emission techniques, the observed spectral features 
can be very broad with the peak height reduced in proportion due to the 

convolution of L-iie-ua cui-al— l-iii£=shap£-vzth^thc=ir-str omental— f — T-n— 

general, as the spectral resolution of the instrument decreases, the contamina- 
tion of a spectral line by adjacent lines of the same gas or of other gases 
becomes more of a problem; the inversion process tends to become more prone to 
errors, needs more spectroscopic data and limits the accuracy of the result (see 
section on 'method of analysis'). 

(b) Instruments and Techniques 

The instruments described below are divided into two types classified by 
the relationship between resolution and throughput. 

’ Strong dependence between throughput and resolution . Radiometers and 
classical grating spectrometers with slits are in this class. Improving the 
resolution leads to a decrease in throughput and therefore in luminosity. 

Filter radiometers, using individual interference filters or circular 
variable filters with resolution typically equal to 2% of the center wavelength, 
make measurements integrated over part or whole of a spectral band. In most 
cases, the radiation is modulated by a chopper. Radiometers measuring 
atmospheric emission are often cooled to liquid nitrogen temperature in order to 
achieve the best performance. This technique is only suitable for monitoring 
species such as ozone, water vapour or nitric acid whose spectral bands are 
relatively free from contamination by other species. 

The main advantages of radiometers are that they are small, have few moving 
parts and have a good light gathering efficiency. 

Although Classical spectrometers used on balloon platforms have given good 
results for almost 20 years, we will only illustrate the state of the art by two 
examples . 

Solar absorption spectroscopy is carried out by the group from Liege 
University using the largest grating spectrometer which ever has been flown 
(ZANDER, 1970). 
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It consists of a 38 cm aperture telescope and a 2.5 m focal length spectro- 
meter operating in a double pass configuration. The spectral resolution is 
usually about 0.03 to 0.05 cm” with signal-to-noise ratios in excess of 100 
in the 4000 cm region. The weight of the instrument is 300 kg. High 
quality spectra are obtained with this instrument but the time available for an 
observation limits the number of species that can be detected during the same 
occultation. 

Emission spectroscopy with an Ebert-Fastie type of grating instrument is 
performed by the Upper Atmosphere Group from the National Physical Laboratory, 
England. The resolution at 1000 cm -1 is 0.25 cm - . It covers a spectral 
region from 600 to 2000 cm - * . 

Figure 3 indicates the configuration of the instrument. A telescope 
focuses radiation on to the input slit of the spectrometer; a large concave 
mirror then collimates the radiation onto the diffraction grating. The 
radiation is dispersed into its spectral components by the grating, only a 
narrow spectral interval of radiation is then focused by the concave mirror onto 
the output slit and detector optics. 


Ebert mirror 



Figure 3. The cooled NFL grating spectrometer 
for emission measurements in the stratosphere. 


Three liquid helium cooled copper doped germanium detectors are used and 
permit simultaneous spectra of different species to be measured. As the 
detectors are photon noise limited, all the optical components are cooled to 
liquid nitrogen temperature. The weight of the instrument is 150 kgs. 

Weak dependence between throughput and resolution . The following instru- 
ments employ techniques of spectral filtering which eliminate, to first order, 
the reciprocal relationship between throughput and resolution. These instru- 
ments usually have a greater resolution- luminosity product than classical 
instruments. This can be used to increase the resolution, reduce the sample 
time (and hence allow more species to be observed) or reduce the size of the 
instrument. However, since the incident flux on the detector is considerably 
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greater than in classical spectrometers, these advantages tend to disappear when 
photon limited detectors are used. 

Four examples will be described. The grille spectrometer is a variant of 
the grating instrument, invented and developed by Girard (GIRARD, 1963). The 
entrance slit is replaced by a grille which is an array with hyperbolic 
structure having alternate opaque and transparent bands. The radiation is 
reflected from the grille (Figure 4) then through a conventional Littrow 
mounting grating system and finally passes through the grille. It is then 
filtered to select the order of the grating. The resolution is determined by 
the narrowest step of the grille while the throughput is determined by the total 
area of the grille. 

A 60 cm focal length grille spectrometer using a 30 cm diameter telescope 
has been used to make stratospheric absorption measurements from balloon plat- 
forms. The resolving power is 20000 and the weight of the instrument is 80 kg. 


PRIMARY MIRROR (0 30 cm} 



Figure 4. Optical scheme of the ONERA grille spectrometer for 
solar absorption measurements. M, off axis (F ■= 600 mm); 

Rj , predisper sing grating; Rg , main grating; G, grille; 

D, detector. 


The Fourier transform interferometer (Figure 5) consists of a bean splitter 
(for example, germanium coated KC1 or KBr), a fixed mirror, a moving mirror and 
a detector. Incoming radiation is split into two beans, each bean is reflected 
back to the bean splitter where they recombine and are focused onto the 
detector. The interf erogram, produced at the detector, is obtained by moving 
one of the mirrors and hence changing the path of one bean relative to the 
other. Monitoring the path difference is usually performed by counting the 
fringes when monochromatic radiation from a laser traverses the paths. The 
fringes can be used to trigger the signal electronics to sample the interfero- 
gran. In most modern interferometers "corner-cube" reflectors are used instead 
of plane mirrors. 

The interferometer combines Felgett's Multiplex Advantage (FELGETT, 1958) 
with the Throughput Advantage of Jacquinot (JACQUIHOT, 1958). All spectral 
elements are analysed simultaneously and the resolution is determined by the 
maximum optical path difference and is independent of throughput to first order. 

The main disadvantage of an interferometer is that a spectrum is not avail- 
able in real time, the interf erogram must be transformed by a computer. A badly 
sampled point in the interf ero gran can severely degrade the whole of the trans- 
formed spectrum. An interferometer is sensitive to vibration and changes in its 
operating environment which can displace the mirrors from their true position 
and therefore distort the interf ero gran. 
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Figure 5. The Michelson interferometer: an example 
shown without collimation. S, source; C, chopper; 
BS, beam splitter; Dj D 2 , detectors; R-, 

"roof top" reflectors. The R 2 reflector xs moved 
along the optical axis in order to produce a 
difference of optical path between the two beams. 


When used for limb-scanning absorption spectroscopy the interf erogram must 
be sampled in a short time so as to keep the change in the position of the 
source to a minimum and therefore obtain the maximum possible height resolution 
in the measured concentration profile. In the instrument described by FARMER et 
al. (1980), the interf erogram is sampled in 100 8 which results in a height 
resolution varying from 1 km in the upper stratosphere to 4 km in the lower 
stratosphere. 

Several instruments are now currently used from balloon platforms (see 
SCHINDLER, 1970 and MUR CRAY et al., 1981). 

Gas correlation instruments use a sample of the constituent itself within a 
cell to produce a filter sensitive only to those wavelengths corresponding to 
molecular transitions of the gas. The spectral bandwidth is comparable with the 
line width of the gas, typically 10 -J cm , giving a resolving power of 
10” in the infrared, quite independent of throughput. In practice radiation 
is collected from many lines and therefore the energy falling on the detector is 
considerably greater than in most high resolution instruments. 

One practical example of this technique, the pressure modulator radiometer, 
has been developed for atmospheric remote sensing by Houghton and colleagues at 
Oxford University (TAYLOR et al., 1972). The principle of operation can be 
understood by referring to Figure 6. A cell is placed between the collecting 
and detector optics; the pressure of gas within the cell is modulated by a 
piston at frequency f. If the maximum and minimum pressures within the cell are 
P 2 and P^ respectively, then the transmission through the cell in the region 
of the spectral line at these two pressures varies between t» and which 
leads to a modulated signal at frequency f proportional to the difference bet- 
weei x -and t (Figure 7). Operating the cell at a high mean pressure 
makes the instrument suitable for studying gases with strong pressure broadened 
lines where information concerning the profile is contained within the wings. 
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Figure 6. The pressure modulator radiometer used by the 
Clarendon Laboratory (from HOUGHTON et al., 1971). 
Radiation traverses the cell before being focused onto 
the detector. Pressure modulation is effected by a 
resonant piston driven electromagnetically . 




Figure 7. Illustration of the principle of the PMR 
instrument, a) Spectral line transmission at 
pressures P,, p 2 , p^, b) Cell transmittance. 
Modulation betveen pressures P^ and P2 leads to 
• the F^-p 2 transmittance curve. Modulation 
between pressures P 2 and Pj leads to the P 2 -Pj 
curve. 
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Despite the large resolving power and light collecting efficiency, a 
pressure modulator radiometer can be mechanically simple and robust. The 
problem with such an instrument is in the containment of gases in the cell for 
long periods without deterioration; many gases of atmospheric interest are 
extremely unstable, corrosive or both. 

Laser heterodyne radiometer. Infrared heterodyne detection is achieved by 
allowing radiation from a coherent local oscillator (frequency v^) and from a 
source (frequency v^) to illuminate a square law detector (photomixer). The 
heterodyne signal has a frequency v. - Vj. It is passed through a radio 
frequency filter of bandwidth B then through a second square law detector and 
time averaged (MUMMA et al., 1978). 

Recent developments in high speed infrared detectors and in tunable lasers 
make infrared spectroscopy very attractive to achieve high spectral resolution 
(0.001 cm ) and so completely resolve the profile of a spectral line. 


The heterodyne technique eliminates the effect of detector and background 
radiation noise but at the expense of having to operate with large bandwidths 
and narrow fields of view. 


This technique is limited by the availability of suitable local 
oscillators, so that only selected spectral intervals can be studied. 

Technical difficulties are still limiting balloon borne applications. 
Nevertheless CIO results have been obtained by solar spectroscopy (MENZIES et 
al., 1981) and have recently been reevaluated using improved spectroscopic data 
(MENZIES, 1983), 

I 

METHOD OF ANALYSIS 


When a set of atmospheric spectra is obtained, either by an emission or by 
an absorption technique, the task remains of "inverting" the results to get an 
altitude concentration profile for the species of interest. This means that, 
from equation (6) which expresses radiation as a function of atmospheric state, 
one has to express the atmospheric state in terms of that radiation. 

The solution to such a problem is not unique since in practice there are 
only a finite number of observations and also, because there is noise in all 
physical measurements. Therefore the problem is one of estimation, finding the 
appropriate criteria which determines the best solution from all possible ones 
which are consistent with the observations. The problem can be expresed 
mathematically in the following way; let us assume that the spectroradiometer 
has measured several atmospheric radiances R(i); these may be measurements at 
several zenith angles or at a constant zenith angle and at several altitudes. 

We want to find the concentration profile, C., so that the calculated radiance 
R.(Cj) is in agreement with the measured radiance to within some constant 
factor; mathematically we need to minimize 

F(R - R(C(z)) ) 


subject to certain constraints on C.. 
sum of the squares. ^ 


The function F in a simple form is the 


(a) Calculation of the Radiances 


Practically the calculation of the radiances, R(i), is made by dividing the 
atmosphere into homogeneous layers. 
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For absorption spectra the transmission function t(v) is usually computed 
line by line using the expression; 


T(V) 


exp 


ds £ C. I k. . 
p j J i J-jn 


n, j, and i respectively refer to the number of layers, absorbing species and 
lines of each species, s is the optical length in each layer and k ; is the 
absorption coefficient. 


ljn 


The starting point is the calculation of the monochromatic absorption 
coefficient which requires data on the individual vibration-rotation lines of 
atmospheric molecules. Several groups have collected data for the atmospheric 
species (ROTHMAN et al., 1983; HUSSON et al., 1982). They are very useful 
tools, but work has still to be done in several cases to improve the precision. 


The way the atmosphere is divided into layers influences the precision of 
the results. Thin layers (say 1 km thick) will lead to the best approximation 
but introduce more steps into the calculation. A compromise has to be found 
between precision and calculation time. 


(1) Layers may be simply defined by choosing boundaries corresponding to the 
tangent heights of the observations. 

(2) In the FSCDATM program of APGL (CLOUGH et al., 1982), the selection of the 
boundaries is determined by an acceptable line width differential through 
the layer, the variation of the mean line width staying lower than a 
factor of two from one layer to the next. 

(3) Mankin' s method (MANKIN, 1979) uses an isothermal layer in which both the 
pressure and absorber concentrations vary along the line of sight. 


In all cases, the average pressure and temperature are defined by the 
Curti s-Godson approximation in each layer. 


Pdm/ 


dm 


T dm/ 


dm 


where dm is the differential air mass. 


The calculation of the optical path in each layer must include refraction 
by the atmosphere at large zenith angles (tangent altitude lower than 25 km 
(SNIDER, 1975)). 

Instrumental effects, such as apodization and phase distortion, have to be 
included. Usually the calculated radiances are finally convolved with the 
instrumental function. 


An original algorithm has been developed by MANKIN (1979) for calculating 
transmittances by Fourier transforms of the absorption coefficients which 
greatly reduces the time required to generate synthetic spectra. 

For solar absorption spectroscopy the intensity of the source I Q (v) is, 
in some spectral intervals, severely contaminated by absorption lines in the 
solar spectrum; one has to simulate correctly this absorption (see GOLDMAN et 
al., 1973). Spectra obtained from high altitudes and low zenith angles can be 
used as a reference. To simulate these spectra one single layer of CO at 4500 K 
and a Doppler line shape can be used with good approximation. A sophisticated 
solar model (MULLER and SAUVAL, 1975) can obviously be used. 
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For emission spectra, which are usually measured with low resolution, line 
by line calculations over large spectral intervals can be very costly on 
computer time. So it has proved desirable to develop techniques for averaging 
over the finest structure. A few line spacings is a sufficient range of 
averaging to smooth out lines without destroying the band contour. This leads 
to the idea of a band model. 

In a band model the lines within a band are replaced by an infinite array 
of lines of uniform statistical properties. An interval of this array is taken 
to have similar properties to those of an interval of the band under considera- 
tion. The best known model for dealing with bands of asymmetric top models is 
due to GOODY (1964); the mean transmittance, T, in any spectral interval Av is 
given by, 


s ,, s li \-l/ 2 
exp - — u (1 + - — ) 

d y irp 


where (s/d) is the mean value of the integrated absorption strengths in Av , S, 
divided by the mean line spacing d, Y is the mean width, u the absorber amount 
and p the pressure. The Goody model assumes that the lines are pressure 
broadened, although corrections to the band model for Doppler effects have been 
considered by GILLE and ELLINGSEN (1968). 

To simplify the evaluation of the radiative transfer equation along a path 
through an inhomogeneous atmosphere, the Curtis-Godson approximation is used to 
represent the mean atmosphere. 

From equation (8), 


R(s) 



B 


dr 

ds 


(s.s 1 ) ds'dv 


If the frequency interval, Av = v^ - v is small, then the Planck function 
can be replaced by an average value B(T,s'). Therefore, 


' R(s) 


s 

B(s') 

s' 


V 2 dx ,. 
dl (s * 8 > 

V 1 


ds 1 dv 


where B(s') is the mean Planck function evaluated at the n*"^ 1 
the mean transmittance from the observer up to and including 
The atmospheric radiance can now be calculated rapidly using 
evaluate t . 


step and t is 
the n fc step, 
a band model to 


More sophisticated methods of evaluating the atmospheric radiance have been 
developed, some such as Lowtran (SELBY et al., 1976) use an empirical model 
while others use multiparameter analytical methods to perform the calculation. 


(b) Inversion Process 


The radiative transfer function, deduced from equation (6) by integrating 
over the frequency range v^ to v» is a nonlinear function of the concentra- 
tion profile c(z). To simplify the calculation, the function R(c) can be 
expanded as a Taylor series about a guessed value, c , of the solution: 


R(c) . E <c g) , |& (. - Cg ) 


3 2 R 

3c 2 


(c 




g’ 

2 


2 
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If the second and higher order terns are ignored, then a set of linear equations 
in c remain. Linear matrix techniques can then be used to determine c(z). The 
procedure can be iterated to obtain the optimum values for c(z). 

It is not always necessary to perform a matrix inversion since it is 
possible to solve the equations sequentially as in the "onion peeling" 
technique, which will serve to illustrate a simple inversion algorithm. With 
reference to Figure 8 , the values R correspond to the atmospheric radiances 
measured at zenith angles 6(0 > 90°); the atmosphere is divided into shells 
whose boundaries are defined by the tangent heights at each zenith angle. An 
initial guess is made at the concentration, c, in the outermost shell, S, from 
which R(Cj) is calculated. The calculation is iterated until, 

R - R(Cj) £ e 

where e is a constant comparable with the experimental error in R. Knowing the 
amount in the first shell, the amount in the second shell is determined and so 
on down to the lowest layer. It may not be apparent that any constraint has 
been put on C.; in fact, we have constrained c. by assuming that the values 
are constant Within each shell. Care has to b^ taken when inverting noisy data 
since it often increases the effect of noise by propagating it downwards. For 
example, if there is a positive error in R., then the value c^ will be over- 
estimated resulting in an underestimation in shell l + 1 , irrespective of the 
noise there. 



Figure 8 . Geometry of a limb remote sensing experiment. altitude 

of the tangent ray; R^ Ro , direction of observation; a., , 

zenith angle; S. s. So S/, different layers in the atmosphere 
with respectively Cj C 3 concentration for gases. 


In order to avoid this drawback, Mill has suggested a method for limb 
inversions which is the reverse of the onion peeling method (MILL, 1977), the 
procedures being performed from the lowest to the highest levels. In general, 
since the emission or absorption profiles decrease with altitude the data with 
the highest signal-to-noise ratio is given in preference, making this method 
suitable for the inversion of noisy spectra. 

Some other examples and numerical solutions are given in references GOLDMAN 
and SAUNDERS , 1979; NIPLE et al., 1980; WANG and GOULAG, 1975. 

RESULTS 

Remote sensing of diatomic species, like HC1, HF, CO and NO is usually 
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performed with the absorption technique since the fundamental 1-0 transitions 
are located in the near infrared (X <5.5 microns) where the blackbody emission 
is too small to be detected using current emission techniques (Figure 1). Other 
species have often been measured by both techniques. Only recent results 
obtained by ranote sensing from balloons are discussed. 

(a) HC1 and CIO 

The 1-0 band of HC1 consists of regularly spaced pairs of lines due to 
" > Cl and 'Cl between 2700 and 3000 cm with the center at 2905 cm . 

A complete analysis of this region has been made by FARMER et al. (1976). 

Using the AFGL line listing supported by extensive laboratory measurements on 
CH^, Farmer pointed out that under stratospheric conditions several lines 
appear to be uncontaminated: 

Rj » R 2 » R 3 » p 5» p 6 for r35c1 
R 1 , R 2 , R 3 for H 37 C1 

The R^ and R^ lines of HC1 have been used exclusively where the 
published profiles have been derived from observations of individual lines. The 
PMR technique has been used by EYRE and R0SC0E (1977), to measure the solar 
absorption by all the lines in the band. In all cases the spectroscopic data of 
TOTH et al. (1970) have been used. 

A compilation by RAPER (1977) has discussed the various results, pointing 
out some discrepancies even when measurements were made at the same latitude, 
(See ACKERMAN et al., 1976; WILLIAMS et al., 1976a; EYRE and ROSCOE, 1977.) 

% 

More recent observations are now available. Measuranents of the latitude 
distribution of HC1 have been made by MANKIN and COFFEY (1983) and GIRARD et al. 
(1982) from aircraft flying at 10-12 km. The total column exhibits a minimum 
close to the equator with an increase of approximately 1.5 between 0° and 30° 
and 1.8 between 30° and 60° latitude. Therefore only profiles measured at the 
same latitude can be compared. 

Recent HC1 balloon observations are shown in Figure 9. They all indicate a 
continuous increase from 20 to 40 km with a volume mixing ratio greater than 
10 above 30 km. Although these profiles have been obtained at different 
latitudes (from 30° to 65°) they do not clearly show the latitude variation 
mentioned above for the total column. More precision is still needed to make a 
valuable conclusion. 

The CIO fundamental vibration-rotation band is located near 850 cm ^ . 

The low CIO concentration in the middle stratosphere and the low strength of the 
band make it difficult to measure by current spectroscopic techniques. CIO 
features in the atmospheric spectrum can only be isolated using high resolution 
spectroscopy and consequently the line positions must be known accurately. 

Recent laboratory measuranents by Maki (MAKI et al., 1982) have been made 
with an accuracy of 2 x 10~3 cm - l which have allowed the solar balloon 
heterodyne spectra obtained in 1979 by MENZIES (1983) to be interpreted. The 
band intensity comes from MARGOLIS et al . (1978). Identification of CIO is made 
difficult because the spectral region also contains weak lines of NO- and 
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The observations and altitude profiles (Figure 10) were made during sunset 
and are compared with the results of a photochemical model. 
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Figure 9. HC1 profiles (3.4 v). 



Mixing ratio 


Figure 10. CIO profile obtained at sunset near 11. Op with a 
laser heterodyne spectrometer with 0.001 cm resolution 
(MENZIES, 1983). 
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This result illustrates the possibility of making very high resolution 
spectroscopic measurenents of the atmosphere from balloon borne platforms. 

(b) HF 

The fundamental band of HF is centered at 4005 cm - *. The main absorption 
features in this spectral region come from HjO, however solar CO lines are 
also present. ZANDER et al. (1977) concludes that the R. line at 4038.975 
cm -1 is the most suitable for measuring the HF concentration since it is the 
strongest line at stratospheric temperature not obscured by any other absorption 
feature. Nevertheless, the P. line at 4000.99 cm -i can also be used if the 
instrumental resolution i6 sufficient to isolate the HF line from the water 
vapour and solar CO lines. All the published HF concentration profiles due to 
remote sensing in the IR come from an analysis of either one of these two lines. 
However the lines at 3877.71, 3920.32 and 3833.68 cm * have beeen considered. 

Because of the importance of HF and HC1 in the CFM problem and because the 
HF/HCl ratio provides information on the ratio of natural and anthropogenic 
sources of atmospheric Cl and F, these two gases are often measured 
simultaneously. The latitudinal variation of these two gases exhibits the same 
trend (MANKIN and COFFEY, 1983; GIRARD et al., 1982). 

Zander first detected HF above 27.5 km in 1975; further results have now 
been obtained with the same instrument. The spectral data used by Zander in the 
analysis of the R^ line is due to CONNES et al . (1967). The results are 
presented as a mean mixing ratio above the measurement altitude; they clearly 
show an increase in HF concentration from 1976 to 1979. 

BUIJS et al. (1980), during a flight over Alaska, measured the HF column 
density with a precision of ±25%. The volume mixing ratio profile which they 
obtain has an error of ±50% near 30 km and ±10% below 25 km. 

From a large set of spectra obtained during sunset, FARMER et al. (1980) 
deduce an HF mixing ratio profile showing an increase with altitude with an 
uncertainty of about ±20%. Spectral parameters are due to Goldman (GOLDMAN et 
al., 1974). 

These profiles are shown in Figure 11. Although the 65°^ profile appears 
distinctly higher than the profile measured at 30°N, in good agreement with the 
predicted latitudinal variation, the total column measured by Zander et al. is 
contradictory to this latitude variation. More measurements are still needed to 
draw a conclusion on the variability of the results. 

(c) Freons 11 and 12 

Chlorof luorome thanes have been extensively measured by in-situ sampling 
techniques. These species can be detected in the lower stratosphere by remote 
sensing techniques. However, the concentration profiles which are deduced are 
less accurate than those measured by in-situ techniques due to the lack of 
spectroscopic data. 

Absorption features due to CFM 11 and CFM 12 at 10.8 and 11.8 microns, 
respectively, were first observed in solar occultation measurements from a 
balloon by WILLIAMS et al. (1976b). These features are very weak at strato- 
spheric altitudes due to the rapid decrease of CFM density above the tropopause 
resulting from photolysis. The profiles of CFM 11 and CFM 12 measured by 
Williams and coworkers are shown in Figures 12 and 13. Also shown are profiles 
determined by POLLITT et al. (1979) using emission spectroscopy. 
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(d) NO and N0 2 

Using solar occultation techniques to determine the concentration of such 
short-lived photochemically produced species is open to question. The variation 
of sunlight during the occultation period introduces ambiguities in interpreting 
the measured absorption in terms of meaningful atmospheric abundances. 

Calculation of the photolysis rates for a spherical atmosphere (equinox, 

45° latitude) with a simple scattering approximation has been made by B0UGHNER 
et al. (1980). They conclude that sunset NO measurements below 30 km have to be 
corrected for variations in twilight. However, the correction is highly 
dependent on the photochemical model used. Therefore, remote absorption measure- 
ments of NO and NO^ are often suspect, especially in the lower stratosphere. 

For these reasons, remote sensing measurements of NO are coupled with N0 2 
measurements in order to get a good estimate of the NO + N0 2 profile. 

The 1-0 fundamental band of NO, located at 1876 cm is difficult to 
measure using emission techniques because the energy radiated from this part of 
the spectrum is very weak. However, a PMR instrument using this technique has 
measured NO successfully (DRUMM3ND and JARN0T, 1978; R0SC0E et al., 1981). 
Although the accuracy of these measurements is relatively low they clearly show 
the diurnal variation of NO and NO^. 

Two other sets of NO and N 0 2 results have been obtained using the solar 
occultation technique in the same spectral region and using the same spectro- 
scopic data (APGL, ROTHMAN et al., 1978). 

In preliminary work, MURCRAY et al. (1980) show only NO column densities 
while NIFLE et al. (1980) obtained a profile. A more complete analysis has been 
made by BLATHERWICK et al. (1980) which included a simplified photochemical 
model. The NO and NOj profiles are shown in Figures 14 and 15, respectively. 

The photochemical correction, quite negligible for N0 2 , becomes important for 
NO below 25 km and reaches a factor of 2 around 21 km. 



P0LLIT et *1. (1979) 
Emission Spectrometer 


pi WILLIAMS et al. <1979b) 

^ PMR 

Figure 13. CFC1 2 profiles at 32 °N. 
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From a more recent balloon flight, LOUISNARD et al. (1983) have obtained 
profiles which do not include any photochemical correction. The profiles which 
are shown in Figures 14 and 15 assume that the mixing ratios remain constant as 
the solar zenith angle varies. They observe a discrepancy below 30 km in the NO 
mixing ratio profiles measured at sunrise and sunset while the discrepancy 
appears above 30 km in the NO^ profiles. In both cases the total column is 
lower at sunrise than at sunset in agreement with the total column measurements 
by LAURENT et al. (1981) from a plane. 


Also in the infrared, the weak + v band of NO* has been observed 
by KENDALL and BUIJS (1983) using a high resolution (0.02 cm -1 ) interfero- 
meter. 


NOj has also been measured by the solar occultation technique in the 
visible near 450 nm. In this spectral region numerous Fraunhofer lines appear 
and also corrections for Rayleigh scattering and absorption due to and 
H^O have to be included in the analysis. Profiles have been obtained by 
GOLDMAN et al. (1978), OGAWA et al. (1981) and EVANS et al. (1982). Observa- 
tions of the diurnal variation in the NO2 total column have been made by 
POMMEREAU (1982) but have not been confirmed by other experiments. 

The original technique of NAUDET et al. (1984) uses a star at the limb to 
make atmospheric absorption measurements of N0 2 at night. By using a star as 
the source instead of the sun not only is the spectral contamination due to 
solar lines removed but also errors in the measuranent due to the rapid 
variation of N0 2 at sunrise and sunset. The concentration of NO2 measured 
by Naudet (Figure 15) is clearly greater than the results obtained by other 
workers even though there is a large variation in the measuranents. 

(e) N 2 0 5 

This species has a complex vibration-rotation band centered at 1247 cm'^, 
the density of lines within the band is so great that they cannot be resolved 
and therefore appear as a continuum. Under stratospheric conditions the band is 
severely contaminated by N ? 0 and CH a making spectroscopic detection very 
difficult. 


MURCRAY et al. (1978) has looked for N2O5 in the 8.1 micron region of 
emission spectra without a positive identification and hence gives an upper 
limit of 1.2 x 10^ cm - ^ above 18 km. EVANS et al. (1978) have a tentative 
identification of N2O5 a few hours after sunset; they obtain 2 ppb at 30 km. 
R0SC0E (1982), using a gas cell of CH^ a g a filter, has observed an "N2O5- 
like" substance having the same day-night variation which the models predict for 
N2O5 . 

(f) C10N0- 


Although CIONO2 has absorption bands in the infrared at 780.2, 1292 and 
1730 cm it is very difficult to detect because these bands are contaminated 
by strong absorption features of other stratospheric gases such as CH^ t N 2 0, 

*^2 anc * O3 . The density of lines within the Q-branches is so great that 
even at a resolution of 0.01 cm -1 no structure can be resolved. 


The only observation of C10N0 2 has been made by _ MUR CRAY et al. (1979) 
using an interferometer with a resolution of 0.02 cm -1 . They obtain a mixing 
ratio of 0.9 ± 0.2 ppbv at 30 km although interference by other species makes 
interpretation of their spectral data extremely difficult. 
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(g) N0 3 

The only balloon measurement of NO3 has been made by RIGADD et al . (1983) 
using visible spectroscopy of the 661 nm absorption feature. The measuranent 
was made at night using the occultation technique with a bright star (Venus, 
Arcturus) as a source. From a balloon borne platform at 38.8 km they were able 
to measure, with good accuracy, the profile of NO3 from 30 to 40 km, the 
region where the concentration is a maximum. The profile is shown in Figure 16. 
Ozone was also measured simultaneously to assist in the interpretation of the 
data. 
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Figure 16. NO. profile obtained near 660 nm by 
RIGAUD et al. (1983) 

(h) Nitric Acid 

Nitric acid has strong vibration-rotation bands at 5.9, 7.5 and 11.3 
microns; individual lines within the bands can be resolved only with high 
resolution spectroscopic techniques. GOLDMAN et al. (1971, 1981) have deter- 
mined the statistical band model parameters for the three IR bands and BROCKMAN 
et al. (1978) have measured the positions and strengths of over 1000 lines in a 
part of the 11.3 micron band. However, the lack of spectroscopic data, 
especially the temperature dependence of the absorption coefficients, makes 
analysis very difficult. 

The 11.3 micron band is relatively free from contamination by other species 
and therefore low resolution spectroscopic and radiometric observations in this 
band allows stratospheric nitric acid concentrations to be determined (P0LLITT 
et al., 1979; EVANS et al., 1982; FISCHER et al., 1980). Contamination from the 
wings of the Freon 11 and Freon 12 bands centered at 10.8 and 11.8 microns 
becomes significant in the lower stratosphere. 

The other bands of nitric acid are contaminated by features due to H2O, 

CH^ and NjO. The Q-branch at 1326 cm” 1 can be resolved with medium 
resolution spectrometers and has been used to determine nitric acid concentra- 
tions (L00ISNARD, 1983). 
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Since nitric acid is a relatively long-lived species, dynamics play a 
significant role in its distribution in the lower stratosphere, both with 
latitude and altitude. Latitude surveys of the total column of nitric acid have 
been made from aircraft by MURCRAY et al. (1975), GIRARD et al. (1982) and 
COFFEY et al. (1981). All these results show a minimum in total nitric acide at 
the equator rising by an order of magnitude towards the poles. 


Recent measurements of nitric acid, made both in emission and absorption, 
are shown in Figure 17. Although all but one of these profiles were measured 
at the same latitude it is not possible to determine how much of the scatter is 
due to natural variability rather than measurement error. 
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(i) CO 


The 1-0 fundamental transition of CO is centered at 2143 cm * . Ozone 
(Vf + v-) and nitrous oxide also have strong absorption bands in this 
spectral region; weak absorption features due to water vapour and COj are 
present. However, the greatest difficulty in using absorption spectroscopy to 
measure CO concentration arises from the solar photsphere; a significant 
fraction of the equivalent width of each line is due to solar CO (see 
earlier section) . The telluric absorption by CO can be obtained from spectra 
measured at high altitudes and small zenith angles where the atmosphere mass 
path is very small. 


Several lines (Rq, Rj, and P 2 ) have been analysed by a 
number of experimenters using J spectroscopic data from the AFGL compilation. The 
most recent results are summarized in Figure 18. 
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FARMER'S (1980) result shows a constant volume mixing ratio of 1 x 10 
between 21 and.30 km with uncertainties of ±0.2 x 10~ at the lower altitude 
and ± 0.5 x 10 at the upper bound. 

ZANDER (1981a) using only zenith angles less than 90°, has deduced values 
of the mean mixing ratio above the balloon altitude: 2.1 ± 0.3 x 10 above 

30.6 km and 3.7 ± 0.4 x 10 above 36.8 km. 

From a balloon at 40 km and using the limbscanning technique. LOUISNARD and 
LADO (1983) have measured CO profiles between 20 and 40 km. A minimum in the 00 
mixing ratio is observed around 22-25 km followed by an increase with altitude. 

There is general agreement between these three results within the quoted 
errors. 

(j) Methane and Nitrous Oxide 

-6 -7 

CH^ and N„0, with volume mixing ratios close to 10 and 10 , 

respectively, have been measured extensively by cryogenic sampling techniques 
with good accuracy. Although strong absorption features appear in numerous 
spectral regions the accuracy which can be obtained by renote sensing is limited 
by the uncertainties in the spectroscopic data; the precision of remote sensing 
techniques is still lower than the precision of in-situ techniques. Neverthe- 
less, profiles of CH^ and N.O have best measured using solar occultation 
spectroscopy and are presented in Figures 19 and 20, respectively. 
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(k) Water Vapour 
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Although water vapour is one of the most abundant of minor species in the 
stratosphere it still remains one of the least understood both in terms of 
dynamics and photochemistry. The vibration-rotation spectrum of water vapour is 
rich and complex; its characteristically random appearance being caused by the 
great difference between its three moments of inertia. 

Like other molecules with a complex spectra, the spectral line data for 
water vapour is far from perfect and serious discrepancies can sometimes be 
found between measured and calculated spectra in the regions away from the band 
centers where high rotational energy states occur. Not all water vapour lines 
are suitable for determining concentration profiles since many of the lines are 
saturated even at low concentrations. 

Remote sensing techniques for the measurement of water vapour over long 
atmospheric paths offer a distinct advantage over in-situ techniques since they 
are less likely to suffer from errors caused by local contamination due to water 
vapour outgassing from the instrument and balloon. Local water vapour 
contamination at float altitude can be determined by studying the variation of 
the water vapour signal with zenith angle for small zenith angles. 

Some recent measurements of stratospheric water vapour profiles are shown 
in Figure 21. Most of them have been obtained from measurenents in the Vj 
band centered at 6.29 microns by absorption or emission. Primary interfering 
gases in this region are CH^ and molecular oxygen. Although there is a large 
spread in the profiles they all show an extranely dry lower stratosphere with 
the mixing ratio increasing with altitude in the upper stratosphere. 
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Figure 21. H^O profiles. 




64 


(1) Ozone 

Stratospheric spectra of ozone have been measured in several spectral 
regions. 

Pollitt (POLLITT et al., 1979) has made emission measurements of the v 
and v-i bands at 9.6 microns, close to the maximum of the Planck function at 
stratospheric tanperatures . These two bands have been extensively studied by 
spectroscopists. 

Absorption spectra at sunset have beem^obtained by Farmer (FARMER et al., 
198IJ) in the v x + V 2 + v 3 hands at 2785 cm and the band at 3041 
cm . The absolute accuracy of the 0^ mixing ratio derived from the 
spectra is limited by the uncertainty in the 0^ line strengths for these bands 
(±25%). Obviously this uncertainty affects the total ozone burden but not the 
shape of the profile. 

Solar absorption spectra in the v and bands have been obtained by 
Louisnard (LOUISNARD et al., 1983). They obtain two different profiles from 
analysis of these bands; this discrepancy is due to uncertainties in the line 
strengths used. However, they point out that the uncertainty in the Lorentz 
halfwidth, and especially its variation with tenperature, is also a contributing 
factor . 

Ozone profiles at night (using stellar occultation) have been obtained by 
Rigaud (RIGAUD et al., 1983) from measurasents in the 640-670 nm spectral region 
where the 0- absorption cross section is well known and is independent of 
pressure and temperature; water vapour is only weakly absorbing in this spectral 
region. 

The profiles obtained by remote sensing techniques are shown in Figure 22 
and are in good agreement with those obtained using in-situ techniques. 

CONCLUSION 

This review of recent stratospheric measurements obtained by remote sensing 
techniques in the visible and infrared illustrates how the measuring techniques 
have developed and brings out certain points which are summarized below. 

Comparison of results obtained by different experimenters shows that it is 
very important to be able to distinguish between the natural variability in 
concentration and the error in the measurement. Profiles measured with greater 
accuracy are required to validate photochemical models of the stratosphere. In 
a recent paper, HARRIES (1982) suggests that key stratospheric species need to 
be measured to an absolute accuracy of 2%; even though such an accuracy has not 
been achieved the trend is in this direction. 

International intercomparison campaigns (CNES, NASA, MAP/GLOBUS) have 
brought together identical and complementary measurement techniques with the aim 
of making simultaneous measurements of the same species within the same air 
mass. The results from these campaigns, still being analyzed in most cases, do 
not appear in this review. These campaigns have encouraged experimenters to pay 
rigorous attention to instrumental calibration and data analysis in order to 
minimize all sources of error. 

In parallel, the work of laboratory spectroscopists has brought about 
improvements in the spectral parameters of key stratospheric species. However, 
many of the heavy molecules (HNO^ , N^O., I ' reons ) have received less 
attention and the lack of spectral data limits the accuracy with which these can 
be measured. 
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The simultaneous measurement of several species with identical or 
complementary instruments has proved to be indispensible for obtaining 
cross-correlated sets of data on key species and for studying their natural 
variations. A Fourier Transform interferometer, being able to multiplex 
spectral information would appear more suitable for making simultaneous 
measurements of groups of species than instruments which carry out spectral 
measurements line by line. However, grating spectrometers can reduce this 
advantage considerably by scanning the grating rapidly and by using different 
orders to analyze several spectral regions in parallel. 

The development of remote sensing instruments must continue in two 
directions to allow new species to be detected and to enable measurements to be 
made with greater accuracy. 

(a) Resolution 

The report of the detection of CIO by 1R heterodyne absorption spectroscopy 
illustrates the possibility of making atmospheric spectral measurenents at very 
high resolution comparable with the Doppler line widths. At present measure- 
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merits are limited to the 8-12 micron region where the technology of CO. and 
N^O laser local oscillators is well developed. However, by using tunable 
diode lasers as a source of local oscillator power it will become possible to 
make very high resolution measuranents anywhere in the infrared spectral region. 

Although laboratory interferometers have been developed for obtaining high 
resolution spectra (0.0025 cm 1 2 3 4 ) it appears unlikely that balloon borne 
infrared interferometers will be able to match the resolution which heterodyne 
systems are potentially capable of achieving. 

All the high resolution instruments that have been reported use absorption 
spectroscopy ; there exists a requirement for the development of a high resolu- 
tion mid infrared mission instrument for studying the diurnal variation of 
photochemical ly active species. 

(b) Signal-to-Noise 

The progress in the field of infrared detector technology promises a real 
improvement in the signal-to-noise ratio and therefore an improvement in the 
accuracy of measurements. Photon limited detectors are most suited for use in 
high resolution cooled grating spectrometers where the photon flux on the 
detector is small; their full potential is not realised in interferometers 
unless the spectral bandwidth (and hence multiplex advantage) is severely 
reduced by using cooled narrowband filters. 

Although there will be an increasing shift towards remote sensing of 
stratospheric species from satellites, remote sensing from balloons will still 
have an important role to play in: 

(1) detecting new species, 

(2) addressing specific points of stratospheric composition, 

(3) testing new instrumentation, and 

(4) validating satellite measurements. 
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INTRODUCTION 

The history of submil limetre atmospheric research began in 1957 with the 
publication by GEBBIE (1957) of the first measurements of atmospheric trans- 
mission, at wave numbers between 10 and 35 cm at a spectral resolution of 
0.2 cin . These data were recorded using a large aperture (300-mm) lamellar 
grating interferometer and a Golay cell detector at the Jungfraujoch mountain 
observatory in Switzerland. The spectrum observed is reproduced in Figure 1 
and shows, for the first time, the existence of several transmission windows 
between 12 and 28 cm . Many of the observed spectral features were 
identified as due to water vapour H 2 0 and others were provisionally assigned 
to ozone 0, . Subsequent laboratory spectroscopy (STONE, 1964; GEBBIE et al., 
1966) confirmed that numerous 0^ lines were present in the observations, along 
with magnetic dipole transitions in molecular oxygen (GEBBIE et al., 1969). 

The development of submil limetre atmospheric research from 1957 until about 
the end of 1976 has been reviewed by HARRIES (1977). Following the early 
activities just described, a number of groups carried out studies of the sub- 
millimetre emission spectrum of the lower atmosphere from aircraft (HARRIES et 
al., 1972; BUSS0LETTI and BALUTEAU, 1974; MARTEN and CHAUVEL, 1975; MANKIN, 
1975). During these studies, new methods such as phase modulation (CHAMBERLAIN, 
1971, 1979) and new heljum cooled detectors were introduced, and spectral 
resolutions of 0.06 cm were achieved. Figure 2 illustrates two strato- 
spheric emission spectra in the 12-27 cm range recorded from an' aircraft at 
an altitude of 12 km, a zenith angle of 75° and with a spectral resolution of 
0.06 cm . These observations showed the existence of highly reproducible 
spectral structure due to Oj and other molecules. 

Experiments following that were mounted on balloon platforms (HARRIES et 
al., 1973, 1976; CLARK and K0JDALL, 1976). At the higher altitudes accessible 
by such methods, lines were narrow (particularly H.0, which causes consider- 
able opacity still at aircraft altitudes), and it became possible to detect weak 
lines of other species and to use limb-sounding methods to derive concentration- 
height profiles throughout the stratosphere for a number of species. Figure 3 
shows an example of an emission spectrum measured from a balloon floating at 35 
km with a zenith angle of 92° (HARRIES et al., 1976). The molecular assignments 
made at the time for many of the emission lines are indicated in the figure. 

Considerable progress was made, also, in the early 1970's in the avail- 
ability of laboratory spectra of stratospheric trace gases, mainly through the 
work of Fleming and co-workers (e.g., FLEMING and WAYNE, 1975; FLEMING, 

1976a, b). Figure 4 shows a beautiful result from FLEMING and WAYNE (1975): a 

laboratory measurement of the transmission spectrum of 18 torr of ozone in a 1-m 
cell (Figure 4b) in comparison with a calculated spectrum (Figure 4a), showing 
excellent agreement between the two. 
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Figure 1. The solar absorption spectrum in the submil lime tre 
region measured using a lamellar grating interferometer at 
an altitude of 3.58 km, a temperature of 2°C, a pressure of 
503 torr, precipitable water of 2+0.5 mm, and a mean 
zenith angle of 49.5°. Recording time 35 min. This 
represents the first observation of the transmission 
windows at 12, 14, 16, 11 and 28 cm"!. From GEBBIE (1975). 



Figure 2. Stratospheric emission spectra in the 12 - 27 cm"l 
region measured from an aircraft at an altitude of 12 km, a 
zenith angle of 75°, and a spectral resolution of 0.06 cm"l . 
From HARRIES et al. (1972). 
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Figure 3. Stratospheric emission spectra in the 34 - 39 cm -1 
region measured from a balloon at an altitude of 35 km, a 
zenith angle of 92 and a spectral resolution of 0.07 cm“l. 
From HARRIES et al. (1976). 


Since the time of the HARRIES (1977) review, progress has occurred both in 
the direction of considerably improved experimental measurements from balloons 
and in the area of theoretical studies of the far infrared spectrum. The first 
aspect will be discussed in the section "Recent Experimental Results" and the 
second in "Theoretical Calculations and Spectroscopic Data". In the light of 
the recent results, out present understanding of the far infrared and submilli- 
metre spectrum is reviewed in the last section. 

RECENT EXPERIMENTAL RESULTS 

Kendall and Clark of the University of Calgary have reported further 
measurements of the stratospheric emission spectrum from balloon altitudes, 
using a rapid-scan Michelson interferometer and a cooled detector (KENDALL and 
CLARK, 1978, 1979a, b, 198J; KENDALL, 1979). These data were measured at a 
resolution up to 0.07 cm over the 30-110 cm range. A number of flights 
were made, details of which may be found in KENDALL (1979). 

In this work the most notable result has been the possible detection of 
OH. The question of OH emission has been investigated in some detail by these 
authors (KQ7DALL, 1979; KBIDALL and CLARK, 1979b), including a calculation of 
the line positions and strengths, and it appears that the spectral resolution of 
their measurements was just on the limit for the detection of OH. Figure 5 
shows an example of an emission spectrum from KENDALL and CLARK (1979b): The 

top curve is the experimental result (see caption for details), and the lower 
curve is a line-by-line calculation including HjO and 0^ only. The arrows 
mark the predicted OH line positions, and it can be seen that features have been 
detected that might be due to OH. 
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Figure 4. Comparison of synthetic (a) and experimental (b) 
transmittance spectra of ozone for laboratory conditions 
(room temperature approximately 2 296 K; path length 1.0 m; 
pressure^l8 torr 0 J + 12 torr 0 ; spectral resolution 
0.05 cm ). The arrows indicate the positions of strong 
H„0 lines, which influence the experimental spectrum. 

From FLEMING and WAYNE (1975). 


GROSSMAN and OFFERMAN (1978) reported rocket-bo^ne measurements of the 
63 pm emission from thermospheric atomic oxygen 01 ( p.- p 2 fine 
structure transition in the ground state). An interesting interpretation of 
the data is given by these authors in terms of the breakdown of local thermo- 
dynamic equilibrium at the height of the measurements (about 80 km), which 
results in much lower emission intensity than previously assumed. The 
significance of this result, of course, is that this emission line represents a 
much less effective radiative cooling mechanism than previous calculations 
(assuming LTE) had indicated. 

CARLI et al. (1980, 1982, 1983) have recently measured the emission 
spectrum of the stratosg^ere in the 7-85 cm interval with an unapodized 
resolution of 0.0033 cm . This significant improvement in resolution, made 
possible by the use of both a high efficiency polarizing Michelson interfero- 
meter (MARTIN and PUPLETT, 1969) and new sensitive bolometers, has provided very 
detailed information on the complex structure of the stratospheric submillimetre 
spectrum. Some of the difficulties encountered in the early assignments have 
been explained and several new constituents have been identified. 

New measurements at higher frequencies have recently ^een obtained by 
Chance and Traub. These measurements cover the 70-260 cm interval with an 
unapodized resolution of 0.032 cm and have made possible the study of a new 
spectral interval and the identification of several lines due to important 
stratospheric species such as HF, HC1 and OH (CHANCE et al., 1980; TRAUB and 
CHANCE, 1981; CHANCE and TRAUB, 1984b). 
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Figure 5. Measured submillimetre stratospheric emission 
spectrum obtained using a Michelson_j.nterf erometer 
at a spectral resolution of 0.08 cm from a 
balloon at approximately 35 km altitude (top curve) 
compared with a calculated spectrum for a standard 
atmosphere of H^O and 0 on which the positions 
of OH lines have been marked by arrows; an 0 ^ 
line is similarly indicated. From KENDALL and 
CLARK (1979b). 


The results of the recent measurements by Carli et al. and by Chance and 
Traub together provide a full coverage at high resolution of the emission 
spectrum of the stratosphere in the far infrared and submillimetre region, and 
will be reviewed in the last section. 


THEORETICAL CALCULATIONS AND SPECTROSCOPIC DATA 

In France, Rabache and Rebours continued their studies of the submillimetre 
spectrum of the lower stratosphere. RABACHE (1977) studied the details of the 
pure rotation spectrum of nitric oxide, NO, and provided a listing of line 
positions and strengths. He concluded, however, that the NO lines made such a 
weak contribution to the spectrum by contrast with lines of 0_ , H.O and 0 2 
that they were essentially undetectable, given the current instrumental per- 
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formances at the time. REBOURS and RABACBE (1979) discussed the method proposed 
by BURROUGHS and HARRIES (1970) for making quantitative estimates of strato- 
spheric ^0 concentration using an "internal calibration" procedure against 
0 . lines. Rebours and Rabache were able to show that the method gave reason- 
ably accurate results (±20%) in the case of HjO which generally decreases in 
density with height in the stratosphere in a monotonic way. However, much 
larger errors (i 100% and more) were found in the case of 0,, which is 
distributed in the atmosphere in the form of a layer: In this case, clearly 

the method breaks down and cannot be used. 


At the National Physical Laboratory, UK, BANGHAM (1978a, b) carried out a 
number of line-by-line spectral emission computations, generally for the case of 
limb sounding geometry with an observer on a balloon or spacecraft. These cal- 
culations were carried out in order to deduce the amount of retrievable 
information contained in the submillimetre spectrum at a resolution of about 
0.01 cm (the best resolution achievable at this time) as well as to analyse 
existing da£a. An example of this work is given in Figure 6, where the spectrum 
near 41 cm is calculated in the case of a high resolution limb-scanning 
submillimetre spectrometer mounted on a spacecraft at 250 km looking at the 
atmospheric limb at a tangent height of 20 km. The solid line is for H^O, 

O2 and Oj only: The broken line shows t^e effect of adding HC1 and HF with 

a relative concentration of 1 part in 10 7 by volume (constant with height). 

This result should be compared with the more recent observational results shown 
in Figure 12. 


Several laboratory measuranents and calculations have been made in order to 
improve the spectroscopic data available in the far infrared. Useful new 
information has been included in the revision of 1982 of the AFGL Atmospheric 
Absorption Line Parameter Compilation (ROTHMAN et al., 1983b) and in the new 
AFGL Trace Gas Compilation that also appeared in 1982 (ROTHMAN et al., 1983a). 
Furthermore, a new catalogue hjis been started at JPL (POINTER and PICKETT, 1980) 
for frequencies up to 100 cm , on the basis of new microwave measurements and 
theoretical calculations. This catalogue, of which two revisions have already 
been produced in 1981 and 1984, now provides a rather exhaustive data base as 
far as constituents, line position and line strength are concerned. 



Figure 6. A calculated spectrum of the submillimetre 
atmospheric emission spectrum in the 42 - 42 cm 
region for a standard atmosphere containing ^0, 

Oj , and 0j, plus HC1 and HF completely mixed at 
a mole fraction of 10 . The calculation was 

undertaken for the same observational geometry and 
spectral resolution as that in Figure 5. From 
BANGHAM (1978a). 
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High resolution measurements in the submil lime tre region have recently been 
carried out at the Instituto Di Ricerca Sulle Onde Elettromagnetiche. The pure 
rotational spectra of several minor atmospheric constituent^, including 0^ , 

SOj , Hj, NOj, and H-Oj, have been measured in the 10-80 cm 
spectral region with a frequency accuracy of the order of 2.10-4 cm 
(BALDECCHX et al., 1984). A full analysis ha6 been completed in the case of 
SO 2 (CARLOTTI et al., 1984a) and ozone (CARL0TT1 et al., 1984b) and is in 
progress for the others. An example of the results obtained is shown in Figure 
7 in the case of nitrogen dioxide (CARLI et al., 1984a). The high resolution 
makes it possible to resolve most of the transitions and provides important 
information on the extension of the microwave data to higher frequencies. 

In the 100-250 cm * spectral interval laboratory measurements of C10N0„ 
(CHANCE and TRAUB, 1982) and H0C1 (CHANCE and TRACB, 1983) have been made at the 
Smithsonian Astrophysical Observatory. These molecules play a major role in the 
ozone chemistry, but are difficult to measure in the stratosphere. 

While these measurements and calculations are rapidly filling the gaps as 
far as the knowledge of line position is concerned, large uncertainties still 
exist in the value of the half width of most transitions of the stratospheric 
spectrum and new measurements of this quantity will be required in the future 
for a more precise data analysis. 

THE FAR INFRARED STRATOSPHERIC SPECTRUM 

(a) The Main Features 

A pictorial representation of the far infrared spectrum of a selection of 
stratospheric constituents is shown in Figure 8 (PARK et al., 1981). This 
figure is not updated and in fact some constituents like SO 2 and HBr have not 
yet been detected in the far infrared stratospheric spectrum while some more 
recent measurements like HCN and OH are not included. Nevertheless, the figure 
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Figure 7. Laboratory mjasuraaents of the NO 2 absorption spectrum^ 
between 37 and 39 cm . The unapodized resolution of 0.003 cm 
makes it possible to resolve most of the individual rotational 
transitions within this Q-branch. 




Figure 8. The logarithm of the line strength at 296 K as a function of the 
frequency for a few atmospheric gases, from PARK et al. (1981). This dia- 
gram provides a useful pictorial representation of the position and strength 
of the transitions that contribute to the stratospheric emission spectrum 
in the 0 - 250 cm-1 . 





Figure 8 continued 
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Figure 8 continued. 

provides some very useful information on the position and the density of the 
lines that compose the stratospheric spectrum. 

Magnetic dipole transitions of molecular oxygen and pure rotational 
transitions of both water vapour and ozone are the dominant features because of 
the relatively high concentration of these constituents. Oxygen is present in 
the entire spectral interval with a few strong and easily resolved lines that 
are grouped in a typical sequence of triplets. At high resolution also, weaker 
transitions due to the isotropic form "O' 16 © become visible (CARLI et al., 
1980). 


Water vapour also is present in the entire spectral interval with several 
features irregularly spaced. The strong water vapour lines can grow into broad 
features and eventually make the full spectral interval opaque when a 
sufficiently high column density is present along t&e observation path. At high 
resolution several sharp transitions due to HDO, °0, H. 0 and water 
vapour excited to the first level of the v» vibrational state become visible 
(CARLI et al., 1980). 

Ozone is present with a sequence of numerous closely spaced transitions 
that extend from the submillimetre region up to about 150 cm . A decreasing 
density of ozone lines can be noticed in the 100-150 cm region. The Q- 
branches of ozone can also be easily detected in very low resolution spectra, 
but a relatively high resolution is needed to detect individual transitions. 

With a resolution of 0.003 cm , most features of the ozone spectrum are 
resolved and several transparent "windows", where transitions of other minor 
constituents can be studied, become visible. Transitions of rotational spectra 
of ozone isotopes 6 0 i6 0 l6 0 and 6 0* 6 0 l6 0 and of vibrationally 
excited ozone to Jtfie first Vj, v' 2 and levels have also been identified 
in the 10-100 cm region, stressing the dominant role that ozone plays in the 
submillimetre stratospheric spectrum (CARLI et al., 1980, 1982). 

(b) Hitric Acid, HNOj 

Nitric acid was first identified in submillimetre spectra measured with a 
resolution of 0.06 cm (HARRIES, 1973) and its detection has been confirmed 
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by measurements at higher resolution. Two features occur in the spectral 
interval shown in Figure 9, but are not sufficiently strong for a positive 
detection, because of the geometry of observation (tangent height of about 28 
km). At lower altitudes the concentration of HNO. increases and several lines 
can be detected in spectra with a tangent height of 16 km (see Figure 10) in the 
spectral interval from 9 cm to about 24 cm . The detection becomes 
difficult at higher frequencies because of both the decreasing intensity of 
HNO^ features and the increasing opacity of the spectrum. 

(c) Nitrous Oxide, N 2 0 

Nitrous oxide is a linear molecule and its equispaced lines have a period 
of about 0.84 cm" . The detection of this molecule was already possible in 
spectra at 0.06 cm resolution (HARRIES, 1973). With a resolution of 0.003 
cm more than 16 lines of nitrous. oxide can be resolved going from the , 

J" ■ 11 transition at 10.053 cm to the J" = 41 transitions at 35.145 cm 
One example is shown in Figure 10. 

(d) Nitric Oxide, NO 

The concentration of nitric oxide undergoes a diurnal cycle so that its 
monitoring throughout a 24 hour period with an emission measurement is important 
for the study of the photochemistry of this constituent. As pointed out by 
RABACHE (1977) nitric oxide shows only some rather weak features in the sub- 
millimetre stratospheric spectrum, but with the sensitivity of more recent 
measurements some small signal has been detected (CARLI et al., 1982). The 
strongest and better resolved feature is the “9.5 transition at 

36.077 cm” as shown in the stratospheric spectrum of Figure 11. 


H,0.0j 



Figure 9. Stratospheric emission spectrum with an 
unapodized resolution of 0.003 cm measured from 
a floating altitude of 38 km with a tangent height 
of 28 km. In this interval the two features of hydro- 
chloric acid are easily resolved, while the emission 
features of nitrous oxide are blended and the emission 
of nitric acid is insufficient for a positive detection. 
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Figure 10. Stratospheric emission spectra measured with an 
unapodized resolution of 0.003 cm from a floating 
altitude of 38 km in the case of a tangent height of 16 km 
(curve (a)), and the the case of a tangent height of 28 km 
(curve (b)). In this interval two features of nitrous 
oxide are resolved and a strong signal due to nitric acid 
becomes visible at low tangent heights. 
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Figure 11. Stratospheric emission spectrum measured with an 
unapodized resolution of 0.003 cm from a floating 
altitude of 38 km with a tangent height of 28 km. In this 
interval a small feature of nitric oxide has been detected. 
The negative signal at_phe centre of the strong water 
vapour line at 36.6 cm is an instrumental artifact due 
both to the use of the zenithal signal as a reference zero 
level and to water absorption within the instrument. 
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(e) Nitrogen Dioxide, N0 2 

As in the case o£ NO, the concentration of stratospheric nitrogen dioxide 
shows some diurnal variation and the detection of this constituent with an 
emission technique is of great interest. From recent high resolution laboratory 
measurements (CARLI et al., 1984), an emissivity of about 3% is expected from 
stratospheric N0 2 both in the interval ^etween 52.657 cm -1 and 52.665 from 
the tjranch and in the 0.003 cm spectral element around 67.345 cm~ 

from the Q, branch. Much lower emissivity i6 expected from the other Q- 
branches with the present resolution. The first feature is unfortunately 
covered by an ozone Q-branch, but the second occurs in a relatively open 
spectral interval. A positive assignment is made difficult by the weakness of 
the signal and by the blending with other small features at 67.3518 and at 
67.3525 cm due to ozone and isotopic ozone respectively, but some evidence 
exists that a signal from stratospheric N0 2 has been detected, see Figure 14 
(CARLI et al., 1983). It is clear that the measurement of stratospheric 
nitrogen dioxide in the far infrared region requires a sensitivity and a 
reso lution^higher than those used up to now and it may be possible only in the 
67.345 cm region. 

(f) Hydrofluoric Acid, HF 

The first transition of the rotational. spectrum of HF is clearly detected 
in the stratospheric spectrum at 41.111 cm , see Figure 12. A small feature 
at 41.114 cm - due to ozone excited to the first v 2 vibrational level and 
the tail of the strong water vapour line at 40.988 cm are superimposed on 
the HF line, but are not sufficient to mask the contribution that clearly 
anerges above the other signals (CARLI et al., 1983). The same does not occur 
in the case of the next two transitions at 82.171 and 123.130 cm^ that are 
covered by strong lines of water vapour at 82.155 and 123.129 cm A , 
respectively, and cannot be resolved. At higher frequencies, three other HF 
features, that are free from blending effects with strong water vapour lines, 
have been measured by TRAUB and CHANCE (1981). 
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Figure 12. Stratospheric emission spectrum measured with an 
unapodized resolution of 0.003 cm from a floating 
altitude of 38 km with a tangent height of 28 km. In this 
interval HF, HC1 and HCN display strong features that can 
be resolved and used for their measurement. 
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(g) Hydrochloric Acid, HC1 

35 37 

The two isotopes of chlorine. Cl and Cl, occur naturally with an 
abundance ratio of about 3 to 1 so the pure rotational spectrum of HC1 displays 
two adjacent sequences of equispaced lines with a perio^ of about 20.8 cm . 

Some of the features that can be measured below 100 cm are shown in Figures 
9, 12 and 15 (CARLI et al., 1980, 1983). Pressure broadening in measurements 
with a low tangent height introduces some slight blending between the J" = 

0 H 3 'C1 transition at 20.874 cm and the ozone line at 20.860 cm (see 
Figure 9), and a serious blending betwem the J" ° 1 H J/ C1 transition at 
41.681 and the ozone line at 41.675 cm (see Figure 12). No blending has 
been observed for the J" = 0 and the J" = 1 transitions of the other isotope and 
for the J" = 2 and J" = 3 transitions of both isotopes. 


Several othex transitions are present about 100. cm In particular, the 
J" ° 4, 5, 6 of H 35 C1 at 104.139, 124.827, 145.440 cm -1 and the J" = 6 of 
H J 'C1 at 145.223 cm have been successfully used by TRADB and CHANCE (1981) 
to measure a profile of HC1 in the stratosphere (see Figure 16). 


(h) Hydrogen Cyanide, HCN 


Hydrogen cyanide is a linear molecule with a large dipole moment which 
displays in the submillimetre a strong rotational spectrum. About 12 
transitions of this molecule can be resolved in the stratospheric spectrum going 
from the J" = 3 transitions at 11.825 cm up to at least the J" ° 21 transi- 
tions at 64.918 cm (CARLI et al., 1982). An example of an HCN line that has 
been resolved in the stratospheric spectrum is shown in Figure 12. 


(i) Carbon Monoxide, CO 


The rotational spectrum of carbon monojp.de displays a sequence of equi- 
spaced lines yith a period of about 3.8 cm and a maximum line strength 
around 50 cm . Among about half-a-dozen lines that have been identified in 
the stratospheric spectrum (CARLI et al., 1980), the J" «* 12 transition at 
49.932 cm (see Figure 13) provides the strongest and best resolved signal of 
CO in the far infrared stratospheric spectrum. 

(j) Hydroxyl Radical, OH 


The measurement of hydroxyl radical is of great importance for the study of 
the ozone chemistry because of the role of this molecule in many stratospheric 
reactions. The possibility of measuring this constituent in the far infrared, 
first suggested by KSJDALL and CLARK (1979b) (see Figure 6), was confirmed by 
the recent high resolution measurements. 

The far infrared spectrum of OH starts at 60 cm * and h88 its strongest 
features in the 150 cm region. The first four weak features, respectively 
at 61.200, 61.303, 83.724 and 83.868 cm , have been clearly identified in 
spectra at high resolution (CARLI et al., 1983). An example of these results is 
shown in Figure 15. The second and fourth transitions are seen as isolated 
lines down to tangent heights lower than about 33 km, as shown in Figure 15. 


Furthermore, a total of twelve transitions of OH have been identified in 
the 80-200 cm in stratospheric spectra recorded by CHANCE and TRAUB (1984b). 
An example of these measurements is shown in Figure 17. 


(k) Other Constituents 


The far infrared spectral region continues to be analyzed in order to 
identify other constituents that can usefully be measured in these regions 
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Figure 13. Stratospheric emission spectrum measured with an 
unapodized resolution of 0.003 cm from a floating 
altitude of 38 km with a tangent height of 28 km. In this 
interval a small feature due to carbon monoxide is clearly 
visible close to one of the few transparent spectral 
intervals that are present in the submil limetre stratospheric 
spectrum. 
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Figure 14. Stratospheric emission spectrum measured with an 
unapodized resolution of 0.003 cm from a floating 
altitude of 38 km with a tangent height of 35 km. A small 
feature, that is tentatively assigned to nitrogen dioxide, 
suggests the possibility of the measurement of this con- 
stituent, given a larger signal-to-noise ratio. 



86 



WAVENUMBER CM-1 


Figure 15. Stratospheric emission spectrum measured with an 
unapodized resolution of 0.003 cm from a floating 
altitude of 38 km with a tangent height of 28 km. In this 
spectral interval, features due to hydrochloric acid and 
hydroxyl radical are clearly resolved. 


Some new results have been obtained very recently. 

First, on the basis of an expected feature of hydrogen peroxide at 112 
cm a firm 2 a upper limit of 0.05 ppb was determined for HjOj strato- 
spheric mixing ratio near 26.5 km by CHANCE and TRAUB (1984a). 

Second, a feature of CIO has been identified at 22.89 cm * (CARLI et al., 
1984b) in recent high signal-to-noise ratio stratospheric spectra recorded with 
J He cooled detectors. 

CON CLUS IONS 

The understanding of the far infrared stratospheric spectrum has improved 
continuously since the early measurenents by H. A. Gebbie. A rather exhaustive 
picture of the constituents that can be detected and of the spectral intervals 
that provide the most useful measuranent is now available. 

The results obtained show that far infrared ranote sounding can provide 
valuable and unique information on stratospheric composition and its evolution. 
Further improvements are expected in the case of a few molecules from the 
measurement of better line parameters (half width in particular), in the case of 
lines that are not fully resolved even at higher resolution and, probably most 
important of all, from measuranents with better signal-to-noise ratio, as is 
possible with high sensitivity detectors associated to cooled narrow band 
f il ter s. 
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Figure 16. Examples of lines of hydroxyl radical that have 
been detected in the 80 - 220 cm 1 interval in the strato- 
spheric emission spectrum obtained by CHANCE and TRAUB during 
the balloon flight of 23 January 1983. The abscissa scale 
is in cm 






88 



Figure 17. Example of lines of hydrochloric acid that have 
been detected m the far infrared stratospheric spectrum 
obtained by CHANCE and TRAUB during a balloon flight on 
20 June 1984. The dotted line is a nonlinear least- 
square fit to the observed spectrum (solid line). 
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ABSTRACT 

The Balloon-borne Microwave Limb Sounder (BMLS) measures atmospheric 
thermal emission from millimeter wavelength spectral lines to determine vertical 
profiles of stratospheric species. The instrument flown to date operates at 
205 GHz to measure C£0, 0^ , and ^2°2' A 63 GHz radiometer will be added to 
test the technique for determining tangent point pressure from the MLS experi- 
ment on the Upper Atmosphere Research Satellite (UARS). Many additional species 
could also be measured by the BMLS. A radiometer, at .270 GHz would provide 
measurements of H0_, NO., HNO^ , N.O, °0 0 6 0, and HCN. With this 

addition the BMLS can test the current theory of 0^ photochemical balance in 
the upper stratosphere. 


INTRODUCTION 


As part of NASA's Upper Atmosphere Research Program, a balloon-borne micro- 
wave limb sounder (BMLS) has been developed to perform measurements of the 
earth's stratosphere. The BMLS followed earlier development of an aircraft 
instrument for the same purpose (WATERS et al., 1979, 1980) and was a step 
towards development now underway of a similar instrument (WATERS et al., 1978) 
for the Upper Atmosphere Research Satellite (UARS) planned by NASA for launch in 
1989. 


The BMLS is the first use of microwave remote sensing techniques for limb 
sounding. Good reviews of the principles and application of microwave remote 
sensing are by STAELIN. (1969), SCHANDA (1976), and NJ0KU (1982). The major 
advantages of limb sounding over other viewing geometries are that (1) it 
provides the largest possible signals from species with small abundances, and 
(2) it has much better vertical resolution than ground-based or aircraft-based 
measurements which obtain vertical resolution from spectral line shape. 
Technology is now available for balloon-based instruments operating at 
frequencies below approximately 300 GHz (wavelengths longer than 1 millimeter). 
Advances in technology are currently being made which can provide instrumenta- 
tion to much higher frequencies. Figure 1 illustrates some of the many species 
and spectral lines which are available for measurement in the 200-300 GHz 
frequency range alone. The species abundances which were used for the calcula- 
tions in Figure 1 are given in Table 1. Spectroscopic data used for these 
calculations are from the JPL catalogue (P0YNTER and PICKETT, 1981) with the 
latest revisions as of March 1984. Figure 2 compares the vertical resolution of 
limb sounding and ground-based or aircraft-based microwave measurements. Other 
general advantages of microwave sounding include (1) the composition measure- 
ments are very much less sensitive to uncertainties in the atmospheric 
temperature profile than are infrared emission measurements, and (2) continuous 
measurements can be made during both day and night. The major disadvantage of 
limb sounding is that the instrument must be above the atmospheric layer being 
measured to fully exploit the technique. 

Development of the BMLS was motivated by the need for measurements of 
stratospheric C$,0 which is one of the best indicators of stratospheric ozone 
depletion by industrial chlorof luocarbons (CFCs) (e.g., MOLINA and ROWLAND, 

1974; STOLARSKI and CICERONE, 1974; CRUTZEN et al., 1978; WUEBBLES, 1983). The 
CIO rotational line at 204.352 GHz was selected for measurement because of its 
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Figure 1. Atmospheric spectral lines in the 200-300 GHz range. 
The vertical scale is the logarithm of the calculated optical 
depth for a limb sound path through the middle stratosphere. 
Species abundances used for the calculation are in Table 1. 
Within each 0.5 GHz interval all spectral lines from a given 
species are summed. The heavy horizontal line is the ap- 
proximate noise level for anission measurements using a state- 
of-the-art radiometer with a few minutes integration and 
spectral resolution to resolve the stratospheric lines. 


freedom from interference by other atmospheric lines and to optimize measurement 
signal-to-noi se . (Line strengths increase with frequency whereas instrument 
sensitivity decreases.) An 0^ spectral line at 206.132 GHz and an H20„ 
line at 204.575 GHz were also selected for simultaneous measurement. The 0^ 
line was chosen because of the insensitivity of its anission to stratospheric 
temperature variations. A stratospheric temperature variation of 10°C changes 
the 206.132 GHz 0^ limb emission by only one per cent (WATERS, 1980). The 
H_0„ line was chosen because it could be measured by the same radiometer 
which measures Ci.0 and 
WATERS et al. (1981). 


. Initial measurement results are described by 



92 


Table 1. Stratospheric species abundances. Abundances given 
here are relative to total number density by volume. Species 
are listed in order of increasing mass. 


Species 

Log (Abundance) 

Species 

log (Abundance) 

h 2 o 

-5.3 

16o1 70160 

-8.7 

HOO 

-9.1 

1 6()1 6q 1 7 o 

-8.3 

HCN 

-9.8 

160180160 

-7.8 

CO 

-7.0 

160160180 

-7.5 

h 2 co 

-10.0 

CtO 

-9.4 

NO 

-8.0 

HOC! 

-9.8 

°2 

-0.7 

37 C 40 

-9.9 

H0 2 

-9 7 

H0 37 Ct 

-10.3 

16 0 18 0 

-3.1 

OCS 

-10.7 

^2^2 

-9.5 

HN0 3 

8.7 

N20 

-7.0 

S02 

-11.0 

NOj 

-8.0 

7 ®Br0 

-11.1 

0, 

-5.2 

81 BrO 

-11.1 

3 






Figure 2. Weighting functions for limb and spectral line shape 
microwave measurements. These functions describe the vertical 
resolution of the technique; the measured signal is proportional 
to the integral, through the observed atmospheric path, of the 
product of this function and the mixing ratio of the species 
being measured. The curves shown here have been normalized 
to unity area. The limb curve includes smearing by the BMLS 
antenna when pointed at 15 km tangent height from a float 
height of AO km; for higher tangent heights this curve is nar- 
rower. The line shape curve is calculated for infinitesimally 
spaced spectral channels with infinitesimal resolution. 
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INSTRUMENT DESCRIPTION 

(a) Summary 

The BMLS consists of (1) an antenna and optics system; (2) a heterodyne 
radiometer with three intermediate-frequency output bands; (3) electronic filter 
banks for spectral analysis; (4) command and data handling systems; (5) power 
supply and thermal control systems; and (6) a ground support system. An overall 
block diagram is shown in Figure 3. The BMLS mass (including batteries for a 24 
hour flight) is approximately 350 kg and it consumes approximately 300 W power. 
Additional mass of the gondola required for a solo BMLS flight is approximately 
300 kg. 

(b) The Antenna and Optics 

The antenna is an offset Cassegrain with a 50 cm diameter main reflector of 
machined aluminum. Vertical scanning of the field of view is performed by a 
flat mirror in front of the main reflectOT. This scanning mirror is pivoted at 
its lower end and discretely driven by a stepping motor-screw drive mechanism in 
0.2° increments. A shaft encoder gives the mirror position. The scan sequence 
is read from a programmable memory. There are 64 scan positions available. 

Figure 4 shows measured 205 GHz antenna patterns at every fourth scan 
position. The measured full width at half power points on the patterns is 
0.30°. For a 40 km balloon float height this width corresponds to a vertical 
resolution at the limb of 3 km for an observation path having 15 km tangent 
height and 1.5 km for a path having 35 km tangent height. 

(c) The Radiometer 

Initial flights with the BMLS were with a radiometer having an ambient 
temperature second harmonic mixer (ZIMMERMAN and KATTAUCH, 1979) and a phase- 



Figure 3. Overall BMLS block diagram 
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Figure 4. Measured BMLS antenna patterns. The circled 
numbers indicate the scan position for that pattern. 


locked 104 GHz klystron local oscillator. The intermediate-frequency (i.f.) 
output of the mixer was split into one band at 4 GHz and another band at 2 GHz. 
GAS-FET amplifiers were used for the first stage amplification. The overall 
single sideband system noise temperature of this radiometer was approximately 
10,000 K. (For microwave radiometry it is customary to measure power in terms of 
the equivalent temperature of a thermal source.) 

To improve sensitivity by a factor of approximately 10, the initial BMLS 
radiometer has been replaced by a radiometer with a liquid nitrogen cooled mixer 
and preamplifier. A block diagram of the cooled radiometer is shown in Figure 
5. The mixer, fabricated at JPL, is a fundamental waveguide mixer using 
Schottky diodes from the University of Virginia. A phase-licked commercial Gunn 
oscillator at 67.7556 GHz is tripled, using a waveguide multiplier developed at 
JPL, to provide the local oscillator at 203.2669 GHz. The local oscillator and 
signal from the BMLS antenna are combined with a folded quasi-optical Fabry- 
Perot whose design principles are discussed by PICKETT and CHI0U (1983). 



Figure 5. Block diagram of the cooled 205 GHz BMLS radiometer 
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The optical layout is shown in Figure 6. This layout also indicates a 
dichroic plate where a bean will be separated for input to a 63 GHz radiometer 
which will be added to the BMLS to test the technique of measuring tangent point 
pressure on the OARS MLS. The layout also indicates the location of a polariz- 
ation grid which can be used to separate the beam for a 270 GHz addition to the 
BMLS which is discussed later. The optical components were fabricated at JPL. 
The dewar was procured commercially; its 5 liter volume gives a measured hold 
time of approximately 80 hours. 

Table 2 gives the noise budget for the radiometer. The rms noise AT for 
a filter bandwidth B and integration time t is obtained from the overall noise 
temperature T^ (given in this table) by the formula 

AT = CT N /(Bt) 1/2 

where C = 1 for total power measurements and C ** 2 for square wave switching. 

(d) The Spectral Analysis System 

The three spectal bands are down-converted in frequency by heterodyning. 
Three identical filter banks then simultaneously perform spectral analysis of 
these bands. 

The heterodyning arrangement used for BMLS flights performed to date is 
shown in Figure 7. For these flights each filter bank contained 35 filters. 



Figure 6. Optics layout of the cooled 205 GHz BMLS radiometer. 






Table 2. Cooled 205 GHz BMLS radiometer noise budget. All 
values are measured and are single sideband. 
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Figure 7 


Heterodyning arrangement used on BMLS flights to date 
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All filters operate simultaneously so that no spectral scanning is required. 

The individual filters were six-pole and constructed of discrete elenents on 
printed circuit boards. Their positions and widths were chosen to fully cover 
and resolve the spectral lines. Three widths (nominally 2, 8 and 32 MHz) are 
used in the filters. Figure 8 shows examples of measured filter response. The 
output from each filter is detected with a square-law diode. The detected 
signal is then passed through synchronous filters which give two outputs propor- 
tional to the power received on the two halves of the switching cycle. 
Differencing of the signals from the two halves of the switching cycle is done 
during data analysis by the BMLS ground support system computer. 

i 

For future flights the BMLS filter banks are being replaced by copies of 
filter banks developed for the UAHS MLS. These filter banks have a center 
frequency of 400 MHz and contain 13 channels. They are composed of both lumped 
element (L-C) and surface acoustic wave (SAW) filters and have better spectral 
characteristics, greater thermal stability, better integrators, and less 
digitization noise than the current BMLS filter banks. Table 3 gives the filter 
characteristics. Figure 9 shows their coverage of the 204 GHz C&0 line. The 
horizontal width of each bar indicates the width for that filter. The vertical 
bar gives the rms instrument noise with a one minute integration for measure- 
ments in total power mode and, equivalently, for a four minute integration for 
measurements with a 50 percent chopping as was the case in previous BMLS 
flights. The amount of chopping required will be between these limits and must 
be determined by tests on the full system after it is completely assembled. 

The power at each filter output is measured by a tunnel diode detector, 
chosen because of its low noise and low thermal drift. An input matching net- 
work provides a nominal -31 dBm power level which is well within the tunnel 
diode square law range, and which gives an output voltage of approximately 1 mV. 




INPUT FMQUINCY/MMi 


Figure 8. Examples of measured BMLS filter response. 
Top: a narrow bandwidth filter; Bottom: a broad 
bandwidth filter. The noise at the dashed line 
level is due to the measurement scheme. 



BRIGHTNESS TEMPERATURE (Kelvin) 
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Table 

3. UARS 

MLS filter 

bank characteristics. 

Channel 

Center 

frequency 

'o 

3 dB Bandwidth Percentage 

Bandwidth 
B 

Filter 

Type 


(MHz) 

(MHz) 

(B/f 0 x 100%) 


1 

209 

128 

61.2 

L-C 

2 

305 

64 

21.0 

L-C 

3 

353 

32 

9.1 

L-C 

4 

377 

16 

4.2 

L-C 

5 

389 

8 

2.1 

SAW 

6 

395 

4 

1.0 

SAH 

T 

398 

2 

0.5 

SAW 

8 

400 

2 

0.5 

SAW 

9 

402 

2 

0.5 

SAW 

10 

405 

4 

1.0 

SAW 

11 

411 

8 

1.9 

SAW 

12 

423 

16 

3.8 

L-C 

13 

447 

32 

7.2 

L-C 

14 

495 

64 

12.9 

L-C 

15 

591 

128 

21.7 

L-C 



Figure 9. Coverage of the 204.352 GHz CIO spectral 
line by the UARS MIS filter bank. 
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A low noise/lov drift Precision Monolithics OP-27E operational amplifier 
provides first stage amplification of the signal from the detector. This signal 
is further amplified by a national LM108A operational amplifier to a level of 
approximately 3 volts. The voltage is then applied to a voltage to frequency 
converter whose output is counted for 1.8 s, and the counter contents read by 
the data handling system. 

(e) The Command and Data Handling System 

A block diagram of the BMLS command system is shown in Figure 10 and a 
block diagram of the data handling system is shown in Figure 11. These systems 
interface the BMLS with the telemetry system of the balloon launch facility. 

(f) Power Supply and Thermal Control Systems 

Electrical power for the BMLS is provided by lithium cell batteries. 
Typically, a flight will include approximately 15 packs of 11 cells per pack. 
This provides approximately 450 ampere hours at 28 volts. The battery voltage 
is converted and regulated to various dc voltages required by BMLS subsystems. 

During flight the instrument is thermally insulated by completely covering 
the gondola and certain subassemblies with thermal insulation. The optical path 
in front of the antenna is covered with 4 cm thick styrofoam which has a loss at 
205 GHz measured to be less than 1%. Polyethylene foam, which is more resistant 
to ultraviolet, is used elsewhere for thermal insulation. Thermal control is by 
regulating the exchange of heat between the instrument and a radiator plate. 

Heat exchange is by circulating ethylene glycol between the heat producing 
elements and the radiator. Regulation is by feedback control of the pumps which 
circulate the coolant fluid. There are three separate thermal control circuits 
which independently regulate (1) the local oscillator assembly, (2) the inter- 
mediate frequency amplifier assembly, and (3) the filter banks and data system 
assembly. 

Figure 12 shows selected BMLS temperatures throughout the flight of 11 May 
1981 from the National Scientific Balloon Launch Facility in Palestine, Texas 
(32°N). Also shown is the balloon height and outside air temperatures . Below 


DATA 
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Figure 10. The BMLS command system. 
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Figure 11. The BMLS data handling system. 


40 km the outside air temperatures were measured by a radiosonde launched from 
Longview, Texas (~100 miles north of Palestine), at approximately 6 a.m. on 12 
May 1981. The outside air temperature at 40 km is from climatology for 30°N May 
(CIRA, 1972). 

(g) The Ground Support System 

The BMLS ground support system processes the data stream received from the 
balloon launch facility telanetry. It produces a magnetic tape record of the 
unprocessed data, performs averaging and initial calibration of radiometric 
data, and provides outputs needed for in-flight decisions. The system is based 
on a Hewlett-Packard HP-1000 computer and is housed in an air-conditioned air- 
ride trailer which is taken to the launch facility for flights. 

During flight several programs are simultaneously run on the ground support 
systan computer. The first program receives the data frame from the telemetry, 
produces a one-line screen output indicating the data were received and the 
antenna scan position and other information on the data frame, and integrates 
the data by either a preset integration time or whenever a configuration change 
(e.g., change in scan position) occurs in the instrument. The second program 
reads a data buffer produced by the first program, performs calibration of the 
data and stores it in a disc file, and creates a log output describing the 
received data. A third program produces a continuous screen display, with 
optional hard copies, of all engineering data. A fourth program accesses the 
disc files of calibrated data and performs various operations on these data such 
as averaging and plotting. 

DATA ANALYSIS 

The first step in data. analysis is to radiometrically calibrate the outputs 
from the BMLS filter channels. This is accomplished by using the linearity of 
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Figure 12. Measured temperatures during the BMLS flight 
of 11 May 1981. The balloon height and outside air 
temperatures are also shown. 


the detected filter output relative to the input power and two calibration 
references. One calibration reference iB the ambient temperature target at the 
radiometer input. The second calibration is the "sky view" at approximately 50 
degrees elevation. Tests on the current BMLS indicate this calibration is 
accurate to approximately 5Z. Flans are underway to add the UARS MLS calibra- 
tion target to the BMLS which should provide IX or better calibration accuracy. 

A second step in data analysis is to calibrate the BMLS pointing. The 
pointing reference is the sharp increase in water vapor emission as the BMLS 
f ield-of-view is scanned downwards through the tropopause. The measured 
increase in emission is compared with that calculated for the water vapour 
distribution from radiosonde measurements made at nearly the same time and place 
as the BMLS measuresents. Figure 13 shows results from the 22 September 1983 
flight in which there was an 0.3° offset between nominal and inferred pointing. 
The theoretical expressions for calculating the water vapor emission are given 
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BRIGHTNESS TEMPERATURE (Kelvin) 

Figure 13. Measured and calculated emission for the BMLS 
field-of-view scanned through the tropopause. On this 
flight the pointing was approximately 0.3 degrees (one 
beamwidth) from the nominal a priori value. Approximately 
15% loss in the optics accounts for the difference be- 
tween measurement and calculation below -6.3 degrees 
elevation. The two calculated curves are for 0 and 15% 
relative humidity between 10 and 18 Von. The radiosonde 
measurement o? humidity stopped at 10 km where 15% 
relative humidity was measured. 


in WATERS (1976); the Van Vleck - Weisskopf line shape was used and the Gaut - 
Reifenstein empirical correction was included. Calculated emission was also 
compared with ground-based measurements of 205 GHz emission versus elevation 
angle at a time when the total column water vapor was measured by an infrared 
solar absorption instrument calibrated against radiosonde measurements. The 
comparison between measurements and calculation is shown in Figure 14. Agree- 
ment is within 10% which translates to a BMLS pointing error of approximately 
0.1 degree or less. 

The calibrated BMLS data are then used in a retrieval process to obtain 
vertical profiles of the species being measured. The retrieval technique now 
being used for the BMLS is one of constrained least squares. Say {p . } are a 
set of parameters to be estimated, e.g., points on a vertical profiliL We make 
an initial estimate of these parameters and from them calculate the brightness 
temperature T f or BMLS channel i 
i 

t b = /// W({p • }, ...) ds dv dfi (1) 

i J 

where the integrations are over the observation path s, solid angle 0, and 
channel frequency interval v. The antenna gain pattern, filter spectral 
characteristics, and effects of radiative transfer and other atmospheric para- 
meters are included in W({p.}, ...). From the preceding equation we 
numerically calculate the partial derivation 
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Figure 14. Comparison of ground-based measurements and cal- 
culations of atmospheric emission at 203 GBz. At the time 
and place of the 205 GHz measurements a calibrated infrared 
solar absorption instrument measured 1.1 ±0.1 gm/cm2 
vertically integrated water vapor. Calculations are given 
for 1.06 and 1.15 gm/cm2 integrated water vapor. 


C. , = 3T 


i/3 Pj 


B 

Expanding (1) in a Taylor series of Ap. about p. and keeping linear terms 
gives ■* 


( 2 ) 


V 


= % + C-A £ 


(3) 


where a single underline indicates a vector and a double underline indicates a 
matrix. We now require that the difference between calculated brightness 
and measured brightness T^ be approximately equal to the measurement noise 


e 2 * (i/N) £<T _ Tj^) 2 (4) 

where £ is the rms measurement noise and N is the number of measurement 
channels. If we require better agreement than this, then we will likely 
retrieve a profile fit to noise which ie undesirable. Equation (4) imposes a 
constraint on the solution for Ap.. As written above it assumes equal noise 
in all channels; modification for J unequal channel noise is straightforward. 


In general there is no unique solution for Ap. due to the smearing 
described by the integrals in (1) and due to uncertainties introduced by 
measurement noise. We must establish a criterion for choosing one of several 
possible solutions which agrees with the measurements to within the noise. One 
reasonable criterion is choosing the "smoothest" profile consistent with (4). 

The smoothest profile as measured by the second differences of the points on the 
profile is obtained by minimizing the scalar quantity 



where superscript T indicates transpose and the square matrix K is given by 
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K = 


0 0 0 0 
-1 2-1 0 

0 - 12-1 
0 0-1 2 


2-10 
-1 2 -1 

0 0 0 


( 6 ) 


whose numberof rows and columns equals the number of elements in the parameter 
vector . 


Minimizing (5) subject to the constraint (6) gives the solution 

A_£> = [C T C + eK T K) _1 . C T • (Tg - T b ) (7) 

which can be integrated if necessary to account for nonlinearities. The para- 
meter C in (7) determines the relative weight given to smoothing versus fitting 
to measurement. If 4 is chosen to be zero no weight is given to the smoothing 
requirement and the standard least squares solution is obtained; as £ is made 
larger less weight is given to the measurement. In practice £ is set to the 
value which makes the difference between calculated and measured brightness 
equal the noise. Selection criteria other than maximum smoothness as defined by 
second differences can be implemented simply by changing the form of K from 
that given in (6). TWOMEY (1977) gives some examples and further discussion of 
this point. RODGERS (1976) gives a more sophisticated and comprehensive dis- 
cussion of retrieval techniques. One feature of microwave limb sounding which 
is of significance to the retrieval process is that only those elements of the 
data vector T^ which correspond to optically thin signals at all points in 
the spectral channel being measured are included in the retrievals. This makes 
the procedure much more linear than would otherwise be the case and, for 
example, means that the partial derivatives (2) need only be calculated once. 

This profile retrieval technique was used to determine the accuracy of 
obtaining CIO and 0. profiles from the BMLS measurements. In this process the 
BMLS measurement is simulated by calculating the emission for a specific pro- 
file. This simulated measurement is then used in the retrieval algorithm to 
obtain a profile which is compared with the original profile. Error bars for 
the retrieved profile points are obtained from the retrieval matrices using the 
expected measurement noise. For simulated retrievals shown below, the species 
mixing ratio profile was represented as a continuous piece-wise linear function. 
Values at the ends of the individual linear segments were inferred by the 
retrieval process. The smoothing parameters c was set to zero in order to 
obtain independence of the different points. 

Results of the C£0 retrieval simulations are shown in Figure 15. The two 
cases shown are, respectively, for mid-day and midnight C£0 profiles as 
predicted by current theory. The simulated measurement was for the 15 UARS MLS 
filters, eight elevation angles, and a typical balloon altitude of 40 km. The 
error bars are for an integration time per elevation angle corresponding to 1 
minute for non-switched (or 4 minute for switched) operation of the radiometer. 
Thus, complete profiles with the accuracies indicated in Figure 15 should be 
possible with time resolution of between" 8 and 32 minutes. 

Results of the 0, retrieval simulations are shown in Figure 16. As for 
the CIO retrievals, the simulated measursnent was for the 15 UARS MLS filter 
channels and eight elevation angles from 40 km height. The error bars are for 
the same integration time as for CLO and for the larger system noise taaperature 
of the 0^ band (see Table 2). These results indicate that very accurate 0, 
profiles can be retrieved from the BMLS measurement. They also danonstrate the 
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Cio VOLUME MIXING RATIO fPPBV) 


Figure 15. Results of simulated CIO retrievals. The solid 
curves are the "true" profiles and the dashed curves the 
retrieved. The horizontal extent of the crosses on the 
retrieved profiles indicate their rms uncertainty for the 
measurement sensitivity of the instrumentation discussed 
in the text with one minute total power integration time 
at each of eight elevation angles. Profile points above 
40 km have degraded accuracy because they are obtained 
from the information in the spectral line shape instead of 
information in the limb scan. 



0 123456789 

0 3 VOLUME MIXING RATIO (PPMV) 


Figure 16. Results of simulated 0. retrievals. The two 
profiles are, respectively, for the KRDEGER-MINZNER (1976) 
standard 0^ profile and this standard profile minus three 
standard deviations as also given by Krueger-Minzner. The 
solid curves are the "true" profiles and the dashed curves 
the retrieved. The horizontal extent of the crosses on the 
retrieved profiles indicate their rms uncertainty for the 
measurenent sensitivity as discussed in the text. 
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capability of the technique to measure vertical structure with 3 km or better 
resolution. 

MEASUREMENTS OF HO,, NO,, UNO, , HCN , N,0 AND 16 0 18 0 16 0 BY ADDING 
A 270 GHz RADIOMETER TO THE BMLS Z 


Understanding of the stratosphere is now so sufficiently advanced that 
simultaneous measurements of several species are extremely important, if not 
essential, for further progress. Tn order to make such measurements practical, 
both logistically and economically, it is desirable to obtain several species 
from single instruments. In the case of the BMLS. many additional species could 
be measured as suggested in Figure 1. However, because of resource limitations, 
compromises must be made. One very attractive compromise is to add to the BMLS 
a single 270 GHz radiometer which simultaneously measures six additional species 
of HO, , NO,, HN0 3 , HCN, N,0 and 16 0 I8 0 16 0. This radiometer also 
simultaneously measures 0^ . 


Some specific justifications for these additional measurements are 
discussed in the following subsections. 


(a) Test of 0^ Photochemical Balance 

The largest percentage depletion of 0. by CFC's is predicted to occur 
between heights of 30 and 45 km. Current theory (WUEBBLES, 1983) predicts a 40% 
0 3 depletion at 40 km if CFC releases continue at their present rate. This 
region should thus provide the best early warning of these effects. The 
chemistry there is expected to be less complex than at lower altitudes, and 
reactions are sufficiently fast that transport processes should not upset photo- 
chemical balances. One expects, therefore, the theory for this region to be 
more reliable than for lower altitudes. However, there is a serious discrepancy 
between measured and calculated 0, abundances above about 35 km. Measurements 
are approximately 50% higher than 'calculated, as shown in Figure 17. The 
discrepancy is not likely due to measurement errors, as several measurement 
techniques have yielded similar results. It is also not likely due to 0 3 
variability, which is measured to be less than approximately 20% at these 
heights. This discrepancy must be resolved in order to have confidence in the 
theory. 


The 0, photochemical balance expected to occur above about 35 km can, 
when second order effects are neglected, be written (e.g., FR0IDEVAUX, 1983; see 
also ANDERSON et al., 1980): 

J 0 ,[0,] = (J 03 [0 3 ]/k 1 [0,][M])(k,[0 3 ] + k 3 [H0 2 ] + k^NO^ + k^ [CIO] ) 

where [A] is abundance of species A. The left side is the 0, production rate, 
limited by the rate Jq, of breaking the 0, bond. The right side is the 0 3 
destruction rate. The first factor on the right is the predicted abundance of 
atomic 0, in equilibrium between 0 3 photolysis at rate Jq 3 and recombination 
with 0, (in the presence of another molecule M) at rate k, . The rate 
constants k,,k 3 ,k^,k 3 refer, respectively, to the 0 3 destruction rate- 
limiting reactions: 0 + 0 3 , 0 + HO,, 0 + NO,, and 0 + CIO. 

All abundances in the above equation are measured by the BMLS with the 270 
GHz addition ([0, ] and [M] are obtained from the 63 GHz pressure measurements 
which will be made when the UARS 63 GHz radiometer is added to the BMLS). This 
provides the possibility of the BMLS performing a needed test of the 0 3 photo- 
chemical balance. Such a test, which has not been done to date, should give 
valuable insight concerning the cause for the discrepancy between measured and 
calculated 0, in the upper stratosphere. 
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Figure 17. Measured and theoretical 0, profiles. The mid- 
latitude Caltech equinox model is shown (see FROIDEVAUX (1983)1. 
the KO and SZE profile is also theoretical, for equinox and 
30°N. The U.S. Standard Atmosphere Profile is that of 
KRUEGER and M1NZNER (1976) based on measurements . The IORI 
profile was measured during an international ozone rocket 
intercomparison campaign. The BMLS measurements made on 22 
September 1982 from Palestine, Texas (32 < *N) are also shown. 

The full widths of the horizontal bars on all measurements 
are two standard deviations. 


The accuracy of the BMLS 0^ photochemical balance test can be estimated 
by considering the relative importance of the various terms in the balance 
equation, and the accuracy with which the BMLS can measure the respective 
abundances. At 35-40 km altitudes the relative importance of the terms in the 
second factor on the right side of the balance equation are (FROIDEVAUX, 1983) 
k 2 l0 3 J k 3 lH0 2 3 ^ 10% ; k^NO-] ^60%; k.tClOJ %20%. The 

appropriate abundances which go into the equation are daytime averages, so a 
reasonable time to allow for the BMLS measurements is one hour. With one hour 
integration the BMLS accuracies for daytime measurements will be: ■vll for 0^ , 
' v '5% for H0 2 , ' v '25% for NO 2 , and ^5% for C 0. Thus the uncertainty in the 
test due to BMLS day time measurement uncertainty will be %20%. Additional 
uncertainty will be introduced by uncertainties in Jq 2 , Jqj and k Q . The 
uncertainty in J_ and Jqj should be less than 10% (Froidevaux, private 
communication; the BMLS measurement gives total 0^ overburden which reduces 
uncertainty in Jqj). The effects of uncertainties in k Q should also be less 
than 10% (these have been measured and the atmospheric temperature profile at 
the time of BMLS flights is also measured which allows the k Q temperature- 
dependence to be accounted for). Thus, an overall accuracy of ^30% is expected 
for the 0^ photochemical balance test from daytime BMLS measurements. This 
accuracy is useful since the discrepancy between measured and calculated 0, in 
the upper stratosphere is ‘'-50%. (It should be possible to further improve tne 
accuracy of the test by using the more accurate BMLS nighttime NO 2 measure- 
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ments to help infer average daytime NO- values; M. Allen, private communica- 
tion.) z 

Since the uncertainty in the BMLS 0. photochemical balance test is 
largely dominated by the uncertainty in the N0 2 measurenent, simultaneous 
N0 2 measurements by other techniques would be very desirable. A technique 
which can provide these measurements both day and night, as well as measurements 
of the important related species NO, is balloon-borne infrared pressure- 
modulated radiometry (PMR), (DRUMM3ND and JARNOT, 1978; R0SC0E et al., 1981). 

A balloon PMR can be flown together with the BMLS. 

It would also be very desirable to have simultaneous measurements of [0] to 
independently test the validity of the first factor on the right side of the 
preceding photochemical balance equation. One measurement has been performed 
which indicates this expression is valid (see WMO, 1908; Figure 1-54) more 
evidence for its validity is needed, since as noted by FROIDEVAUX (1983), the 
measured 0- abundances in that test were substantially lower than most 0^ 
measurements (and in atypical good agreement with theory). 

(b) H0 2 

Microwave spectroscopy is the only demonstrated method for remotely 
measuring H0 2 (DE ZAFRA et al., 1984) which, with OH, is intimately involved 
in all of what are now thought to be the major catalytic cycles of stratospheric 
0j destruction. 

HO x strongly influences the NO x catalytic destruction of by 
reactions such as: 

H0 2 + NO -* N0 2 + OH 


and 


N0 2 + OH + M -*■ HN0 3 + M 

where the HNO^ diffuses downward and is removed by rain. 

HO increases the CM> x catalytic destruction of CU by converting 
inactive? chlorine in HC2. to active atomic chlorine by fne reaction 

OH + HCZ H 2 o + C l 

H0 x itself catalytically destroys 0 3 by reactions such as: 

OH + 0 3 ■* H0 2 + 0 2 


and 


H02+ 0 -*■ OH + 0 2 

and by a similar chain where H0 2 reacts with 0 3 . 

Thus, knowledge of the distribution and variation in stratospheric H0 x is 
key to general understanding of stratospheric 0 3 destruction. Below about 
40 km H0 2 is expected to be the most abundant stratospheric species of H0 X 
but only a few stratospheric H0 2 measurements have been made to date. 

Chemical models of the stratosphere are especially sensitive to the 
0H/H0 2 ratio and, for a specific test of stratospheric chemistry theory, 
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simultaneous OH measurements are needed with HO-. Such measurements could be 
made from a single balloon carrying the BMLS ana submillimeter or far-infrared 
instruments. It would also be desirable to simultaneously measure H^O and 
other chemically coupled species; this could also be done from the same balloon. 

(c) N0 2> HNO 3 , N 2 

These species have key roles in the N 0 x catalytic destruction of O, 
which is predicted to account for approximately 70% of total 0^ destruction. 
Understanding of this catalytic cycle is crucial. 

N 0 2 is a key radical m the cycle, for which the rate limiting reaction 
is 


N0 2 + 0 NO + 0 2 
which, coupled with 

NO + 0 3 -> N0 2 + 0 2 

destroys 0,. However, there are only a limited number of N0» measurements 
and most of these were taken during sunset or sunrise which does not observe the 
diurnal variation in N0 2 . The BMLS N0 2 measurements throughout a complete 
diurnal cycle could provide this essential information. These will be 
especially valuable when made simultaneously with NO measurements. 

HNO 3 removes active N 0 X from the lower stratosphere by downward 
diffusion and rainout. Above 25 km, theoretical models of the stratosphere have 
predicted significantly more HN0- than measured, by as much as a factor of 3 
at 30-35 km (e.g., WMO, 1981). very recent models, e.g., FROIDEVAUX (1983) may 
have reduced this discrepancy, but more reliable measurements of connected 
species are needed. A specific test of model predictions would be obtained by 
simultaneous measurements of HNO^ , N0 2 and OH, all of which could be per- 
formed throughout a complete diurnal cycle from a single balloon carrying the 
BMLS and far infrared or submillimeter instruments for OH. 

N 2 o, upwardly transported from the troposphere and oxidized by O(^D), 
provides the major source of stratospheric NO . The measurement of the 
decrease in stratospheric N 2 0 mixing ratio gives information on this source 
and transport processes. 

(d) HCN 

Stratospheric HCN has been measured from aircraft-based infrared solar 
absorption by COFFEY et al. (1981), and from balloon-based submillimeter thermal 
emission by (ARLI et al. (1982). Only very recently (CICERONE and ZELINER, 

1983) has a stratospheric model been developed which includes HCN chemistry. 
Predictions of HCN from this model agree with the Coffey et al., measurements in 
the lower stratosphere, but are a factor of 2 below the Carli et al . , measure- 
ments at 40 km. The BMLS measurements of HCN will extend to 40 km or higher, 
and will provide additional data needed to test the model. 

(e) 16 0 18 0 16 0 

Balloon-borne mass spectroscopic measurements of MAUERSBERGEK (1981) show a 
pronounced altitude-dependent enchancement of heavy ozone (O 3 , containing one 
atom). The measured enhancement is largest at 32 km altitude where it is 
a factor of 1.4 above the expected 0/ 0 ratio. Although enhanced photo- 

dissociation of heavy 0, in the upper stratosphere is predicted (CICERONE and 
MCCRUMB, 1980; BLAKE et al., 1983), reactions between 18 0 and 3 Z 0 2 are 
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sufficiently fast that no such enhancement in the abundance of or heavy 
ozone is expected (FROIDEVAUX, 1983; KAYE and STROBEL, 1983). Thus, the 
Mauersberger heavy ozone measurements are inconsistent with current theory and 
it is important that they be checked. 

(f) o 3 

The BMLS 270 GHz radiometer also measures a strong ozone line at 276.923 
GHz. This line provides an independent check on the current BMLS ozone measure- 
ment using the 206.132 GHz line. Such a check is important to help verify that 
ozone measurements can be performed by the BMLS and UARS MLS to the expected 
accuracy of approximately 1 per cent. 

EXPECTED 270 GHz SIGNALS AND MEASUREMENT ACCURACY 

The signals expected for the 270 GHz measurements can be determined by 
using the expected stratospheric abundances of the respective species in an 
atmospheric radiative transfer calculation (see, e.g., WATERS, 1976). The 
calculations done here break the atmosphere into layer approximately 0.5 km 
thick and integrate the radiative transfer equation through the limb path, 
accounting for spherical atmospheric geometry. The full Voigt spectral line 
shape and the 1976 U.S. Standard Atmosphere model of temperature and pressure 
are used. Strengths of the respective lines, and their temperature dependence, 
were taken from the JPL catalogue (POYNTER and PICKETT, 1981), with latest 
revisions as of March 1984, and are believed accurate to 1% or better (all 
appropriate dipole moments are known to 0.1% and uncertainties in the transition 
matrix elements do not exceed 1%). Line width uncertainties introduce less than 
5% uncertainty in retrieved profiles, based on simulated retrievals; this is due 
to height information being obtained principally from limb scan and not line 
shape. Frequencies of all the lines have been measured to the required 
accuracy. 

Figures of calculated signal levels given in the following subsections also 
include a cross showing the noise level of the 270 GHz instrumentation 
(described later) for an indicated integration time. The vertical extent of 
this cross is 2 standard deviations for the instrument operated in non-switched 
mode and 1 standard deviation for square-wave switching. 

(a) H0 2 

Expected abundances of H0 2 are shown in Figure 18. The recent ground- 
based measurements of DE ZAFRA et al. (1984) have now demonstrated that H0 2 
can be measured from microwave emission lines near 270 GHz. In fact, this is 
the only demonstrated method by which H0 2 can be remotely measured. 

Figure 19 shows calculated emission from the 265.770 GHz H0 2 line which 
is one of the triplet of H0„ lines near 265.7 GHz measured by DE ZAFRA et al. 
(1984). All three of the HD 2 triplets would be simultaneously measured by the 
BMLS 270 GHz radiometer. A pressure-broadening coefficient of 3.1 MHz/mb at 
T=300 K with T -0- ' 5 temperature dependence was used for the calculation, 
following DE ZAFRA et al . (1984). The curves in Figure 19 indicate the expected 
range of daytime H0 2 signals. The lower profile 2 is consistent with the de 
Zafra et al., ground-based measuranent. 

Since there is a significant mesospheric contribution to the H0 2 signal, 
it is difficult to estimate, from only measuranent signal-to-noise, the accuracy 
with which stratospheric H0 2 profiles can be retrieved from BMLS measurements. 
To determine this accuracy we used the profile retrieval procedure described 
earlier. Results are shown in Figure 20. The simulated measuranent for these 
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Figure 18. Stratospheric HOj abundances. From WMO (1981). 



Figure 19. Calculated limb emission from the 265.767 GHz EO^ 
line. These calculations are for an elevation angle of -1 
degree from 40 km. The tvo curves are for dashed profiles 2 
and 3 of Figure 18, respectively. Also shown is the rms noise 
and the spectral resolution for the instrumentation described 
in the text with 4 minute integration. 
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HC> 2 VOLUME MIXING RATIO (PPBV) 

Figure 20. Results of simulated H0 2 retrievals. The dashed 
lines give the retrieved profile and the solid lines are the 
"true" profile. The horizontal width of the retrieved profile 
points gives the rms uncertainty due' to the instrument noise 
as discussed in the text. 


retrievals used the same 8 elevation angles as for the CIO retrievals described 
earlier. Channels 2 MHz wide within 20 MHz of line center were used. The 
initial estimate for the retrievals was dashed curve 1 of Figure 18 and 
retrievals were performed on the simulated data vector calculated separately for 
dashed profiles 2 and 3. The rms measurement noise was that for the 270 GHz 
instrumentation described with 4 minute integration per elevation angle for non- 
switched radiometer operation (equivalent to 16 minutes per elevation angle for 
switched operation). 

The results shown in Figure 20 indicate that useful HO^ accuracy can be 
obtained to altitudes as low as approximately 25 Von. In fact integration times 
a factor of four shorter than those assumed for the retrieval simulations done 
here should provide useful measurements. This means that useful H0 2 profiles 
between 25 and 40 km could be obtained with time resolution of approximately 1/2 
hours or better. 

(b) N0 2 

Expected stratospheric abundances of NO. are shown in Figure 21. The 
symbols are measurements made near sunset (WHO, 1981) and the solid curves are 
predictions by the Caltech model (L. Froidevaux, private communication) for noon 
and midnight. Figure 22 shows calculated emission by the 265.560 GHz N0 2 line 
for the noon and midnight theoretical profiles of Figure 2.3-5 and for elevation 
angles of -1° and -3° from 40 km. A pressure-broadened linewidth parameter of 
3.0 MHz/mb with T * temperature dependence was used for the calculations. 

This line is one of three NO, lines which could be simultaneously measured by 
the 270 GHz instrumentation described below. The BMLS NO, signals are 
stronger at night because of the nighttime increase in N0 2 . A 15 minute inte- 
gration at midnight will give an expected signal to noise of approximately 10 in 
the 30-40 km altitude region. During midday a one hour integration is required 
for a signal-to-noise of approximately 10 at 30 km; one hour integration gives a 
signal to noise of approximately 2 at 40 km. Data from all three N0 2 lines, 
which have approximately equal strengths, can be combined to improve signal-to- 
noise by a factor of approximately 1.7 over these values. 





NOj VOLUME MIXING RATIO (ppbv) 


Figure 21. Stratospheric NO,, abundances. From WMO (1981! 


no 2 



FREQUENCY FROM 265.560 GHz (MHz) 


Figure 22. Calculated limb emission by the 265.560 GHz line 
at elevation angles of -1 degree and -3 degrees from 40 km 
altitude. These angles correspond to observation paths 
having tangent heights of 39 and 31 km, respectively. Not 
the difference in vertical scale for the noon and midnight 
signals. The instrument noise for 20 MHz spectral resolu- 
tion and integration times of 60 minutes at noon and 15 
minutes at midnight are shown. 
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(c) HOT 

COFFEY et al. (1981) detected stratospheric HOT from aircraft infrared 
absorption and obtained an average HOT mixing ratio above 12 km of 0.17 ppbv. 
CAUL I et al. (1982) detected several submillimeter HCN lines from balloon 
measurements of thermal limb emission and obtained a mixing ratio of 0.13 to 
0.26 ppbv in the 20-40 km altitude range. CICERONE and ZELLNER (1983) calculate 
an HOT profile which has 0.15 ppbv at 15 km decreasing to 0.06 ppbv at 50 km. 

To determine the expected strength of the 265.886 GHz HOT line in stratospheric 
emission we used a constant HOT mixing ratio of 0.1 ppbv between 15 and 50 km. 
This profile is intermediate between the theoretical and measured values in the 
middle stratosphere. A collision linewidth parameter of 3.0 MHz/mb quoted_by DE 
ZAFRA et al. (1984) for HOT was used with a temperature dependence of T 
Results of the calculation are shown in Figure 23. It can be seen that this 
line is very strong, having a brightness temperature of nearly 50 K for -1 
degree elevation angle. Measurement of this line with a signal-to-noise of 100 
can be achieved with a few seconds integration. Expected HOT profile accuracies 
are a few percent. 



Figure 23. Calculated limb emission by the 265.886 GHz HCN 
line for -1 degree elevation angle from 40 km altitude. 

The instrument noise for a 2 second integration is indicated. 

(d) hno 3 

Measured stratospheric abundances of HN0, are shown in Figure 24. We 
used the dashed profile in the mid-range of these abundances to calculate HNO^ 
emission at 269.2 GHz. Results of the calculations are shown in Figure 25. 

HNOj emission near this frequency is from 'vlO separate lines which are merged 
together by pressure broadening. A pressure broadening parameter of 3.0 MHz/mb 
at T=300 K based on the 2.2 Debye dipole moment of HNO, and the plot of line- 
width versus dipole moment by DE ZAFRA et al. (1984) was used. The HNO^ 
signal is very strong and signal-to-noise of 100 for a mid-range profile point 
is achieved with an integration time of a few seconds. An integration time of 4 
minutes will provide an expected signal-to-noise of approximately 10 at 40 km 
altitude. 

(e) N 2 0 

Stratospheric N 2 0 profiles are shown in Figure 26. The lower dashed line 
in this figure was used to calculate the N 2 0 limb emission shown in Figure 27. 




n i i T riiTi 


n — i i iiiiii 


i 

»- 

z 

o 


40r 


30 r 


20 r 


FISCHER |»980) 

O MAY 1979 3)°N 
A FEB 1979 31 C N * 

ARNOLD ET AL 1 1980} 

♦ NOV 1977 45°N ____ 

EVANS ET AL <1978 ' * ’| 

■ JUL AUG 74 76 51 °N 

MARRIES ET AL (1976) 

♦ SEPT 1974 AS'-N K 

FONTANELLA ET AL <1975) 

* JUL 1973 48 N 

MURCRAY ET AL (P»tkxuI communication) 

O OCT 1979 32 i N 

LA2RUS AND GANDRUD <19741 

♦ SPRING 1971 I 

• SPRING 1972 32“N 

• SPRING 1973 


s 


HNO, 


S g C ... 


* 




TTTTT 


VERTICAL COLUMN 


(8 6 i 4 01 ■ 10“ cm * 


J I I 


0.1 


I LUti. 


J i limn i i ii mi 


10 


hno 3 volume mixing RATIO (PPBV) 


Figure 24. Stratospheric HNO^ abundances. From WHO (1981). 



FREQUENCY FROM 269.210 GHz (MHz) 


Figure 25. Calculated HNO^ limb emission near 269.2 GHz 
for several elevation angles from 40 km altitude. This 
spectral region could be covered by the OARS MLS filter 
bank with filter positions and resolution as indicated in 
Figure 9. The instrument noise for an intermediate reso- 
lution filter with 2 seconds integration time is indicated. 





Figure 26. Stratospheric midlatitude profiles. From 
WHO (1981). The points are measured and the dashed and 
solid lines are model predictions. 



FREQUENCY FROM 276.328 GHz (MHz) 


Figure 27. Calculated limb emission by the 276.328 GHz NjO 
line for several elevation angles from 40 km altitude. This 
spectral region could be covered by the UARS MLS filter bank 
with filter positions and resolution as indicated in Figure 9. 
The instrument noise for an intermediate resolution filter 
with 2 seconds integration time is indicated. 
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A pressure-broadening parameter of 2 MHz/mb was used. The NjO line is very 
strong and a signal- to-noise of 100 is achieved with a few minutes integration. 

(f) 16 o 18 o 16 o 

. .Figure 28 shows calculated emission from the 26JS.623 GHz line of 

0 0 b 0. The vertical profile of b 0 18 0 ib 0 used for the calcula- 

tions wa6 the standard KRDEGER-MINZNER (1976) 0, profile scaled at all 
altitudes by the factor 2.04 x 10 for the cosmic .abundance of °0 relative 
to rb O (e.g., LEIGHTON, 1959). This abundance of lb 0 i8 0 Ib 0 is con- 
sistent with the ground-levgl mass spectroscopic measurements of MAUERSBERGER 
(1981) which gave 6.3 x 10 relative abundance for total heavy ozone, since 
the expected amount of heavy ozone in the b 0 16 0 lb 0 isomer is one-third 
(e.g., FROIDEVAUX, 1983). A signal- to-noise of approximately 10 throughout 
the 20-40 km altitude range is achieved with a few minutes integration. Any 
enhancement in stratospheric heavy ozone, such as that measured by Mauersberger , 
will increase the signals from those shown here. 

Also shown in Figure 28 is the NOj emission line at 265.643 GHz. 



Figurjg28g .galculated limb emission by the 265.623 GHz line 
of 0 1 0 l 0 with no enhancement of heavy ozone. The 
instrument noise for 4 minutes integration time with 5 MHz 
spectral resolution is shown. 


(g) o 3 

Figure 29 shows calculated limb emission by the 276.923 GHz 0. line using 
the KRUEGER-MINZNER (1976) standard 0^ profile. Measurement signal-to-noise 
of 100 is achieved with a few seconds integration. 

270 GHz INSTRUMENTATION 

A block diagram of the 270 GHz radiometer which could be added to the BMLS 
is shown in Figure 30. This radiometer is very similar in design to the 
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FREQUENCY FROM 276.923 GHx (MHz) 


Figure 29. Calculated limb emission by the 276.923 GHz 0^ 
line for several elevation angles. This spectral region 
could be covered by the UARS MLS filter bank with filter 
positions and resolution as indicated in Figure 9. The 
instrument noise for an intermediate resolution filter with 
2 seconds integration time is indicated. 
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Figure 30. 270 GHz radiometer block diagram. 
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existing cooled 205 GHz BMLS radiometer described earlier. It uses a Schottky 
diode waveguide mixer, and local oscillator obtained by tripling a phase-locked 
Gunn oscillator at 89.33 GHz. Gunn oscillators at 90 GHz are commercially 
available with at least 30 mW output power. The tripler has an efficiency of at 
least 2 percent which will provide at least 0.6 mW at the local oscillator 
frequency of 267.99 GHz. This local oscillator frequency was chosen so that 
there are no known significant spectral interferences in either sideband with 
the HO , N0„ , HCN, HNO- , N-0, and 0- signals being measured. It also 
placed a 259.532 GHz HOCi. line at 8.46 GHz intermediate frequency (i.f.), 
which is within the passband of the 9 GHz i.f. amplifier for N^O and 0-; 
this line is also a candidate for measurement by the BMLS. The Fabry-Perot 
diplexer which combines local oscillator and signal has less than 0.3 dB loss in 

the signal path and less than 1.5 dB loss in the local oscillator path. 

A layout of the 270 GHz radiometer optical components and dewar is shown in 
Figure 31. The 270 GHz signal is split from the BMLS antenna optical path by a 
polarization grid. The location of this grid in the optical path of the current 

BMLS is shown in Figure 6. All BMLS measurements at 63, 205 and 270 GHz can be 

performed simultaneously with this arrangement. The Fabry-Perot which combines 
signal and local oscillator has 3.5 GHz free spectral range and is transmissive 
to all spectral lines being measured. Figure 32 shows its characteristics. 

The signal from the mixer is split into three i.f. bands. The first band 
from 1.0 to 1. 5 .Qlx contains HNO- . The second, from 2.0 to 2.5 GHz contains 
HO , NO-, lb 0 18 0 lb 0 and HCN . The third, from 8.7 to 9.2 GHz, contains 
NjO and 0^ . Signals in these bands are amplified and then converted to 
lower frequency bands for spectral analysis. Refer to Figure 33. 



Figure 31. Optical layout for 270 GHz radiometer. The 
optics producing the local oscillator beam are identi- 
cal to those shown in Figure 6 for the 205 GHz radio- 
meter. 
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Figure 32. Characteristics of Fabry-Perot for 
coupling 270 GHz local oscillator and signal. 



Figure. 33. Spectral band containing HOj , NO 2 , HCN, 
and signals. This spectrum is calculated 

for emission from 40 km at -1 degree elevation angle. 
The estimated instrument noise for 15 minutes inte- 
gration time with 2 MHz resolution is indicated. 
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Spectral analysis of HNO_ , N.O and 0. signals can be performed by 
copies of the UARS filter banks , described earlier. The HO^ and NO, 
signals, consisting of several lines of narrower width than those 01 the other 
species, require a different spectral analysis system. A spectrometer with 
500 MHz bandwjidtli Mid 1-2 MHz resolution would simultaneously measure the HO^ , 
NOj , HCN and signals. Such a spectrometer could be readily 

implemented. 

An estimated noise budget for the 270 GHz radiometer is given in Table A. 


Table A. Estimated noise budget of the 270 GHz radiometer. 
All values are single sideband. 



0.8-1.$ GHz 
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8.0-9. 2 GHz 
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band 
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(HN0 3 ) 

(N 2 0, 0 3 ) 

loss before mixer (dB) 

0.5 

0.5 

0.5 

mixer conversion loss (dB) 

9.0 

9.0 

11.0 

mixer noise temperature (K) 

800 

800 

800 

l.f. triplexer loss (dB) 

1.0 

2.0 

3.0 

l.f. noise temperature (K) 

30 

30 

200 


overall noise temperature, T ^ (K) 1300 1400 6500 
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ABSTRACT 

Recent progress in the technology of stratospheric ion mass spectrometry is 
described. 

The state of the art as well as some typical measuring problems and needs 
for future developments are briefly reviewed. 

INTRODUCTION 

Penetration into the atmosphere of high energy particles and short wave 
radiation such as solar UV light, X rays and cosmic rays, results in formation 
of charged particles through direct ionization processes. At sufficiently high 
gas densities, ions and electrons thus formed make multiple collisions with 
neutral particles, and electron attachments as well as ion-molecule reactions 
occur, modifying the identity of the charged particles. In the stratosphere the 
main ionization source is cosmic radiation and the lifetime of ions is of the 
order of 100 to 3000 s between 15 and 40 km. The subsequent rapid ion-molecule 
reactions involve different trace gases so that the terminal ions have a totally 
different nature than the primary ones. 

Until 1977, the experimental data on stratospheric ions were limited to 
total ion density and mobility measurements (PALTRIDGE, 1965; BRAGIN et al., 
1966; ROSE et al., 1972; WIDDEL et al., 1977; BRAGIN, 1967; MORITA et al., 1971; 
MITCHELL et al . , 1977). 

Information on positive and negative ion composition below 50 km was based 
on modelling (FEHSENFELD and FERGUSON, 1969; MOHNEN, 1971; FERGUSON, 1974; REID, 
1979), and extending the framework formerly developed for D-region ion 
chemistry. Therefore it remained merely speculative and incomplete. 

Nevertheless, such information is important for several reasons: 

- Knowledge of the nature of stratospheric ions is essential for our under- 
standing of physical processes in the field of atmospheric electricity. The 
conductivity of the atmosphere, for instance, depends upon ion mobility and 
thus on ion composition. 

- Also, stratospheric ion mass spectrometry may become an important tool for 
detection of trace gases with very low concentrations (ARNOLD et al., 1980; 
ARIJS et al., 1983a, b; VIGGIANO and ARNOLD, 1983). 

- Furthermore, stratospheric ions may play an important role in gas to particle 
conversion and nucleation of stratospheric aersols (MOHNEN and KIANG, 1976; 
ARNOLD, 1980, 1982), which in turn could influence the Earth's radiation 
budget and climate. 

- Finally, in situ mass spectrometr ic investigations of stratospheric ion 
composition can provide experimental thermochemical and kinetic data on 
ion-molecule reactions taking place under conditions difficult to simulate 
in the laboratory. 

Determination of the nature and abundances of stratospheric ions is there- 
fore a desirable objective and balloon-borne ion mass spectrometers have been 
and are developed at several places (ARIJS et al., 1978; ARNOLD et al., 1978; 
OLSON et al., 1978; CUNNINGHAM and HOFMAN, 1982; BALLENTHIN et al . , 1983). 
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The purpose of this paper is to describe the experimental and technological 
aspects of such ion mass spectrometers as far as they were available in litera- 
ture. The scientific results obtained with these instruments have been reviewed 
elsewhere (ARNOLD, 1980; ARIJS, 1983). 

EXPERIMENTAL PROBLEMS AND PAYLOAD DESCRIPTION 

In situ stratospheric ion identification requires the development of an 
instrument, wherein ions are sampled through a small orifice, focussed into an 
appropriate mass spectrometer and detected after filtering according to their 
mass to charge ratio. This instrument has an associated electronic package 
taking care of driving the mass filter and focussing device, delivering high 
voltage to the detector, recording signals and interfacing between the 
instrument and the ground-borne experimenter. 

For ion composition measurements in the D-region, rocket-borne mass 
spectrometers have been developed by several groups (NARCISI and BAILEY, 1965; 
GOLDBERG and BLUMLE, 1970; KRANKOWSKY et al., 1972; ZBINDEN et al., 1975). 

Use of similar instruments for stratospheric purposes (ARNOLD et al . , 1977) 
however, poses some severe problems among which short measuring times and ion 
breakup due to induced shock waves are the most important. For stratospheric 
measurements it i6 better to make use of more convenient platforms, such as high 
altitude balloons which are able to lift heavy and large instruments up to 45 km 
altitude. Present techniques even allow continuous measurements in an altitude 
range from 45 down to about 15 km by using valve controlled balloons. 

Although the philosophy in designing a balloon-borne mass spectrometer for 
in situ sampling of stratospheric ions is similar to the one followed in 
constructing a rocket-borne payload, some important differences evolve from the 
nature of the problem itself and from the use of a different platform. 

A major difference, for example, i6 the pumping system required to maintain 
the mass filter in vacuum. Due to the high ambient pressure and the long 
duration of a balloon flight, a pumping system with a long standing time 
(several hours), a large pumping capacity and a high pumping speed will be 
required for a balloon-borne instrument. 

A general view of a typical experimental configuration is shown in Figure 
1. The gondola consists of three different parts: a cryopump with the mass 

spectrometer, an electronics compartment and a mechanical support structure. In 
the configuration of Figure 1, the electronics compartment consists of a light- 
weight aluminium cylinder, sealed by a large aluminium flange, on which the 
cryopump is fixed. 

This vessel is pressurized with dry air to avoid high voltage breakdown, 
and contains several electronic modules mounted on platforms parallel to the 
sealing flange. Hermetically sealed connectors are also provided on the latter 
for interfacing to telemetry, to telecommand and to devices located outside and 
for testing the payload. All sealings are realized with special silicone 
O-rings qualified for a -55 to +200°C temperature range. The mechanical support 
structure and crash pad are made completely of metal (aluminium) to avoid 
charging up problems. Its geometry is adaptable to specific needs, such as 
flying multiple experiment payloads. 

In constructing the payload, the use of heavily degassing materials is 
avoided as much as possible, in view of contamination problems discussed further 
on. 
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Figure 1. General view of ion mass spectrometer gondola. 

A. cryopump; B. opening device; C. electronics container; 
D. mechanical support structure. 


DETAILED INSTRUMENT DESCRIPTION 
(a) High Vacuum Part 

The high vacuum part of a typical payload as developed at our institute is 
pictured in Figure 2. It consists of a high speed cryopump, a sampling hole, a 
remote opening device and a quadrupole mass filter with an ion focussing device 
and an ion detector. 

In order to prevent scattering of sampled ions by background molecules, to 
insure proper functioning of the quadrupole^mass filter and to avoid H.V. break- 
down in the ion detector, a pressure of 10 Torr or less is required inside 
the mass spectrometer housing. Since the outside pressure ranges from 1 to 100 
mbar, a basic problem of the experiment was designing a balloon-borne pump with 
high pumping speed. To overcome this difficulty, a liquid helium cryopump was 
developed, described in detail elsewhere (INGELS et al . , 1978) and, therefore, 
explained here only briefly. 

The pump body of the liquid helium cooled cryopump is made of stainless 
steel and all flanges are sealed with copper or polyimide O-rings. Apart from a 
glass fibre superinsulation, thermal insulation is achieved through two copper 
radiation shields, cooled by cold evaporating helium gas passing through an 
attached spiral. An inner blackened chevron baffle reduces radiation from the 
ion focussing lens, quadrupole mass filter and ion detector. 
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Figure 2. Detailed view of high vacuum part of the mass spectrometer. 

I. spring loaded opening device; 2. protechmcal cable cutter; 

3. liquid helium container; 4. superinsulation; 5. liquid helium 
filling port; 6. electrical feedthrough; 7. spiraltron ion detector; 
8. chevron baffles; 9. quadrupole mass filter; 10. ion lens; 

II. polyimide plunger. 


The liquid helium content of the reservoir is about 2.3 liter, for a 
standing time of about 12 hours. Measured pumping speed is about 1,200 liter/ 
sec. OLSON et al . (1978) have designed a somewhat different cryopump for a 
balloon-borne ion mass spectrometer; whereas, recently Arnold and colleagues 
(VIGGIANO et al , , 1983) have introduced the use of liquid neon pumps for this 
kind of experiment. Due to the larger heat of evaporation of liquid neon, such 
pumps have a much longer standing time and, therefore, a simpler, more light- 
weight design can be used. However, in order to pump light gases such as 
He and Ne efficiently, a charcoal absorption layer must be added to the inside 
of the cryosurf aces . 

On the downward looking input flange of the cryopump, a smaller flange is 
mounted, electrically, of insulated from the pump body. This allows the 
application of a draw-in potential to attract ambient ions towards the 
instrument. A sampling aperture is drilled in the central part of the inlet 
flange, where thickness is reduced to 0.1 nm. The diameter of the sampling hole 
is a compromise between two conflicting requirements: maximum ion signal and 

minimum inside pressure. Most of our data were obtained with a sampling hole of 
0.2 mm diameter. 

Prior to balloon launching the sampling hole is covered with a polyimide 
plunger, which can be removed by remote control as soon as the gondola has 
reached the desired working altitude. This opening system has been described in 
detail previously (INGELS et al., 1978). 
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In the configuration of Figure 2, the quadrupole mass filter and ion 
focussing device are mounted rigidly on the sampling flange to avoid centering 
problems. In previous flights, quadrupoles having a rod radius of 6.3 mm and a 
length of about 120 mm have been used (type: Finnigan and Vacuum Generators). 

OLSON et al. (1978) describe the use of a similar mass filter, although 
application of smaller quadrupoles (rod diameter = 4.5 mmm) has been reported 
(ARNOLD and QIB, 1984). 

In our setup, the ion focussing device consists of a single cylindrical 
lens element placed in front of an ionizer delivered by the quadrupole manu- 
facturer. The elements of this ionizer can be put at different potentials, 
according to whether it is used as an electron impact ionizer or as an element 
of the ion focussing lens. Switching from lens to ion source configuration can 
be done by remote control. The ionizer has proven to be very useful for in 
situ calibration of the instrument, which allowed an unambiguous ion mass 
determination (ARIJS et al., 1980). At present it is merely used for preflight 
testing. 

Although the use of the ion lens contributes to an increase in the ion 
signal, it also hampers the expansion of the gas beam, which carries ions 
through the sampling hole. As a result, a rather high pressure exists just 
behind the inlet aperture. Since in this region ions are accelerated by 
electric fields the ion abundance data can be falsified to a high degree by 
field-induced ion dissociation. 

To overcome this ARNOLD and QIU (1984) recently used a configuration where 
no focussing device was used. Ions were simply guided into the quadrupole by a 
50 V potential on the mass filter field axis (pole bias). 

After passing through the mass filter, the ions are collected by an ion 
detector. In view of the lower signal, pulse counting techniques have to be 
applied and, therefore, an appropriate detector needs to be chosen. In our case 
a continuous dynode electron multiplier (Spiraltron SEM 4219 - Galileo Optics) 
was used, although the application of other types (Johnston Laboratories MM-1) 
was also reported (OLSON et al . , 1978). 

(b) Electronics Package 

The various modules located in the pressurized part of the payload are 
shown in the block diagram of Figure 3. First there is the lens supply unit, 
furnishing the electrostatic ion lens and the sampling aperture with appropriate 
voltages of selectable polarity, referenced respectively to the sampling plat 
and to the metallic structure. It can operate either in an ion sampling or in 
an ion source mode, the latter mode only being used during in-flight mass scale 
calibration or for payload testing. 

The quadrupole supply unit produces the RF and DC voltages required to 
drive the quadrupole mass filter. It is designed so that RF and DC voltages can 
be controlled independently, allowing a great variety of operating modes. The 
mass range covered by the quadrupole depends upon the selected RF frequency, the 
maximum attainable RF voltage, the quadrupole dimensions and the desired mass 
resolution. In the current payload configuration a mass range from 0 to 330 amu 
at high resolution is realized with 1500 V peak-to-peak RF excitation at 2 MHz 
and a quadrupole rod diameter of 6.3 mm. 

The spiraltron ion detector, located behind the quadrupole mass filter, 
requires high voltages, delivered by a separate module. About 3 kV between 
input cone and anode is needed by the spiraltron for proper operation in the 
saturated pulse counting mode. Furthermore, the module allows for input cone 
voltage switching between -3 kV or +1.5 kV, in order to accommodate for positive 
or negative ion acceleration. 
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Figure 3. Block diagram of onboard electronics. 


In the pulse buffer, which is AC coupled to the anode of the detector, 
collected charge pulses are converted into voltage pulses with mean amplitudes 
of a few hundred mV. Instead of using the classical charge amplifier approach, 
a very simple circuit, based on a high speed, low input capacitance buffer 
amplifier (NS LH0033), is applied here. The charge cloud leaving the spiraltron 
is merely integrated on the RC network, formed by the total wiring capacitance 
(about 20 pf) seen by the anode and by a 1 kn leak resistor. 

Since collisions with neutrals must be avoided in the mass filter assembly, 
it must be operated at pressures below about 10 Torr. Therefore, inside 
pressure is continuously monitored by a Penning gauge in order to verify the 
cryopump's operation. Another high voltage unit supplies this gauge with about 
2.3 kV. The current consumed by it, which is related to pressure, is measured 
by a logarithmic electrometer. This value is included in the telemetry format 
of the payload. 

In order to enhance the performance of the balloon-borne mass spectrometer, 
a microprocessor (INTEL 8080) based control and data management unit has been 
included in the instrument. The design objectives and its hard- and software 
implementation have already been described before (NEVEJANS et al., 1982). The 
unit executes the following main tasks: it controls the polarity of the lens 

voltage and the spiraltron high voltage; it determines automatically the optimum 
potential of the sampling aperture; it produces the ramping voltages for the 
quadrupole supply; it counts the pulses leaving the pulse buffer and stores the 
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measured values in a spectrum memory; it acquires important parameters of the 
ambient; it does renote control interpretation and organizes transmission of 
spectrum and housekeeping data in analog and digital (PCM) formats. All of 
these duties are governed by software stored in onboard nonvolatile memory. 

One of the key features of this microprocessor based unit is the selection of an 
appropriate mass scanning mode. In our instrument, all scan-defining parameters 
can be optimized: resolution mode, extent of the scan range and number of 

spectrum channels. Two major resolution modes are implemented (NEVEJANS et al., 
1982). The constant resolution mode, wherein ions are filtered with a mass 
independent, constant, finite resolution (m/Am), can be used either to enhance 
throughput of the mass filter at reduced resolution or to access ion masses, 
beyond the normal mass range of the instrument. This mode is most often applied 
during valve controlled balloon flights, when relative abundance profiles of 
known stratospheric ion species are measured. 

An offspring of the constant resolution mode is the integral mode. Here, 
only RF is supplied to the quadrupole. This results in a shift of the rising 
edge of all mass peaks towards the low end of the mass scale. Consequently, all 
mass peaks overlap and, therefore, total ion abundance can be measured at a 
fixed low level, RF excitation. When the RF voltage is scanned however, the 
mass filter behaves as a high pass filter. 

A second mode is the constant peak-width mode. The resolution (m/Am), 
controlled by the ratio of the DC and RF supply voltage of the quadrupole mass 
filter, becomes a linear function of mass so that peak-width (Am) remains 
constant all over the mass range. The constant peak-width mode is used 
primarily for the unambiguous determination of mass numbers with a 0.1 to 0.8 
amu uncertainty. 

(c) Telemetry, Telecommand and Data Reduction 

The described assembly of electronic modules is connected to an independent 
telemetry-telecommand (TM-TC) system. During flights over Southern France a 
CNES supplied package is used (SITTEL TM-TC), which has nine analog channels 
(IRIG 5, 6, 7, 9, 11, 12, A, C and E) and a maximum of 27 latching relays (three 
different addresses with 9 relays each) available to the user. Up to three 
fully separately powered and integrated mass spectrometers can be interfaced 
with one such package without any difficulty. Each instrument has then three 
analog telemetry channels at its disposal: one for the digital serial data 

stream generated by the on-board microprocessor (PCM rate between 700 and 1300 
bits /sec) and two for transmitting analog copies of technological data and mass 
spectra. 

Each mass spectrometer utilizes eight latching relay contacts to form a 
digital telecommand input and one additional for power on-off switching. The 
digital information is used by the microprocessor software to select one of the 
stored measuring programs or to perform service tasks. 

The digital data, transmitted over an analog telemetry channel, is handled 
by a ground borne real-time computer (HP 1000 series), equipped with PCM bit and 
frame synchronizers (EMR models 720 and 727) and running multitasking software 
(HP RTE-M III operating system). Housekeeping data about the ambient and 
internal conditions (pressure, temperatures, voltages, vacuum, etc.) of the 
payload are displayed and continuously refreshed on a console, while spectrum 
data are stored on digital tape. The latter data can also be plotted on the 
console or on a printing plotter. 

CALIBRATION OF THE INSTRUMENT 

To calibrate the mass scale of our instrument, to check its sensitivity and 
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to determine the appropriate parameters, which are stored in the memory of the 
microprocessor control unit, extensive laboratory studies have been performed. 

The stratospheric plasma is simulated in a vessel, fitting on the mass 
spectrometer, by ionizing a nitrogen - nitric oxide gas mixture (typically 1% NO 
in N^) with UV irradiation (the 129.6 nm line of a xenon discharge). The 
pressure of this ga6 mixture can be stabilized over a wide range (0.1 to 100 
Torr) and the ratio NO/^ as well as the xenon lamp intensity are kept 
constant simultaneously. 

In this way a reproducible ion population is created having a number 
density similar to the stratospheric situation, but with a different 
composition. 

The same setup has been used previously (ARIJS et al., 1982a) to investi- 
gate the effect of the ion lens field on proton hydrate cluster breakup. 

It should be emphasized that so far only relative abundance ratios of 
different ion species were determined and no absolute ion densities have been 
deduced from stratospheric ion composition measurements. The derivation of 
absolute number densities would require a special calibration system, which has 
not been realized yet, due to the difficulty of independently measuring ion 
densities. 

Another remaining problem is mass discrimination. Since, in most of the 
present applications of ion composition measurement, moderate resolution spectra 
have been used, these discrimination effects are believed to be of minor 
importance. If in the future, however, high resolution spectra will be used, 
this phenomenon has to be analyzed in the laboratory. 

MEASUREMENTS 

(a) State of the Art Review 

Since 1977 many measurements have been performed with balloon-borne ion 
mass spectrometer, mainly by two groups (Max Planck Institut fur Kernphysik, 
Heidelberg and Belgian Institute of Space Aeronomy, Brussels). The results of 
most of these measurements as well as typical ion spectra have been reported in 
the literature (see reviews of ARNOLD (1980) and ARIJS (1983) and references 
cited therein). Therefore, we will only give a concise description of the state 
of the art and focus more on future advances in the measurement technique, and 
on problems related to possible misinterpretation of data. 

The major positive ions defected in the altitude region 20 to 45 km turned 
otjt to be proton hydrates - H,0 (H»0) - and nonproton hydrates of the 
H Xj, (H,0 ) , where X most probain.y n is acetonitrile (CH^CN) (ARNOLD, 1980; 
ARIJS, 1983). From the relative ion abundance measurements; the mixing ratio of 
CH_CN was derived in the altitude region 45 to 20 km (HENSCREN and ARNOLD, 

1981a; ARIJS et al . , 1983a). Although recently an in situ formation mechanism 
for CH^CN has been proposed (MURAD et al., 1984), the CH-CN mixing ratio 
profiles as deduced from ion mass spectra seen to suggest a release of 
acetonitrile at the Earth's surface, followed by upward diffusion and 
destruction by OH in the stratosphere (BRASSEUR et al., 1983). 

Apart from the major positive ions, many other positive ions have been 
detected in the stratosphere (HENSCHEN and ARNOLD, 1981b; ARIJS et al., 1982a). 
It was shown that the measuronents of minor mass peaks offer the possibility of 
trace gas detection of species such as CH^OH, NELj, a.s.o. It should be 
noted, however, that signal instabilities, which are very pronounced at low ion 
count rates and possible contamination effects, to be discussed hereafter, 
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hamper Che usefulness of the method as a full grown analytical tool. 


The major negative ions observed in the stratosphere are 
NO, (HNO,) and HSO, (BLSO, )„ (HN0-) cluster ions (ARNOLD, 
19&); Ar2j§, 1983). 4 2 4 * 3 m 


From the relative ion abundances, both HNO^ and H^SO^ concentration 
were derived. The derivations of nitric acid vapour mixing ratios, however, are 
strongly influenced by cluster breakup. The derivations of sulfuric acid ion 
concentrations, on the other hand, are not subject to these problems and are of 
great importance for our understanding of the sulfur chemistry and aerosol 
formation (ARIJS et al., 1981, 1982b, 1983b, d; ARNOLD et al., 1981a, 1982; 

ARNOLD and BUHRKE, 1983; VIGGIANO and ARNOLD, 1981, 1983; QID and ARNOLD, 1984). 
More laboratory data concerning the appropriate ion-molecule chemistry, however, 
are needed to exploit the data to a full extent. 

The detection of minor mass peaks in negative stratospheric ion mass 
spectra (ARIJS et al., 1982b; MCCRUMB and ARNOLD, 1981) has also opened a way 
for detecting trace gases so far not measured (such as HC1, H0C1, a.s.o.), but 
again the same remarks as for positive ions concerning signal instabilities and 
contamination are valid. 


A problem encountered in all in situ analytical techniques and strongly 
influencing stratospheric ion mass spectrometry data is contamination. It is 
evident that a technique capable of detecting trace species in the ppt range and 
lower is extremely sensitive to this problem. It has been established that 
contamination is mostly influencing the data during measurements obtained in the 
ascent phase of the balloon flight or at float altitude (VIGGIANO and ARNOLD, 
1983; ARIJS et al., 1983a, c). Since contaminants originate from outgassing of 
the balloon and gondola, the problem is also more acute in daytime measurements. 
An illustration of the effect of degassing is given in Figure 4, which shows two 
positive ion spectra obtained by our group at 45.7 km altitude during daytime. 
Both spectra were recorded at float altitude when the balloon was oscillating 
between 45.6 and 45.8 km. Spectrum A, containing mass peaks belonging mainly to 
the natural ions, was obtained during a descending phase of the altitude 
oscillation, whereas spectrum B was taken during an ascending motion of the 
payload. It is clear that in spectrum B, contaminants falsify the results con- 
siderably. This is explained by the fact that in spectrum B the gondola moves 
through the wake of the balloon, containing gases released by the balloon and 
some payload construction materials. These gases are mainly organic vapours, 
with a high proton affinity and thus are participating in the positive ion 
chemistry (ARIJS et al., 1983a). 

It is also evident that these effects even if not as dramatic as those in 
Figure 4, may cause misinterpretations, especially in the detection of minor 
mass peaks. 


To overcome these contamination problems, most of the data of our group 
used so far were obtained at nighttime or during descent in flights with valve 
controlled balloons. Another remedy which might be envisaged in the future is 
the use of reel down mechanisms, such as the one recently developed by ANDERSON 
and HAZEM (1983). 

One of the most appealing applications of stratospheric ion mass spectro- 
metry is the detection of trace gases. At present this has been realized by two 
methods, described to full extent by ARNOLD et al. (1980), namely the steady- 
state method and the equilibrium method. In the latter one, the partial 
pressure p(B) of a trace gas B is derived from 
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Figure 4. Two typical positive ion spectra obtained around 
45 km altitude. Spectrum A is recorded during a descending 
motion of the balloon, spectTum B during an ascending motion. 


p(B) 


[A 1 B 1 
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[A* B ] 

1 n J 
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where square brackets denote number densities and K is the equilibrium constant 
in atm . The core ion ^ can be positive or negative; major core ions 
observed in the stratosphere are H^O , H + CN^CN, NO^ and HSO^ . 
The main ligands B detected are H.0 , CH^CN for positive ions and HNO, and 
HjS^i for negative_ clusters. For some cluster families, such as 
H l° * H 2 0) n and N0 3 -<HNO ),, equilibrium constants are known from 
laboratory measurements and thus p(B) can be derived from the ratio of the 
fractional ion abundances. 


In general, however, the ligands B are weakly bound to the core ions. 
Immediately behind the sampling aperture of the mass spectrometer, where ions 
are accelerated by the focussing lens or by the field axis potential of the 
quadrupole and where pressure is still rather high, the charged particles make 
many collisions with neutral gas molecules. As a result, field-induced 
collisional dissociation of a fraction of the ions occurs and some of the 
ligands B are "boiled orf." Therefore, the measured fractional ion abundance 
ratio can differ substantially from reality. 

A typical example of such measurements is shown in Figure 5, where we have 
plotted the ratio H-l 0 + (J, 0 ) 3 /H-0 + as measured in a balloon flight 
performed with a 1 ,000,01)0 m 3 Winzen balloon on 23 September 1982 over 
Southern France. 
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ABUNDANCE RATIO OF H 3 0*(Hj 0) 3 / HjO'fHjOj 


Figure 5. Relative abundance ratio of signal at 
mass 73 to signal at mass 55 as measured m 
a high altitude balloon flight. 


For comparison, we have also shown the ion abundance derived by 
[H 3 0 + (H 2 0) 3 ]/lH 3 0 + (H 2 0) 2 ] - KlB^O] 


and 


- RT In K = AH - TAS 

where the values of LAU et al. (1982) were used for AH and AS and a constant 
water vapour mixing ratio of 3 ppm was assumed throughout the whole altitude 
region considered here. 

The shaded area on Figure 5 represents the results of a calculation for 
standard atmosphere temperature conditions (lower limit for summer and upper 
limit for spring-fall). Use of the in situ measured temperature profile leads 
to the horizontal bars, the width of which is proportional to the temperature 
measurement error. As can be seen, the measured ion abundance ratio below 45 km 
is + much lower thjm expected, which is due to the dissociation of 
H3O ^ nto H 3° ^ H 2°^2 * n the i nstrument * 

Above 45 km the measured ion abundance ratio is larger than expected. 

These measurements were taken at float altitude during daytime; whereas the data 
below 45 km were recorded during descent, due to the cooling of the balloon 
after sunset. 

The too-high ion abundance ratio above 45 km is ascribed to contamination, 
resulting from balloon and payload out gas sing caused by solar heating. One of 
the major products of this desorption is water, which causes a shift of the 
proton hydrate distribution towards higher masses. The high signal at mass 73 
H 3 o + 040)3), responsible for the too high ion abundance ratio, partly 
results from cluster breakup of mass 91 (^0 (HjO)^). 
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Figure 5 is at the same time a good illustration of contamination and 
cluster breakup. It should be pointed out that the effect of cluster breakup 
is strongly dependent on ambient pressure (altitude) and instrument configura- 
tion (diameter of sampling hole, focussing potentials a.s.o.) and may, there- 
fore, vary from flight to flight. For the measurements of Figure 5, a sampling 
hole of 0.4 mm diameter was used. 

It is clear that in order to use ion mass spectrometric data for analytical 
applications, the cluster breakup effects should be avoided or very well under- 
stood so that they can be corrected (ARNOLD et al., 1981b; ARIJS et al., 1982a). 

In any case, care should be taken in the interpretation of fractional 
abundance data, especially if they are used for deriving unknown thermochemical 
values (AH and AS) of some clustering reactions. 

It should be pointed out that the derivation of CH^CN number densities 
from positive ion data and of H^SO^ vapour concentrations from negative ion 
measurements are believed not to be influenced by cluster breakup. In these 
cases, the fractional ion abundances of two distinct ion families (proton 
hydrates and nonproton hydrates; NO, - and HSO^ ions) are compared. 

Cluster breakup only results in a modification of the distribution within the 
same ion family and is unlikely to cause the conversion from one family to 
another. 

CONCLUSIVE REMARKS 

During the past decade considerable progress has been made in the 
technology of detection and identification of ambient ions in the stratosphere. 
The measurements performed hitherto have provided an insight in the major 
processes governing the ion composition in the altitude region 20 to 45 km. To 
get a more detailed understanding and to exploit the data to a further extent, 
however, some developments will be necessary. More extensive calibrations of 
existing instruments are necessary to be able to convert ion count rates to ion 
number densities, to correct for cluster breakup and to understand possible 
mass discrimination effects. New instruments or improvements to existing ones 
need to be developed to avoid cluster breakup and to increase the sensitivity 
of mass spectrometers. Since the latter is at lower altitudes mainly limited by 
ion scattering by residual gas molecules in the instrument, more powerful 
pumping systems will be part of the new developments. On the other hand, the 
use of other types of mass filters, such as magnetic mass sector instruments and 
simultaneous detection techniques may open new perspectives. The altitude range 
of existing measurements need to be extended and the latitudinal variation of 
ion compositions should be investigated. To extent the altitude range the use 
of other measurement platforms such as parachute borne rocket payloads, airplane 
borne instruments (HEITMANN and ARNOLD, 1983) and the use of reel down 
mechanisms may become necessary. Finally, in order to interpret the data, more 
laboratory work on the appropriate ion-molecule reactions is urgently needed. 
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ABSTRACT 

The history of measurements of the stratospheric aerosol layer is briefly 
reviewed. Balloon-borne measurements played a major role both in the discovery 
of the layer and in the subsequent study of its nature. Following this 
introduction, stratospheric aerosol sources and balloon-borne measurement 
techniques are reviewed. Finally, the current status of this research area is 
outlined with data obtained mainly by balloon-borne techniques. It is con- 
cluded that a great deal has been learned m the past five years, due mainly to 
the increase of the occurrence of major volcanic eruptions, which are the 
predominant source of stratospheric particles, but also due to the development 
of instrumentation and techniques which have seen considerable application on 
the balloon platform. 

INTRODUCTION AND BACKGROUND 

While knowledge of suspensions of small particles in the upper atmosphere 
through twilight observations dates back to at least the 18th century, the study 
of the stratospheric aerosol layer did not become a science until Junge's 
pioneering measurements in 1959 (JUNGE et al.» 1961). Using balloon-borne 
particle impactors, they determined that a stratospheric relative maximum in the 
particle concentration existed at about 20 km altitude for particles with radii 
of about 0.1 pm. They also measured condensation nuclei (r 'v 0.01 pm) using 
expansion cloud chambers and found them to be predominantly a tropospheric 
phenomenon with only minor concentrations in the stratosphere. 

ROSEN (1964) developed a balloon-borne optical individual particle counter 
which was a forerunner of an instrument now known as a dustsonde used 
extensively over the past twenty years in balloon applications all over the 
world. About the same time, laser radar (lidar) was applied to stratospheric 
aerosol studies (FIOCCO and GRAMS, 1964). In the seventies, balloon-borne 
particle impactors (BIGG et al., 1970; BIGG, 1975, 1976; BROWNLEE et al., 1976; 
GRAS and LABY, 1978, 1981) and filter samples (LAZRUS et al., 1971; LAZRUS and 
GANDKUD, 1974, 1977; CASTLEMAN et al., 1974) were used to study the nature of 
the aerosol. 

Extensive volcanic activity during the 1980-1982 period resulted in a 
renewed interest and accelerated research on stratospheric particulates by 
numerous techniques. It is now clear that the cpnc_^3tration of optically active 
(r > O.lp m) aerosol varies from the order of 1 cm during volcanically 
quiescent periods to at least ten times this value following large volcanic 
eruptions such as that of Agung in 1963 and El Chichon in 1982. Condensation 
nuclei (r ^ 0.01 pm) typically have tropospheric concentrations of the order of 
10“ cm but fall off rapidly in the stratosphere and become part of the 
normal stratospheric aerosol size distribution. 

Under nonvolcanic conditions, the stratospheric aerosol size distribution 
may be typified by a log normal distribution with a m^an radius of about 
0.08 pm, a total number concentration of about 10 cm and a o of about 1.8 
at an altitude of about 20 km (PINNICK et al., 1976; HOFMANN and ROSEN, 1981; 
ROSEN and HOFMANN, 1977). The global mass of stratospheric aerosol is estimated 
to be about 0.5 Tg (1 Tg = 10“ metric tons) under normal conditions and as 
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large as 10 Tg following major volcanic eruptions (HOFMANN and ROSEN, 1983a; 
MCCORMICK and SWISSLER, 1983). The quiescent layer appears to be at a slightly 
higher altitude in equatorial regions than at the poles but is clearly a 
continuous global phenomenon (ROSEN et al., 1975). 

While stratospheric aerosol contains many trace elements, the overwhelming 
component is sulfate (LAZRUS et al., 1971). The chemical structure appears to 
be sulfuric acid with an approximate 75% H-SO, , 25% H^O composition 
(ROSEN, 1971). Recent measurements suggest that the acid fraction probably 
varies from 60-80% depending on tanperature and water vapor concentration 
(HOFMANN and ROSEN, 1983a). That the aerosol exists as a liquid droplet ha6 
been suggested from its appearance following impaction and also from the 
generally good agreement obtained in comparing in situ measurements with lidar 
backscatter measurements (NORTHAM et al., 1974; PINNICK et al., 1976; SWISSLER 
et al., 1982; HOFMANN et al., 1983), a comparison which relies on the particles 
being spherical droplets of aqueous sulfuric acid. 

STRATOSPHERIC AEROSOL SOURCES 

While >98% of the aerosol which persists in the stratosphere is composed of 
sulfuric acid under both background and volcanically disturbed conditions, the 
precursor sulfurous gas appears to be quite different in the two cases. Most of 
the volcanic gas is SOj with possibly some H 2 S which is rapidly converted to 
SOj. Thus, following volcanic eruptions, H2S0^ vapor is formed through 
the chemical oxidation of S0„. Recent studies suggest that the mo6t plausible 
chain of reactions involves the hydroxyl radical OH but in a catalytic manner so 
that OH is conserved in the process (MCKEEN et al., 1983): 

S0 2 + OH + M -► H0S0 2 + M 
hoso 2 + 0, -»• SO, + ho 2 

HO + NO OH + N0 2 

SOj + H 2 0 HjSO^ /_ 

net: S0 2 + 0 2 + NO + HjO -*■ HjSO^ + N0 2 

The lifetime for this process in the stratosphere (25 km) is estimated to be 
about 40 days. SO 2 decay and aerosol growth observations at 25 km following 
the El Chichon eruption are m agresnent with this estimate (HEATH et al., 1983; 
HOFMANN and ROSEN, 1983b). 

The lifetime of SO 2 in the troposphere is considerably shorter owing to 
its water solubility and rapid washout resulting in the quiescent SO 2 level in 
the lower stratosphere being only about 10-100 parts per trillion by volume 
(JAESCHKE et al., 1976; GEORGII and MEIXNER, 1980; INN et al., 1981). Thus, 
nonimpulsive anthropogenic emissions of S0 2 at the surface, which by far 
exceed eruptive volcanic emissions, are not effective as stratospheric aerosol 
precursor gas sources. It is believed that the surface sulfur source for the 
background stratospheric aerosol is naturally produced and anthropogenic 
carbonyl sulfide (OCS) which has a relatively long lifetime (''-500 days) in the 
troposphere, diffuses into the stratosphere where it is photo-dissociated, 
eventually forming S0 2 (TURCO et al., 1979). 

Most of the very large particles (greater than 1 pm radius) in the strato- 
sphere appear to be composed of aluminum oxide (BROWNLEE et al., 1976) which are 
apparently deposited there through the exhaust of solid fuel rockets. Due to 
their small number (approximately 1 m above a radius of 1 pm before 1976), 
they do not contribute substantially to the total mass distribution but warrant 
continued observations as the use of solid fuel propellants increases. 
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BALLOON-BORNE MEASUREMENT TECHNIQUES 

Table 1 summarizes and references techniques which have been used to study 
the stratospheric aerosol and its precursor gas (H^SO^) from the balloon 
platform. While the reference list may not be exhaustive, the table includes 
all techniques which have contributed to the understanding of this phenonmenon. 
In the following we will briefly identify the advantages and disadvantages of 
each. 

(a) Impactors 

This technique involves the physical impaction of the particles (or 
droplets) on a surface which is then subjected to analysis by optical or 
electron microscopes or by spectroscopic techniques. Advantages include size 
and shape discrimination and some compositional information. Disadvantages 
include possibility of particle loss or contamination prior to analysis, poorly 
known efficiency for the smaller particles and manpower intensive analysis 
requirements . 

(b) Filters 

This method of sampling involves the trapping of the aerosol in a filter 
matrix with subsequent removal and analysis using standard chemical techniques. 

Table 1. Balloon-borne aerosol measurement techniques. 


METHOD 

TYPE 

REFERENCES 

Impactors 

Cascade/Optical Microscope 
Electron Microscope 

JUNGE et al., 1961 
BIGG et al. , 1970 
BIGG, 1975 

BROWNLEE et al., 1976 
GRAS and LABY, 1978 

Filters 

IPC Cellulose 

LAZRUS et al., 1971 
LAZRUS and GANDRUD, 
1974, 1977 

SEDLACEK et al., 1983 

Optical 

Particle 

Counters 

White Light 

ROSEN, 1964 

GRAS and LABY, 1981 

HOFMANN and ROSEN, 1982 


Laser 

MIRANDA and FENN, 1974 

Condensation 

Nuclei 

Counters 

Expansion Chambers 

Growth Chamber /Optical 
Particle Counter 

JUNGE et al., 1961 
KASELAU et al., 1974 

ROSEN and HOFMANN, 1977 

Ion Mass 
Spectrometers 


ARNOLD and FABIAN, 1980 
ARIJS et al., 1981 
ARNOLD and BUHRKE, 1983 

Remote 

Sensing 

Solar Occultation 

PEPIN, 1969, 1970 
THOMPSON, 1981 
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Advantages include accurate chemical composition determinations and inherent 
simplicity in obtaining a sample. Disadvantages include lack of size distribu- 
tion information and possibility of contamination. 

(c) Optical Particle Counters 

These instruments pass an air sample through a beam of light and measure 
the amount of light scattered by individual particles or droplets in the sample. 
Advantages include large amounts of data easily collected, relatively good size 
specificity given the aerosol composition and shape (e.g., for sulfuric acid 
droplets) and simple determination of volatile or nonvolatile state of aerosols. 
Disadvantages include necessity for accurate calibration, lack of compositional 
information other than degree of volatility, degradation of size information for 
unknown distributions of nonspherical particles and a lower size limit of about 
0.1 pm due to the inability of discriminating light scattered by smaller 
particles from that due to air molecules. 

(d) Condensation Nuclei Counters 

While these instruments also generally rely on a particle's ability to 
scatter light, they are distinguished from an optical particle counter because 
they incorporate a growth chamber wherein particles normally too small 
(<0.1pm) to be detected by light scattering, are caused to grow by passing them 
through a region of supersaturation of a working fluid (water or ethylene 
glycol). Super saturation is either caused by expansion or by a thermal gradient 
or transition. Following growth to optically active sizes, the particles are 
counted either individually by light scattering (using photography or optical 
particle counters) or collectively by the amount of light they absorb. 

Advantages include the small sizes (''-0.01 Pm) which may, in principle, be 
detected, the possibility of varying effective size by varying the degree of 
supersaturation and the ability to handle high concentrations. Disadvantages 
include poor overall size discrimination, lack of composition discrimination and 
poor resolution of low particle concentrations for instruments that use photo- 
graphic detection or collective light absorption. In addition, any instrument 
which attempts to measure particles smaller than 0.1 m, must contend with loss 
of particles due to preferential diffusion of small particles to the walls of 
the instrument. 

(e) Ion Mass Spectrometers 

These devices are capable of obtaining the mass spectrum of charged mole- 
cular clusters and thus address the measurement of aerosols in their earliest 
stage of growth. Advantages include the ability to obtain the vapor 

concentration and the relative mass distribution for cluster masses less than 
several hundred amu. Disadvantages include the inability of obtaining the ion 
density so that concentrations of clusters cannot be determined and the 
susceptibility to ion fractionation at the inlet to the instrument. 

(f) Remote Sensing 

This technique involves the measurement of the extinction (absorption and 
scattering) of sunlight by the stratospheric aerosol layer through a tangent 
atmospheric pathlength during sunrise or sunset. Advantages include a 
relatively large signal and aerosol/molecular resolution through wavelength 
discrimination. Disadvantages include the reliance on an optical model of the 
aerosol for data interpretation, poor horizontal resolution due to the long 
atmospheric pathlength (^200 km) and a need for high pointing accuracy to obtain 
desired vertical resolution (<1 km). 
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CURRENT STATUS 

The last five years (1978-1983) have seen a revolution in the study of the 
distribution of particulate matter in the upper atmosphere. Most of what is now 
koown about the sources and detailed formation processes of the quiescent and 
disturbed stratospheric aerosol layer was pieced together through observation 
and theory during this period. The first two years of this period (1978-1979) 
were characterized by minimal major volcanic activity and consequent low back- 
ground conditions for the stratospheric aerosol layer. Thus HOFMANN and ROSEN 
(1981), through numerous balloon soundings using the.dustsonde, were able to 
identify a background aerosol level of about 0.5 cm for radii greater than 
0.15 pm at about 20 km. Comparisons with Junge's measurements in 1959, also a 
period of minimal volcanic activity, suggested an increase of about 9% per year. 
Total sulfate measurements by SEDLACEK et al. (1983) between 1973 and 1979 (two 
periods thought to be relatively free of major volcanic eruptions) suggested a 
6.2% per year increase and when compared with Junge's measurements, an 8% per 
year increase. Such an increase rate is not surprising in view of the increase 
in the burning of fossil fuel and the concomitant release of OCS, a known source 
gas for the quiescent stratospheric aerosol layer. 

The period 1980-1982 saw the eruption of four volcanos which caused major 
stratospheric perturbations. These include Mt. St. Helens (46.2° N) in the 
state of Washington on May 18, 1980, Alaid (50.8°N) in the Kurile Islands on 
April 28, 1981, an apparently unobserved and thus unidentified eruption in early 
January, 1982, and the series of eruptions of El Chichon (17.3°N) in Mexico on 
March 28, April 3 and April 4, 1982. These eruptions, due either to their 
accessible location or sheer magnitude, provided the opportunity to study the 
process of gas to particle conversion in the stratosphere. The nucleation of 
new sulfuric acid droplets from the vapor phase and their subsequent growth to 
optically active sizes were observed and studied carefully. 

In the following we will summarize the current status of stratospheric 
aerosol measurements. A great majority of this information has been gathered 
through the use of balloon-borne sensors. Indeed, the height of the El Chichon 
volcanic injection (25-30 km) dictated the use of balloons for Jn situ 
measurements as currently available research aircraft have ceiling altitudes in 
the vicinity of 21 km. 

(a) Precursor Gases 

Figure 1, from ARNOLD and BUHRKE (1983), shows H-SO^ and HSOj 
measurements made by balloon-borne ion-mass spectrometers in southern France 
(44°N) in 1982 and include the possible identification of a sulfuric acid layer 
at about 25 km, on June 10, 1982, believed to be due to the eruption of El 
Chichon. No in situ measurements of the sulfuric acid vapor precursor, S0„, 
exist above 21 km as no reliable balloon-borne technique has been developed ana 
currently available cryogenic sampling techniques used on aircraft (INN et al., 
1981; VEDDER et al., 1983) have not been used on balloons. However, ground- 
based and aircraft remote sensing of the S0_ overburden indicate large 
enhancements following volcanic eruptions (EVANS and KERR, 1983). 

(b) Condensation Nuclei 

Figure 2, from HOFMANN and ROSEN (1983c), shows profiles of condensation 
nuclei (r 2.0*01 ym) at Laramie (41 °N) before and after the large increase at 
high altitude observed on January 28, 1983. These data were obtained with a 
balloon-borne condensation nuclei counter employing a growth chamber and optical 
particle counter. Typically, one sees high variable concentrations in the 
troposphere (below about 12 km in Figure 2), low stable concentrations in the 20 
km region and the annual (wintertime) increases at about 30 km. The latter have 
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Figure 1. Concentration profiles of HjSO^ and HSO, 

vapor obtained by ion-mass spectrometers on balloons in 
southern France (from ARNOLD and BUHRKE, 1983). 



Figure 2. Profiles of condensation nuclei (r>0.01 pm) 
obtained with balloons at Laramie, Wyoming before 
(December 1982) and after the large increase at high 
altitude that began on January 28, 1983 (from HOFMANN and 
ROSEN, 1983c). 
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been observed every year since 1979 and are believed to be due to nucleation of 
new very small sulfuric acid droplets during stratospheric warming episodes in 
the polar region and subsequent transport to lower latitude during the incursion 
of arctic air at high altitude. It is possible that the high temperature 
associated with the warmings cause sulfuric acid aerosol to evaporate creating 
vapor which then recondenses due to cooling during transport. 

Figure 3 shows the condensation nuclei mixing ratio versus time at Laramie 
in the region of the high altitude events (25-35 km) and in the stable 20 km 
region. The large increase at 20 km in 1982 was due to the nucleation of new 
aerosol following the eruption of El Chi chon (HOFMANN and ROSEN, 1983a). The 
general increase in the background level at 25-35 km since 1980 is believed to 
be due to the increase in volcanic activity since this time, the latter 
supplying enhanced levels of both sulfuric acid aerosol and vapor which are 
transported to the polar regions during winter. Instruments launched by the 
University of Wyoming group in October 1983, on the largest balloon ever flown 
in Antarctica, revealed the presence of the El Chichon aerosol over McMurdo 
Station (78°S) and also indicated the presence of enhanced condensation nuclei 
at high altitude, presumably due to the Southern Hemisphere analog of the events 
observed in the Northern Hemisphere since 1979. 



Figure 3. Condensation nuclei (r _> 0.01 V m) mixing 

ratios versus time in the region of the annual increases 
(25-35 km) and at 20 km as measured with balloons at 
Laramie, Wyoming. Times of major volcanic eruptions which 
perturbed the larger aerosol concentrations are indicated. 
An unidentified eruption in early 1982 is indicated by a 
question mark (from HOFMANN and ROSEN, 1983c). 
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(c) Total Sulfate 

Figure 4, from LAZRUS and GANDRUD (1977), shows the average total sulfate 
ion mixing ratio and volume concentration profiles as determined from filter 
samples on balloon flights at 33 °N during 1976. These measurements indicate 
that measurable amounts of sulfate are present to altitudes as high as 37 km. 

The advantage of balloon-borne measuranents here are obvious as aircraft 
altitudes bring one barely beyond the maximum in the vertical distribution. 
Figure 5, from MROZ et al. (1983), shows the large increase in sulfate following 
the eruption of El Chichon in April 1982. These data were obtained from 
balloon-borne filter samples also at 33°N. Again, aircraft samples were 
ineffective in obtaining data on this eruption during the first six months since 
nearly all the enhanced sulfate remained above 21 km during this period. 

SULFATE CONCENTRATION l/igm' 3 ) 



Figure 4. Average sulfate mixing ratio and concentration 
profiles measured by balloon-borne filter samplers at 
33 °N from four seasonal flights in 1976 (from LAZRUS 
and GANDRUD, 1977). 


(d) Optically Active Aerosol 

Figure 6, from HOFMANN and ROSEN (1983c), shows the aerosol mixing ratio of 
two size ranges versus time as measured at Laramie and in southern Texas with 
balloon-borne dustsondes. Times of major volcanic eruptions which perturbed the 
stratospheric aerosol at 41 °N are indicated in the figure. The El Chichon 
eruption stands out as the most significant both because of the absolute 
magnitude of the increase and the large size of the aerosol which may be 
inferred from the relative magnitudes of the increase in the two size ranges in 
Figure 6. 

Figure 7 shows aerosol profiles in six integral size ranges obtained from a 
complement of three balloon-borne particle counters at Laramie following the El 
Chichon eruption. Data such as in Figure 7 may be used to obtain the average 
radius of the aerosol which is plotted versus time after the El Chichon eruption 
in Figure 8. Two layers developed following this eruption with larger particles 
forming at higher altitude. The gradual merging of the two curves is due to 
gravitational settling of the particles in the upper layer thus decreasing their 
average size while increasing the average size in the lower layer. The particle 
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Figure 5. Sulfate concentration profiles obtained by 
balloon-borne filter samplers at 33 °N before and 
after April 4, 1982 eruption of El Chichon (from) 
MROZ et al., 1983). 



Figure 6. Maximum stratospheric aerosol mixing ratio for two size ranges 
measured by balloon-borne dustsondes as a function of time at Laramie, 
Wyoming and in southern Texas. Times of volcanic eruptions which caused 
stratospheric aerosol perturbations at Laramie are indicated (from HOFMANN 
and ROSEN, 1983c). 
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Figure 7. Aerosol concentration profiles for six size ranges as measured by 
three balloon-borne particle counters at Laramie, Wyoming on January 28, 
1983. The arrow marks the position of the tropopause (from HOFMANN and 
ROSEN, 1983c). 


size in the lower layer is typical of previous volcanic eruptions while that in 
the upper layer has never been observed before. These enhanced sizes are 
believed to be due to the large amount of sulfuric acid vapor which developed in 
the stratosphere following this eruption. The aerosol ^growth curves_^n Figure 8 
are consistent with an H ? S0, concentration of about 10 molecules cm 
and a lifetime of 25-50 days (HOFMANN and ROSEN, 1983b). 

By heating the inlet to a particle counter during a slow balloon descent in 
southern Texas in October 1982, HOFMANN and ROSEN (1983a) were able to determine 
the vapor pressure curve of the aerosol. The results are given in Figure 9. 
These data suggest that the aerosol from the El Chichon eruption was composed of 
an 80% H-SO^, 20% H2O solution in the upper layer at about 25 km and about 
60-65% HjSO^ in the lower layer at about 18 km. More than 98% of the 
aerosol was volatile at 130°C. The distribution of the minute remaining non- 
volatile fraction could be resolved and indicated that it peaked at about 16 km, 
just above the tropopause indicating again that the silicate fraction of a 
volcanic eruption leaves the stratosphere very rapidly and is not important in 
determining any long term climate effect. 

By constructing size distributions from data such as shown in Figure 7, and 
assuming a composition as suggested in Figure 9, the column mass above a given 
altitude may be determined (HOFMANN and ROSEN, 1983b). These data are shown 
plotted versus time in Figure 10 for balloon soundings at Laramie and in 
Southern Texas following the eruption of El Chichon. The gradual increase . 
observed at Laramie, which peaked in early 1983, is due to delayed meridional 
transport at 25 km which kept the upper layer of larger particles south of about 
30°N until the fall of 1982. The total mass scale on the right in Figure 10 
assumes longitudinal homogeneity and latitudinal extents as indicated on the 
mass scales. It may be used to estimate the global mass if one assumes that the 
measurements are representative of average conditions after approximate global 
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Figure 8. Aerosol mean radius versus time after the April 4, 1982 eruption of 
El Chichon in two altitude regions as obtained from balloon-borne particle 
counter data at Laramie, Wyoming and in southern Texas (from HOFMANN and 
ROSEN, 1983b). 


equilibrium has been reached. While the latter condition probably never really 
prevails, the data are consistent with at least 10 Tg present in the strato- 
sphere in early 1983. 

(e) Large Particles 

Evidence for growth of sulfuric acid aeroBol to micrometer sizes was 
obtained by balloon-borne sensors after the eruptions of Mt. St. Helens and 
Alaid. These large droplets were found to be important in calculating expected 
lidar backscattering and in the total aerosol mass (HOFMANN and ROSEN, 1982; 
HOFMANN et al., 1983). Following the eruption of El Chichon, a bi-modal size 
distribution with a mode near 1 ym radius was a dominant feature of the strato- 
spheric aero so 1 (HOFMANN and ROSFN, 1983a). These large droplets are thought to 
form, following major sulfur-rich eruptions, through growth of the preexisting 
aerosol which typically has a radius of about 0.1 ym (HOFMANN and ROSEN, 

1983d). 

Impactor samples collected by balloons have revealed the existence of 
aluminum oxide particles in the 3 to 8 ym diameter range during the period 
1970-1976 (BROWNLEE et al., 1976). The total mass of these particles appears 
to be consistent with the burning of solid fuel by rockets. The latter use 
powdered aluminum as an additive and tests have shown that a byproduct of the 
burnt fuel is aluminum oxide spheres. This manmade aerosol dominates the 3-8 
ym size range while sulfates dominate the submicrometer size range. Above 10 
ym, the particles are mostly extraterrestrial in origin with a composition 
similar to that of primitive meteorites (BROWNLEE et al., 1976). 

SUMMARY 

While this review does not claim to be exhaustive as far as stratospheric 
aerosol measuranents are concerned, it does indicate the considerable contribu- 
tion made by balloon-borne sensors. This was especially true following the 
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Figure 9. Aerosol boiling point temperature ranges (horizontal 
bars) as a function of altitude as obtained from a heated 
intake particle counter on a slowly descending balloon in 
southern Texas. The smooth curves are vapor pressure curves 
for various weight percentages of sulfuric acid and water. 

The ambient temperature, as measured during the balloon flight, 
is given and the theoretical sulfuric acid weight percentages, 
as dictated by the temperatures and by several representative 
ambient water vapor mixing ratios, are indicated by dashed 
lines (from HOFMANN and ROSEN, 1983a). 


great eruption of El Chichon in 1982 when sulfuric acid aerosol formed at 
altitudes accessible only to balloon-borne sensors. In addition, the ability to 
obtain accurate vertical distributions is invaluable in assessing aerosol 
formation processes. Although the cost of large balloons is considerable, 
instruments such as the dustsonde, which weighs less than 20 kg, can be carried 
to 30 km by relatively inexpensive balloons which are easily launched under most 
conditions, thus facilitating numerous measurements at varied locations. It is 
clear that the balloon platform will see continued use in the study of the 
aerosol distribution and its gaseous precursors in the stratosphere. 
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DAYS AFTER ERUPTION 


Figure 10. Stratospheric aerosol mass loadings above several altitudes as 
determined from balloon-borne particle counter measurements at Laramie, 
Wyoming and in southern Texas. The total mass scale on the right assumes 
longitudinal homogeneity and latitudinal extents as indicated (from HOFMANN 
and ROSEN, 1983b). 
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